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I .  INTRODUCTION 


The  objective  of  the  investigation  was  to  utilize  the  SANX  structural 
finite  element  computer  program  to  generate  the  input  data  for  ADINA  and 
NASTRAN  computer  codes.  The  SANX  program^  is  a  structural  finite  element 
code  which  has  been  developed  for  the  purpose  of  approximate  three-dimensional 
analysis  of  non-axisymmetric  configurations  in  cylindrical  coordinates.  The 
ADINA  and  the  NASTRAN3  codes  are  general  finite  element  codes  capable  of 
handling  structural  analysis  of  multi-dimensional  problems.  However,  the 
input  procedure  for  these  codes  requires  a  great  deal  of  time  on  the  part  of 
the  user.  On  the  other  hand,  the  input  procedure  for  the  approximate  three- 
dimensional  program  SANX  is  relatively  simple  since  it  is  similar  to  the  pro¬ 
cedure  for  two-dimensional,  axisymmetric  codes  such  as  SAA.S.4  After  the 
input  data  for  SANX  is  created  then  the  program  itself  generates  a  great  deal 
of  input  information  internally  which  contains  almost  all  the  information 
necessary  for  either  ADINA  or  NASTRAN  programs  except  this  information  is  not 
in  the  correct  form  for  these  codes . 

The  result  of  the  present  investigation  is  a  development  of  a  modified 
version  of  SANX  which  does  not  perform  stress  analysis,  but  rather  uses  its 
own  internal  finite  element  information  to  generate  the  input  data  for  ADINA 
or  NASTRAN  codes.  Actually,  two  different  versions  of  this  program  have  been 
developed,  one  for  ADINA  and  one  for  NASTRAN.  Developing  two  different  versions 
was  dictated  by  the  fact  that  the  input  procedures  for  ADINA  and  NASTRAN  codes 
are  of  completely  different  form  and,  consequently,  it  is  much  more  convenient 
to  generate  the  input  separately  for  each  program. 


1.  Zak,  A.  R. ,  Craddock,  J.  N.  and  Drysdale,  W.  H. ,  "Approximate  Finite- 
Element  Method  of  Stress  Analysis  of  Non-Axisymmetric  Configurations," 
International  Journal  of  Computers  and  Structures,  Vol.  9,  pp.  201-206, 
1978. 

2.  Bathe,  K.  J.,  "ADINA-A  Finite  Element  Program  for  Automatic  Dynamic 
Incremental  Nonlinear  Analysis,"  MIT  Acoustics  and  Vibration  Laboratory, 
Mechanical  Engineering  Department,  Report  82448-1,  1977,  Cambridge, 
Massachusetts . 

3.  McCormick,  C.  W.,  "MSC/NASTRAN  Basic  Training  Manual,"  MSR-46,  15  July 
1977,  (MacNeal-Schwendler  Corporation,  7442  N.  Figueroa  Street,  Los 
Angeles,  CA  90041). 

4.  Sawyer,  S.  G.,  "BRLESC  Finite  Element  Program  for  Axisymmetric,  Plane 
Strain,  and  Plane  Stress,  Orthotropic  Solids  with  Temperature -Dependent 
Material  Properties,"  BRL  Report  No.  1539  (March  1971).  AD  #727702 
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II.  MODIFIED  VERSIONS  OF  SANX 


Starting  with  the  original  SANX  program,  two  modified  versions,  one 
for  ADINA  and  one  for  NASTRAN,  were  developed.  The  two  programs  are  called 
ADIGEN  and  NASGEN,  and  the  flow  charts  for  these  codes  are  shown  in  Figures  1 
and  2,  respectively.  It  may  be  noted  that  the  two  flow  charts  are  similar 
except  ADIGEN  contains  the  Subroutine  ADINA  and  NASGEN  has  a  Subroutine 
NASTRN  in  the  same  place. 


The  first  step  in  generating  these  programs  was  to  truncate  the  SANX 
code  to  remove  the  portion  which  does  the  stress  analysis  and  retain  the 
parts  which  generate  nodal,  element,  and  loading  information.  Referring  to 
Figures  1  and  2,  it  can  be  seen  that  this  includes  all  the  subroutines  in 
these  flow  charts  except  ADINA  and  NASTRN.  The  latter  two  are  the  only  new 
subroutines  in  ADIGEN  and  NASGEN,  respectively.  After  truncating  the  SANX  code, 
the  only  changes  which  were  made  to  the  modified  program,  besides  adding  the 
extra  subroutine  to  each  version,  were  in  the  MAIN  program.  The  remaining 
subroutines  left  from  SANX  have  no  changes.  Descriptions  of  ADIGEN  and 
NASGEN  programs  are  given  in  the  following  sections. 


III.  PROGRAM  ADIGEN 


The  program  ADIGEN  is  designed  to  generate  three-dimensional  data  for 
ADINA  program.  The  listing  for  the  ADIGEN  program  is  given  in  Appendix  A. 

In  the  MAIN  portion  of  the  program  ADIGEN,  there  are  two  basic  changes  from 
SANX.  The  first  change  is  the  addition  of  six  temporary  files  on  which  the 
generated  data  for  ADINA  is  stored.  The  reason  for  six  tapes  is  that  SANX 
does  not  generate  all  the  nodal  data,  element  data,  etc.,  sequentially,  but 
rather  various  portions  of  this  data  are  generated  by  looping  through  each 
segment  of  the  structure.  However,  the  ADINA  input  data  has  to  have  the  data 
in  the  specific  order.  The  procedure  is,  therefore,  to  store  various  com¬ 
ponents  of  data  on  separate  files  and  then  combine  these  in  a  Procedure  File 
which  also  runs  the  program. 

The  six  temporary  files  and  their  functions  are  as  follows: 


TAPE2, 

ADNPD; 

nodal  point  data 

TAPE3, 

ADELD; 

element  data 

TAPE4, 

ADMAX; 

material  axes  orientation 

TAPE7, 

ADNPX ; 

extra  nodal  information  necessary  for 
material  axes  orientation 

TAPE8, 

ADMPC; 

material  property  cards 

TAPE9, 

ADECC; 

element  control  card 
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The  second  basic  change  in  the  MAIN  is  the  loading  information  of  the 
temporary  files  ADMPC  and  ADECC.  There  are  also  some  minor  changes  such  as 
addition  of  common  statement  TOTAL  and  dimensional  statement  which  contains 
two  new  quantities  which  are  later  set  in  the  program. 

The  Subroutine  ADINA  is  a  new  subroutine  and  essentially  creates  the 
remaining  files  ADNPD,  ADELD,  ADMAX,  and  ADNPX. 

After  the  ADIGEN  program  is  executed,  then  the  above  temporary  files 
are  assembled  into  one  file  called  the  ADI  which  represents  the  input  data 
for  ADINA  program.  The  file  ADI  is  obtained  by  assembling  the  temporary 
files  in  the  following  order: 

ADNPD 

ADNPX 

ADECC 

ADMPC 

ADMAX 

ADELD 

The  above  order  is  dictated  by  the  input  requirements  of  the  ADINA  program. 

In  order  to  illustrate  how  this  procedure  works  the  program  ADIGEN  was  applied 
to  a  test  example  and  the  results  will  be  presented  and  discussed  later  in 
this  report. 


IV.  PROGRAM  NASGEN 


The  program  NASGEN  is  similar  to  ADIGEN  in  the  sense  that  most  of  the 
changes  are  in  the  MAIN  program  and  one  new  subroutine,  NASTRN ,  has  been 
added.  The  listing  for  the  NASGEN  is  given  in  Appendix  B.  However,  one  big 
difference  between  ADIGEN  and  NASGEN  is  that  the  latter  program  contains  only 
one  additional  file,  TAPE1 ,  NAS,  as  opposed  to  six  for  the  ADIGEN.  The  file 
NAS  contains  the  generated  NASTRAN  data.  The  reason  why  the  N'ASTRAN  data  can 
be  created  with  one  file  only,  as  opposed  to  six  for  ADINA,  is  that  in  the 
input  procedure  for  NASTRAN,  the  data  cards  do  not  have  to  follow  any  specific 
order.  As  may  be  noted  from  the  Subroutine  NASTRN,  each  data  card  is  started 
with  a  label  which  identifies  the  type  of  data  on  that  card.  The  NASTRAN  does 
its  own  sorting  of  cards. 

Similarly  to  ADIGEN,  the  program  NASGEN  generates  some  of  the  data  in 
the  MAIN  program  and  the  rest  in  Subroutine  NASTRN.  In  order  to  illustrate 
how  the  NASGEN  program  is  used,  a  numerical  test  example  will  be  given  and 
discussed  later  in  this  report. 
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V.  INPUT  PROCEDURE  FOR  ADIGEN  AND  NASGEN  PROGRAMS 


The  programs  ADIGEN  and  NASGEN  accept,  directly,  the  input  data  for 
the  elastic  version  of  the  SANX  program.  The  procedure  for  developing  the 
input  data  is  described  in  Reference  1.  In  order  to  illustrate  the  procedure 
of  using  these  programs,  two  numerical  examples  were  executed  and  the  cor¬ 
responding  input  data  for  ADINA  and  NASTRAN  was  generated.  The  numerical 
examples  are  discussed  below. 


VI.  NUMERICAL  EXAMPLE  -  ADIGEN 


The  ADIGEN  program  was  applied  to  a  numerical  example  which  corresponds 
to  approximately  60°  section  of  a  non-axisymmetric  cylinder.  The  actual  con¬ 
figuration  has  no  particular  significance  except  that  the  data  was  available 
from  previous  investigations  on  the  SANX  program.  The  schematic  of  the  physical 
configuration  used  in  the  numerical  example  is  shown  in  Figure  3.  The  confi¬ 
guration  is  referred  to  a  cylindrical  coordinate  system  as  required  by  the 
original  SANX  program. 

For  the  purpose  of  the  analysis  the  configuration  is  divided  into  four 
segments  in  the  r-6  plane.  This  is  illustrated  in  Figure  4.  The  first  three 
segments  are  identical  and  the  fourth  is  different.  The  first  three  segments 
are  represented  by  segment  Type  No.  1  and  the  fourth  by  Type  No.  2.  The 
finite  element  models,  in  the  r-z  plane,  for  the  two  types  of  segments  are 
illustrated  in  Figure  5  in  the  I-J  grid  configuration.  The  I-J  grid  is  used 
in  the  SANX  program  and,  therefore,  in  ADIGEN  as  the  basis  for  numbering  of 
the  nodal  points  and  the  finite  elements. 

The  above  configuration  was  used  in  the  ADIGEN  program  to  generate  the 
ADINA  input  file  ADI.  The  listing  for  the  ADI  file  is  given  in  Appendix  C. 

It  is  useful  to  identify  various  portions  of  this  file  which  correspond  to 
the  various  parts  of  the  ADINA  input  data.  In  order  to  facilitate  this  dis¬ 
cussion  the  start  of  each  part  of  the  data  is  identified  in  the  right  hand 
side  by  the  original  name  of  the  temporary  file  as  listed  in  the  first  part 
of  this  report. 

Referring  to  Appendix  C,  the  file  ADI  is  divided  into  the  following 
parts: 

ADNPD  Nodal  Point  Data 

This  data  gives  nodal  point  numbers  and  coordinates  for  the  nodes 
which  are  later  used  to  define  finite  elements.  The  coordinates  are 
in  the  cylindrical  coordinate  system  and,  therefore,  the  letter  X 
appears  in  front  of  each  line.  The  scheme  for  generating  numbers  is 
as  follows.  Referring  to  Figure  4  it  can  be  seen  that  the  four  segments 
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are  separated  by  five  planes.  Each  of  these  planes  is  in  the  r-z 
coordinate  plane.  The  nodes  in  each  of  these  planes  are  numbered 
starting  from  multiples  of  200.  For  example,  for  the  first  plane  the 
nodes  start  from  201,  the  second  from  401,  and  the  fifth  from  1001. 

In  any  plane  the  nodes  are  numbered  according  to  the  I-J  grid  coordinate 
illustrated  in  Figure  3.  The  nodes  are  numbered  in  multiples  of  25 
depending  on  the  value  of  J.  For  example,  in  the  first  plane,  for  J=1 
the  numbering  starts  from  201,  for  J=2  from  226,  and  so  on.  Using  this 
scheme  it  is  easy  to  identify  the  actual  position  of  the  node  by  its 
number  relative  to  the  r-z  plane  and  the  I-J  grid. 


ADNPX  Extra  Nodes  for  Defining  Material  Axes  Orientation 

For  each  finite  element,  which  appears  later  in  the  ADI  file, 
there  are  three  extra  nodes  which  define  the  material  axes  orienta¬ 
tions.  This  is  done  by  using  number  7000,  8000,  and  9000  and  adding 
to  them  the  element  number.  For  example,  for  element  201  the  three 
nodes  are  7201,  8201,  and  9201. 


ADECC  Element  Control  Card 

This  card  is  self  explanatory  according  to  the  ADINA  User's  Manual. 


ADMPC  Material  Property  Cards 

There  are  three  materials  used  in  the  example  and  for  each  of  these 
there  are  two  material  property  cards  as  required  by  ADINA. 


ADMAX  Material  Axes  Orientation 

For  each  element  one  card  is  generated  which  contains  an  element 
ID  number,  later  repeated  in  the  element  data  card,  and  the  three  nodes 
defining  the  material  orientation  according  to  the  scheme  required  by 
ADINA. 


ADELD  Element  Data 

For  each  element  there  are  two  cards  generated.  The  first  card 
defines  the  material  information,  according  to  ID  number  from  ADMAX 
data,  and  the  second  card  contains  eight  nodal  points  for  the  three- 
dimensional  finite  element.  It  may  be  noted  that  the  elements  are 
numbered  sequentially.  This  is  done  by  starting  with  segment  No.  1 
and  using  the  scheme  from  SANX  in  which  the  elements  are  numbered  first 
in  the  I  direction  followed  by  variation  in  the  J  direction. 
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VII.  NUMERICAL  EXAMPLE  -  NASGEN 


The  NASGEN  program  has  been  applied  to  a  numerical  example  which  cor¬ 
responds  to  a  90°  section  of  a  non-axisymroetric  cylinder.  The  geometry  of 
the  cylinder  used  in  this  numerical  example  is  shown  in  Figure  6.  The  con¬ 
figuration  is  referred  to  a  (r,  z,  6)  cylindrical  coordinate  system. 

For  analysis  purposes  the  configuration  was  divided  into  five  segments 
in  the  r-0  plane.  The  division  into  the  segments  is  shown  in  Figure  7.  The 
first  and  fifth  segments  are  identical,  the  second,  third,  and  fourth  are 
different  but  identical.  The  first  and  fifth  segments  are  represented  by 
segment  Type  No.  1,  whereas  the  second,  third,  and  fourth  segments  are  represented 
as  Type  No.  2.  The  finite  element  models  in  the  r-z  plane  for  the  two  different 
types  of  segments  are  illustrated  in  Figure  8  in  the  I-J  grid  configuration. 

The  loading  on  the  cylinder  is  assumed  to  be  composed  of  an  axial  and  radial 
pressure  of  30  x  10^  psi  which  is  illustrated  as  acting  on  the  different  segments 
shown  in  Figure  8.  The  axial  pressure  is  assumed  to  cause  axial  acceleration 
which  has  been  calculated  to  match  the  external  pressure  load  assuming  aluminum 
material . 

The  NAS  data  generated  for  this  example  is  shown  in  Appendix  D.  It  is 
useful  to  identify  the  various  types  of  data  in  this  file  by  referring  to 
Appendix  D.  The  NAS  data  file  contains  the  following  generated  input  cards 
for  the  NASTRAN  program: 


GRID  cards; 

These  cards  define  the  cylindrical  coordinates  for  the  grid  points. 

The  scheme  used  to  number  the  grids  is  as  follows.  Referring  to  Figure  7  it 
can  be  observed  that  the  five  segments  are  separated  by  six  planes  in  the 
r-z  plane.  The  grid  points  in  each  of  these  planes  are  numbered  starting 
from  multiples  of  2000.  For  example,  starting  with  the  plane  to  the  left  of 
Segment  No.  1,  the  first  grid  point  will  be  2001.  Taking  these  planes  in 
geometrical  order,  the  first  grid  point  in  the  plane  between  Segment  No.  1 
and  No.  2  will  be  4001,  and  so  on.  In  addition  to  this  scheme,  in  each 
plane  the  grid  points  are  also  numbered  in  the  multiples  of  25  depending  on 
the  value  of  J.  For  example,  in  the  first  plane,  for  J=1  the  numbering 
starts  with  2001,  for  J=2  the  first  grid  is  2026,  and  so  on.  Using  this 
scheme  it  is  easy  to  identify  the  position  of  each  grid  point. 

CHEXA  cards; 

These  cards  define  the  three-dimensional,  isoparametric,  finite  elements. 
Each  CHEXA  has  a  continuation  card  identified  by  the  NASTRAN  procedure.  Each 
pair  of  cards  identifies  the  element  number,  material,  and  the  eight  connecting 
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grid  points.  The  scheme  of  numbering  the  elements  is  according  to  the 
segment  number.  Starting  with  Segment  No.  1  of  Figure  7,  the  first  element 
is  numbered  2001,  in  the  Segment  No.  2  with  4001,  and  so  on  with  multiples 
of  2000.  In  the  segment  the  elements  are  numbered  sequentially. 

PLOAD  cards; 

These  cards  define  the  pressure  load  cards  relative  to  four  grid  points. 
MAT1  card; 

This  card  defines  the  material  properties. 

CORD 1C  card; 

This  card,  together  with  three  additional  GRID  cards  numbered  101,  102, 
and  103,  defines  the  cylindrical  coordinate  orientation  relative  to  the  basic, 
Cartesian  coordinate  system. 

GRAV  card; 

This  defines  the  axial  acceleration  load. 

PSOLID  card; 

This  card  controls  the  output  parameters  for  the  CHEXA  elements. 


VIII.  CONCLUDING  REMARKS 

This  report  describes  two  programs,  ADIGEN  and  NASGEN,  which  are  designed 
for  automatic  generation  of  input  data  for  the  general  purpose  finite  element 
structural  codes  ADINA  and  NASTRAN  respectively.  The  two  programs,  ADIGEN  and 
NASGEN,  have  been  developed  by  modifying  SANX  structural  program  which  was 
previously  developed  at  the  University  of  Illinois.  Using  these  two  programs, 
numerical  examples  were  executed.  The  result  of  these  examples  is  input  data 
for  ADINA  and  NASTRAN  programs.  The  data  for  the  NASTRAN  program  was  checked- 
out  by  executing  the  NASTRAN  program.  Similar  check  was  not  performed  for  the 
ADINA  program  since  this  program  is  not  available  at  the  University  of  Illinois. 
In  conclusion,  it  may  be  observed  that  the  two  new  programs,  ADIGEN  and  NASGEN, 
generate  a  great  deal  of  input  data  for  the  ADINA  and  NASTRAN  codes.  In  each 
case  this  data  was  generated  from  approximately  20  to  30  input  cards  required 
by  the  SANX  code.  Consequently,  the  automatic  input  data  generation  has  a  great 
potential  of  reducing  the  necessary  effort  in  using  the  ADINA  and  NASTRAN 
programs . 
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Figure  1.  Flow  Chart  for  ADIGEN  Computer  Program 
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Figure  2.  Flow  Chart  for  NASGEN  Computer  Program 
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Figure  3.  Cylindrical  Configuration  Used  in  the  Numerical  Fxai™le 


17 


\  I 

\  // 

\  / 


Figure  4.  Division  of  Cylindrical  Configuration  into  Four  Segments 
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Figure  5.  The  I-J  Grid  Used  in  the  Finite  Element  Model 
for  Segment  Types  No.  1  and  No.  2 
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Figure  6.  Nonaxisymmetric  Cylindrical  Configuration 
Used  in  the  Numerical  Example 
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Figure  7.  Division  of  Cylindrical  Configuration  Into  Five  Segments 
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Figure  8.  The  I-J  Grid  for  the  Segment  Types  1  and  2 
and  the  External  Pressure  Load 
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APPENDIX  A 


Listing  of  the  ADIGEN  Program 
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nnnno 


PROGRAM  ADIGEN ( INPUT , OUTPUT , ADNPD , ADNPX , ADECC , ADMPC , ADMAX , ADELD 

1 , ADPRL, TAPE5=INPUT, TAPE6=OUTPUT, TAPE2=ADNPD, TAPE7=ADNPX , TAPE9= 

2  ADECC , TAPE8  = ADMPC , TAPE4 = ADMAX , TAPE3 = ADE LD , TAPEl = ADPRL ) 

INTEGER  CODE 


COMMON /NPDATA/R (1000) ,CODE(1000) ,XR(1000) ,Z(1000) ,XZ(1000) , 

1NPNUM (25,80) ,T(1000) ,XT(1000) 

COMMON/ARG/RRR ( 5 ) ,ZZZ(5) ,RR(4) ,ZZ(4) ,S(15,15) ,P(15) ,TT(6) , 

1H (6,15) ,CRZ (6,6) ,XI(10) , ANGLE (4) ,SIG(18) ,EPS(18) ,N 
COMMON /ELDATA/BETA (1000) ,EPR(1000) ,PR(200) ,SH(200) ,IX(1000,5) , 

1 IP ( 200) , JP ( 200) , IS ( 200) , JS (200) ,ALPHA(1000) ,IT(200) , JT(200) , 

2ST (200) , INP (1000) ,JNP(1000) 

COMMON /BASIC/ACELZ , ANGVEL , ANGACC , TREF , VOL , NUMNP , NUMEL, NUMPC , NUMSC , 
lNUMST 

COMMON/NXMESH/THETAN ( 4 ) ,NST(4) ,NUMS(4,5) ,NPC(8,8) 
COMMON/ANSl/NUMELS (4) ,NUMNPS(4) 

COMMON/NXDATA/NTP , NTYPS, NTS , NTOTS 
COMMON/NONAXI/S1 (30,30) , PI (30) , THETA, BS1 (6,30) 

COMMON/ SOLVE/X (4428) ,Y(4428) ,TEM(4428) , NUMTC , MB AND 

COMMON/TD/IMIN (100) ,IMAX(100) ,JMIN(25) ,JMAX(25) , MAXI ,MAXJ, NMTL, NBC 

C OMMON/CONVRG/ I DONE 

COMMON/PLANE/NPP 

COMMON/RESULT/BS ( 6 , 15)  ,D(6,6)  , C  ( 6 , 6 )  ,AR,BB(6,9)  ,CNS(6,6) 
COMMON/MATP/RO ( 6 ) ,E(12,16,6) ,EE(16) ,AOFTS(6) 

COMMON/TOTAL/NENL 
DIMENSION  TITLE (20) , FTHETA ( 8 ) 


WRITE  ADINA  DATA  BLOCK  HEADINGS 


*********************************************** 


WRITE (1,5001) 
WRITE (2,5002) 
WRITE (3 , 5003) 
WRITE(4 ,5004) 
WRITE (7 ,5007) 
WRITE(8, 5008) 
WRITE(9, 5009) 

5001  FORMAT( "****** 

5002  FORMAT ( n ****** 

5003  FORMAT("****** 

5004  FORMAT("****** 

5007  FORMAT!"****** 

5008  FORMATC****** 

5009  FORMAT! "****** 
c*  ********* 


3/d  pressure  load  data 

NODAL  POINT  DATA 
3/D  ELEMENT  DATA 
MATERIAL  AXES  ORIENTATION 
EXTRA  NODES  FOR  MATERIAL  AXES 
MATERIAL  PROPERTY  CARDS 
ELEMENT  COTROL  CARD 
**************** 


C  READ  AND  WRITE  CONTROL  INFORMATION 

C*  ************************* 


READ (5, 3000)  NTYPS, NTOTS 
DO  150  1=1, NTYPS 

150  READ (5,3001)  THETAN ( I) ,NST ( I) 
DO  151  1=1, NTYPS 
J2  =NST(I) 


*****«) 
******  it ) 

******  ii ) 

******«) 

******!•) 

******••) 

******!•) 

*  *  *  *  * 

*  *  *  *  * 


* 


* 


* 


* 


*  *  * 

*  *  * 
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151  READ (5,3000)  (NUMS ( I , J) , J=1 , J2) 

DO  152  1=  1, NTOTS 

152  READ (5,3002)  (NPC ( I , J) , J=1 , 8) 

3000  FORMAT ( 8110) 

3001  FORMAT (F10 .5,110) 

3002  FORMAT (8110) 

NFACE=NTOTS+l 
FTHETA ( 1) =0.0 

DO  4100  1=1 , NTOTS 
DO  4000  J=1 ,NTYPS 
K2=NST ( J) 

DO  4000  K=1,K2 

4000  IF(NUMS ( J,K) .EQ. I) GO  TO  4100 
4100  FTHETA ( 1+1) =FTHETA ( I) +THETAN ( J) 

WRITE (6,4200) ( (FTHETA ( I) ,1) ,1=1, NTOTS) 

4200  FORMAT (  "  1 "  ,  "  ANGLE  IN  DEGREES  ",5X,"  BEGINNING  SEGMENT  #"/ 

1  (E20 . 10 , 110) ) 

WRITE(6, 3010) 

3010  FORMAT ( "1" , "SEGMENT  DATA  FOR  NONAX I SYMMETRIC  PROBLEM") 

WRITE(6, 3011)  NTYPS, NTOTS 

3011  FORMAT  ("  NUMBER  OF  TYPES  OF  SEGMENTS  =  ",I5,//, 

1  "  NUMBER  OF  TOTAL  SEGMENTS  =",I5) 

DO  153  1=1, NTYPS 
WRITE (6,3012)  I , THETAN ( I) ,NST ( I) 

3012  FORMAT ( "  ",///,"  SEGMENT  TYPE  =  ",I5/,"  THETA  =  ",Fl0.5/, 

1  "  NUMBER  OF  SEGMENTS  OF  THIS  TYPE  =  ",I5) 

J2  =  NST(I) 

WRITE(6, 3013) (NUMS(I,J) ,J=1,J2) 

3013  FORMAT ( "  " , "  SEGMENT  NUMBERS  IN  GLOBAL  SYSTEM  ARE  ",5I5) 

153  CONTINUE 

DO  154  1=1, NTOTS 

154  WRITE (6, 3014) I, (NPC (I, J) , J=l,8) 

3014  FORMAT ("  " , "CONNECTING  NODES  FOR  SEGMENT" , 15 , "  ARE", 8 15) 

NENL=0 

DO  950  NTP  =  1, NTYPS 
THETA=  THETAN (NTP)  /57. 295780 

50  READ ( 5 , 1000  ) TITLE, NNLA, NUMTC , NUMMAT, NUMPC, NUMSC , NUMST, TREF 
1 , INERT , NLINC , INCI , INCF , IPLOT 

WRITE ( 6 , 2  0  0  0 ) TITLE , NNLA , NUMTC , NUMMAT , NUMPC , NUMSC , NUMST , TREF , INERT , 
lNLINC 
NPP=0 

C*  *********************************** 
C  GENERATE  FINITE  ELEMENT  MESH 

C*  *********************************** 

100  CALL  MESH 

NUMELS (NTP)  =  NUMEL 
NUMNPS (NTP)  =  NUMNP 

C*  *********************************** 

C  READ  AND  WRITE  T5MPERATURE  DATA 

C*  *********************************** 

103  IF ( NUMTC . EQ . 0 )  GO  TO  440 

I F ( NUMTC . GT . 0 )  READ(5,1001)  (X ( I) , Y ( I) , TEM ( I) , 1=1 , NUMTC) 
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o  o  o 


I F ( NUMTC . EQ . - 2 )  CALL  TEM2 ( NUMNP ) 

IF (NUMTC . EQ . -2 )  GO  TO  440 
MPRINT=0 

DO  210  1=1, NUMTC 
IF (MPRINT. NE. 0)  GO  TO  200 
WRITE(6, 2001) 

MPRINT=59 

200  MPRINT=MPRINT-1 
210  WRITE (6,2002)  X ( I ) , Y ( I ) , TEM ( I) 

MPRINT=0 

DO  230  N=l, NUMNP 
IF (MPRINT. NE. 0)  GO  TO  220 
WRITE (6,2003) 

MPRINT=59 

220  MPRINT=MPRINT-1 

CALL  TEMP ( R (N) ,Z(N) ,T(N) ) 

230  WRITE (6,2004)  N,R(N) ,Z (N) ,T(N) 

440  MPRINT=0 

DO  460  N=1 ,NUMEL 
IF (MPRINT. NE . 0)  GO  TO  450 
WRITE(6, 2008) 

MPRINT=5  9 

450  MPRINT=MPRINT-1 
II=IX (N, 1) 

JJ=IX (N, 2) 

KK=IX (N, 3) 

LL=IX (N, 4) 

TEM  IS  TEMPORARY  STORAGE  FOR  ELEMENT  TEMPERATURES 

TEM(N)  =  (T  ( II)  +T(  JJ)  +T(KK)  +T(LL)  )/4.00 
460  WRITE(6, 2009)  N, ( IX (N, I) , 1=1 , 5) ,BETA( N) , ALPHA (N) ,TEM (N) 

DO  470  K=1,NUMEL 
470  T (K) =TEM (K) 

C*  *********************************** 
C  READ  AND  WRITE  MATERIAL  PROPERTIES 

C*  *********************************** 
500  CONTINUE 

DO  510  M=1 ,NUMMAT 

READ (5, 1004)  MTYPE,  (NT, RO (MTYPE)  , AOFTS (MTYPE)  ) 

WRITE(6, 2010)  MTYPE, NT, RO (MTYPE) 

READ  (5,1005)  (  (E(I,J, MTYPE)  ,J=1,14)  ,1=1, NT) 

IF (AOFTS (MTYPE) .NE.l.)  WRITE ( 6 , 20 11) ( (E ( I , J , MTYPE) , J=l,13) ,1=1, NT) 
IF (AOFTS (MTYPE) .EQ.l.)  WRITE (6 , 2012) ( (E (I, J, MTYPE) ,J=1,13) ,1=1, NT) 
DO  510  I=NT, 12 
DO  510  J=1 , 16 

510  E (I, J, MTYPE) =E(NT,J, MTYPE) 

DO  900  NL=1 ,NLINC 
WRITE (6,2030)  NL 
ACELZ=0 .00 
ANGVEL=0 . 00 
ANGACC=0 .00 
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noon 


IF( INERT  . EQ.  0)  GO  TO  511 

IF (NL  .NE.  1  .AND.  INCI  .EQ.  0)  GO  TO  511 

c*********************************************************************** 

C  READ  AND  WRITE  DYNAMIC  FORCES 

C*****************************  *******************************  *********** 

READ (5,1030)  ACELZ ,  ANGVEL,  ANGACC 
WRITE (6, 2031)  ACELZ,  ANGVEL,  ANGACC 
511  CONTINUE 

£>*  *********************************** 
C  READ  AND  WRITE  PRESSURE  AND  SHEAR  BOUNDARY  CONDITIONS 

£*  *********************************** 
IF ( NL  .NE.  1  .AND.  INCF  .EQ.  0)  GO  TO  700 
600  IF (NUMPC. EQ. 0)  GO  TO  630 
MPRINT=0 

DO  620  L=l, NUMPC 
IF (MPRINT.NE. 0)  GO  TO  610 
WRITE(6, 2013) 

MPRINT=58 

610  MPRINT=MPRINT-1 

READ (5,1006)  IP ( L) , JP ( L) , PR ( L) 

620  WRITE(6, 2014)  IP ( L) , JP ( L) , PR ( L) 

630  I F ( NUMSC . EQ . 0 )  GO  TO  701 
MPRINT=0 

DO  650  L=l, NUMSC 
IF (MPRINT. NE. 0)  GO  TO  640 
WRITE(6, 2015) 

MPRINT=58 

640  MPRINT=MPRINT-1 

READ (5,1006)  IS (L) , JS (L) ,SH (L) 

650  WRITE(6, 2014)  IS (L) , JS (L) ,SH (L) 

701  IF (NUMST. EQ. 0)  GO  TO  700 
MPRINT=0 

DO  680  L=l, NUMST 
IF (MPRINT. NE . 0)  GO  TO  670 
WRITE(6, 2025) 

MPRINT=58 

670  MPRINT=MPRINT-1 

READ  (5,1006)  IT  ( L)  ,  JT  (L)  ,  ST  (L) 

680  WRITE (6, 2014) IT(L) ,JT(L) ,ST(L) 

************************************ 

DETERMINE  BANDWIDTH,  INITIALIZE  ELASTIC-PLASTIC  RATIO, 

AND  CONVERT  BETA  FROM  DEGREES  TO  RADIANS 
************************************ 

700  J=0 

DO  710  N=1 , NUMEL 
IX (N , 5) =IABS ( IX (N, 5) ) 

DO  710  1=1,4 
DO  710  L=1 , 4 
KK=IABS  (IX(N,I)-IX(N,L)  ) 

IF(KK.GE.J)  J=KK 
710  CONTINUE 

MBAND=3  * J+3 
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non  n  on 


IF (NL. GT. 1)  GO  TO  721 
DO  720  N=1,NUMEL 
EPR(N) =1. 

ALPHA (N) =ALPHA (N) /57. 295780 

720  BETA(N) =BETA(N) /5 7.295780 

721  CONTINUE 
900  CONTINUE 

CALL  ADINA (NFACE , FTHETA) 

950  CONTINUE 

WRITE  ELEMENT  CONTROL  ON  FILE  9 
NZ=0 

WRITE (9, 4400)  NENL, NZ , NZ , NZ , NUMMAT, NZ , NENL, NZ 
4400  FORMAT (3X, "3", 314 , 11X, "8 " , 20X, 14 , 7X, "2" , 314 , 4X, 14) 

WRITE  MATERIAL  PROPERTIES  ON  FILE  8 

DO  960  1=1, NUMMAT 
WRITE  (8,4500)  I,RO(I) 

4500  FORMAT ( 15, F10. 5) 

WRITE  (8, 4  600)  (E  ( 1 ,  J,  I)  ,  J=2 , 5)  ,  E  ( 1 , 7 , 1)  ,  E  ( 1 , 6  ,  I)  , 
IE (1,8,1) , E ( 1 , 10 , I ) 

WRITE(8, 4600)  E(1,9,I) 

960  CONTINUE 

4600  FORMAT ( 3F10 . 0 , 3F10 . 3, 2F1 0.0) 

910  CONTINUE 

1000  FORMAT ( 20A4 /6 15, F5. 0,515) 

1001  FORMAT (3F1 0.0) 

1004  FORMAT  (2l5,2Fl0.0) 

1005  FORMAT (7 FI 0 . 0) 

1006  FORMAT  (2I5,F10.0) 

1030  FORMAT (3F1 0.0) 

2000  FORMAT  (2H1  ,20A4/ 


1  33H0  NUMBER  OF  APPROXIMATIONS - 14/ 

2  33H0  NUMBER  OF  TEMPERATURE  CARDS - 14/ 

3  33H0  NUMBER  OF  MATERIALS - 14/ 

4  33H0  NUMBER  OF  PRESSURE  CARDS - 14/ 

5  33H0  NUMBER  OF  SHEAR  CARDS - 14/ 

6  33H0  NUMBER  OF  TORSION  CARDS - 14/ 

7  33H0  REFERENCE  TEMPERATURE - E12.4/ 

8  33H0  NUMBER  OF  INERTIA  CARDS - 14/ 

9  33H0  NUMBER  OF  LOAD  INCREMENTS - 14/) 

2001  FORMAT  ( lHl , 13X , 1HR, 14X , lHZ , 14X , 1HT) 

2002  FORMAT  (3F15.3) 

2003  FORMAT  (35H1  N  R  Z  T) 

2004  FORMAT  ( 15 , 2F10 . 4 , F10 . 0) 


2008  FORMAT  (74H1  EL  I  J  K  L  MATERIAL  ANGLE  BETA 
1LPHA  TEMPERATURE) 

2009  FORMAT  ( 15 , 414 , 18 , Fll . 1, 2F13 . 3) 

2010  FORMAT  ( 1H1 , "MATERIAL  IDENTIFICATION  NUMBER  =",12/ 
llH  , "NO.  OF  MATERIAL  TEMPERATURE  CARDS  =",12/ 


t 

ANGLE  A 
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o  n 


2 1H  , "MASS  DENSITY  =",E15.7) 

2011  FORMAT  ( lH  /"TEMPERATURE  =",E15.7/ 

11H  /"MODULUS  OF  ELASTICITY-EN  =",E15.7/ 

21H  /"MODULUS  OF  ELASTICITY-ES  =",E15.7/ 

3 1H  /"MODULUS  OF  ELASTICITY-ET  =" ,E15.7/ 

41H  /"POISSON  RATIO-NUNS  =",E15.7/ 

51H  ,"POISSON  RATIO-NUNT  =",E15.7/ 

61H  /"POISSON  RATIO-NUST  =",E15.7/ 

7 lH  , "SHEAR  MODULUS-GNS  =",E15.7/ 

8 1H  ,"SHEAR  MODULUS-GST  =" ,E15.7/ 

91H  ,"SHEAR  MODULUS-GTN  =",E15.7/ 

11H  /"COEFFICIENT  OF  THERMAL  EXPANSION-AN  =",E15.7/ 

2 lH  /"COEFFICIENT  OF  THERMAL  EXPANSION-AS  =",E15.7/ 

3 lH  /"COEFFICIENT  OF  THERMAL  EXPANSION-AT  =",E15.7/) 

2012  FORMAT  ( lH  /"TEMPERATURE  =",E15.7/ 
llH  /"MODULUS  OF  ELASTICITY-EN  =",E15.7/ 

2 lH  /"MODULUS  OF  ELASTICITY-ES  =",E15.7/ 

3 lH  /"MODULUS  OF  ELASTICITY-ET  =",E15.7/ 

41H  /"POISSON  RATIO-NUNS  =",E15.7/ 

51H  /"POISSON  RATIO-NUNT  =",E15.7/ 

6 lH  ,"POISSON  RATIO-NUST  =",E15.7/ 

7 lH  ,"SHEAR  MODULUS-GNS  =",E15.7/ 

8 lH  ,"SHEAR  MODULUS-GST  =" ,E15.7/ 

91H  ,"SHEAR  MODULUS-GTN  =",E15.7/ 

11H  ,"FREE  THERMAL  STRAIN-FN  =",E15.7/ 

21H  ,"FREE  THERMAL  STRAIN-FS  =",E15.7/ 

3 lH  ,"FREE  THERMAL  STRAIN-FT  =",E15.7/) 

2013  FORMAT  (30H1  PRESSURE  BOUNDARY  CONDITIONS/2  OH  I 

2014  FORMAT  (2I5,F10.1) 

2015  FORMAT  (27H1  SHEAR  BOUNDARY  CONDITIONS/17H  I  ^ 

2016  FORMAT  (26H  THE  SYSTEM  CONVERGED  IN  I2,llH  ITERATIONS) 

2017  FORMAT  (33H  THE  SYSTEM  DID  NOT  CONVERGE  IN  12 , 11H  ITERATIONS) 

2024  FORMAT  (43H0  THE  AXISYMMETRIC  OPTION  NAS  BEEN  SELECTED) 

2025  FORMAT ( 30H1  TORSION  BOUNDARY  CONDITIONS/17H  I  J  SHEAR) 

2030  FORMAT ( 1H1 , "LOAD  STEP=",I4) 

2031  FORMAT ( 1H0  , "AXIAL  ACCELERATION  =",E12.4/ 

11H0  /"ANGULAR  VELOCITY  =",E12.4/ 

21H0  /"ANGULAR  ACCELERATION3 " , E12 . 4 ) 

920  STOP 


PRESSURE) 
SHEAR) 


END 

SUBROUTINE  ANGLE  ( R, Z , RC , ZC , ANG) 

FIND  ANGLE  OF  INCLINATION  BETWEEN  0  AND  2*PI 

****************************4,^^^^^^^ 

PI=3. 1415927 
Dl= ( Z-ZC) 

D2= (R-RC) 

IF(ABS(R-RC) .GT.l.E-8)  GO  TO  100 
ANG=PI/2 . 

IF(Dl .GT. 1. E-8)  RETURN 

ANG=-ANG 

RETURN 

c ************************************ 
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C  ALLOW  CIRCLE  TO  CROSS  AXIS 

C*  *********************************** 
100  ANG= ATAN2 ( Dl , D2 ) 

RETURN 

END 

SUBROUTINE  CIRCLE ( ANGl , DELPHI , RSTRT, ZSTRT, RC, ZC, I , J) 

INTEGER  CODE 

COMMON /TD/IMIN (100) ,IMAX(100) ,JMIN(25) ,JMAX(25) , MAX I , MAX J , NMTL , NBC 
COMMON/NPDATA/R( 1000)  ,CODE(1000)  ,XR(1000)  ,Z(100  0)  ,XZ(1000)  , 

1NPNUM (25,80) ,T(1000) ,XT(1000) 

DIMENSION  AR(25, 80) ,AZ (25,80) 

EQUIVALENCE  ( R ( 1) , AR) , ( Z ( 1) , AZ) 

C*  *********************************** 

C  FIND  INTERSECTION  OF  LINE  AND  CIRCLE  =  NEW  R  AND  Z 

0*  *********************************** 
ANG1=ANG1+DELPHI 

RR=SQRT ( (RSTRT-RC) * *2+ ( ZSTRT-ZC) **2) 

AR ( I , J ) =RC+RR*COS (ANGl) 

AZ (I, J) =ZC+RR*SIN ( ANGl) 

RETURN 

END 

SUBROUTINE  MESH 
INTEGER  CODE 

DIMENSION  AR (25,80) ,AZ(25,80) , NCODE ( 25 , 80) 

COMMON/TD/IMIN (100) ,IMAX(100) ,JMIN(25) , UMAX (25) , MAX I, MAX J, NMTL, NBC 
COMMON/NPDATA/R( 1000) ,CODE(1000) ,XR(1000) ,Z(1000) ,XZ(1000) , 

1NPNUM (25,80) ,T(1000) ,XT(1000) 

COMMON/E LDATA/BETA (1000) ,EPR(1000) ,PR(200) ,SH(200) ,IX(1000,5) , 

IIP ( 200)  , JP ( 2  0  0 )  , IS ( 200)  , JS ( 200)  ,ALPHA(1000)  ,IT(200)  ,JT(200)  , 
2ST(200) ,INP(1000) ,JNP(1000) 

EQUIVALENCE  (R ( 1 )  , AR)  ,  (Z ( 1)  , AZ)  ,  ( IX ( 1 , 1)  , NCODE) 

0*  *********************************** 
C  MESH  CONTROL  INFORMATION 

0*  *********************************** 
READ  (5,1000)  MAXI, MAXJ,NSEG, NBC, NMTL 
WRITE(6, 2000)  MAXI , MAX J, NS EG , NBC , NMTL 
0*  *********************************** 
C  INITIALIZE 

0*  *********************************** 
ISEG=-1 
PI=3. 1415927 
DO  110  J=l,100 
DO  100  1=1,25 
NCODE (I, J) =0 
AR(I, J) =0. 

AZ (I, J) =0. 

JMAX (I)  =0 
100  JMIN ( I) =MAXI 
IMIN ( J) =MAXJ 
110  IMAX ( J) =0 

0*  *********************************** 
C  LINE  SEGMENT  CARDS 
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C*  **************************** 

150  ISEG=ISEG+1 

159  IF ( ISEG. EQ. NSEG)  GO  TO  400 

READ (5, 1001)  II, Jl,Rl,Zl,I2, J2 , R2 , Z2 , 13 , J3 , R3 , Z3 , IPTION 

WRITE (6, 2001) Il,Jl,Rl,Zl,l2,J2,R2,Z2,l3,J3,R3,Z3, IPTION 

IPTION=IPTION+l 

AR ( 11 , Jl) =Rl 

AZ (II, Jl) =Zl 

NCODE(Il, Jl) =1 

CALL  MNIMX ( 11 , Jl) 

GO  TO  (150,200,200,300,300,200,200),  IPTION 
C*  *********************************** 
C  GENERATE  STRAIGHT  LINES  ON  BOUNDARY 

0*  *********************************** 
200  DI=  ABS (FLOAT ( 12-11) ) 

DJ=  ABS ( FLOAT ( J2-J1 ) ) 

AR(I2,J2)=R2 
AZ (I2,J2) =Z2 
NCODE ( 12 , J2 ) =1 
CALL  MNIMX (I 2, J2) 

ISTRT=Il 
ISTP=I2 
JSTRT= Jl 
JSTP=J2 

DIFF=MAXl (DI , DJ) 

ITER=DIFF-1 . 

I INC=0 
JINC=0 

IF(I2.NE. II)  IINC=( 12-11) /IABS( 12-11) 

IF(J2.NE. Jl)  JINC=(J2-J1)/IABS(J2-J1) 

KAPPA=1 

IF (12 .NE. II . AND. J2 .NE. Jl . AND. IPTION. NE. 3)  KAPPA=2 
IF (KAPPA. EQ . 2)  DIFF=2 . *DIFF 
RINC= (R2-R1) /DIFF 
ZINC=(Z2-Z1)/DIFF 

WRITE(6, 2002)  DI , DJ, DIFF, RINC , Z INC , ITER, I INC , JINC , KAPPA 
CHECK  FOR  INPUT  ERROR 

IF (KAPPA. NE. 2 . OR. DI . EQ . DJ)  GO  TO  210 
WRITE(6, 2003) 

GO  TO  150 

INTERPOLATE 


210  1=11 
J=J1 

WRITE(6, 2004) 

DO  230  M=1 , ITER 

IF(ITER. EQ.0. AND. IPTION. EQ. 2)  GO  TO  230 
IF(ITER.EQ. 0. AND. IPTION. EQ. 6)  GO  TO  230 
IF(ITER. EQ.0. AND. IPTION. EQ. 7)  GO  TO  230 
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IF (KAPPA. EQ. 2)  GO  TO  220 

IOLD=I 

1=1+1 INC 

JOLD= J 

J= J+JINC 

AR( I , J) =AR( IOLD, JOLD) +RINC 
AZ (I, J) =AZ (IOLD, JOLD) +ZINC 
WRITE (6,2005)  I,J,AR(I,J) ,AZ(I,J) 

CALL  MNIMX ( I , J) 

NCODE ( I , J) =1 
GO  TO  230 

220  CONTINUE 

IF(Il .GT. 12 . AND. IPTION.EQ. 7)  GO  TO  221 
IF(I1 .LT. 12 .AND. IPTION.EQ. 6)  GO  TO  221 
IOLD=I 
1=1+1 INC 

AR ( I , J) =AR( IOLD, J) +RINC 
AZ ( I , J) =AZ ( IOLD, J) +ZINC 
WRITE(6, 2005)  I,J,AR(I,J) ,AZ(I,J) 

NCODE ( I, J) =1 
CALL  MNIMX (I, J) 

JOLD= J 
J= J+JINC 

AR ( I , J) =AR ( I , JOLD) +RINC 
AZ ( I , J) =AZ ( I , JOLD) +ZINC 
NCODE ( I, J) =1 

WRITE(5 , 2005)  I,J,AR(I,J) ,AZ(I,J) 

CALL  MNIMX ( I, J) 

GO  TO  230 

221  JOLD= J 
J= J+JINC 

AR ( I , J ) =AR(I, JOLD) +RINC 
AZ (I, J) =AZ (I, JOLD) +ZINC 
NCODE ( I, J) =1 

WRITE(6, 2005)  I,J,AR(I,J) ,AZ(I,J) 

CALL  MNIMX ( I, J) 

IOLD=I 
1=1+1 INC 

AR ( I , J) =AR ( IOLD , J)  +RINC 
AZ (I, J) =AZ (IOLD, J) +ZINC 
NCODE ( I, J) =1 

WRITE(6, 2005)  I,J,AR(I,J) ,AZ(I,J) 

CALL  MNIMX ( I, J) 

230  CONTINUE 

IF (KAPPA. EQ. 1)  GO  TO  150 

IF(Il.GT. 12. AND. IPTION.EQ. 7)  GO  TO  231 

IF(Il.LT. I2.AND. IPTION.EQ. 6)  GO  TO  231 

IOLD=I 

1=1+1 INC 

AR ( I , J) =AR ( IOLD , J) +RINC 
AZ (I, J) =AZ ( IOLD, J)  +ZINC 
GO  TO  232 
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231  CONTINUE 
JOLD=J 
J= J+JINC 

AR  ( I ,  J) =AR ( I , JOLD) +RINC 
AZ (I, J) =AZ (I, JOLD) +ZINC 

232  CONTINUE 
NCODE ( I ,  J) =1 

WRITE(6, 2005)  I,J,AR(I,J) ,AZ(I,J) 

CALL  MNIMX (I, J) 

GO  TO  150 

C*  *********************************** 
C  GENERATE  CIRCULAR  ARCS  ON  BOUNDARY 

C*  *********************************** 
300  AR(I2,J2)=R2 
AZ (12, J2) =Z2 
NCODE ( 12, J2)  =  1 
CALL  MNIMX ( 12, J2) 

IF ( IPTION. EQ . 5)  GO  TO  320 

FIND  CENTER  OF  CIRCLE 

AR ( 13 , J3 ) =R3 
AZ (13, J3) =Z3 
NCODE ( 13, J3) =1 
CALL  MNIMX ( 13, J3) 

SLAC= ( Z2-Z1) / (R2-R1) 

SLBF=-1 . /SLAC 
SLCE= ( Z3-Z2) / (R3-R2) 

SLDF=-1 . / SLCE 

CHECK  FOR  INPUT  ERROR 

IF(ABS(SLAC-SLCE) .GT. .001)  GO  TO  310 
WRITE (6 , 2006)  Rl,Zl,R2,Z2,R3,Z3, SLAC , SLCE 
GO  TO  150 

310  R4=Rl+ (R2-R1) /2 . 

Z4=Zl+(Z2-Zl)/2. 

R5=R2+(R3-R2)/2. 

Z5=Z2+(Z3-Z2)/2. 

BBF=Z4-SLBF*R4 
BDF=Z5-SLDF*R5 
RC= (BBF-BDF) / (SLDF-SLBF) 

ZC=SLBF*RC+BBF 
WRITE(6, 2007)  RC , ZC 
KAPPA=1 
GO  TO  330 
320  KAPPA=2 
RC=R3 
ZC=Z3 

330  ISTRT=Il 
ISTP=I2 
JSTRT=J1 
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JSTP=J2 

RSTRT=Rl 

RSTP=R2 

ZSTRT=Z1 

ZSTP=Z2 

340  CALL  ANGLE (RSTRT, ZSTRT, RC, ZC, ANGl) 

CALL  ANGLE ( RSTP, ZSTP , RC, ZC , ANG2 ) 

IF ( ANG2 . LE. ANGl)  ANG2=2 . 0*PI+ANG2 

FIND  ANGULAR  INCREMENT 

DI=  ABS ( FLOAT ( ISTP-ISTRT) ) 

DJ=  ABS (FLOAT (JSTP-JSTRT)  ) 

I INC=0 
JINC=0 

IF(ISTRT.NE. IS TP)  IINC= ( ISTP-ISTRT) /IABS ( ISTP-ISTRT) 
IF ( JSTRT. NE . JSTP)  JINC= ( JSTP-JSTRT) /IABS ( JSTP-JSTRT) 
LAMDA=1 

IF( IINC.NE. 0. AND. JINC.NE. 0)  LAMDA=2 
DIFF=MAXl (DI,DJ) 

ITER=DIFF-1 . 

I F ( L AMDA . EQ . 2 )  DIFF=2 . *DIFF 

DELPHI= ( ANG2-ANG1 ) /DIFF 

WRITE (6,2008)  ANGl , ANG2 , DIFF, DELPHI 

CHECK  FOR  INPUT  ERROR 

IF ( LAMDA. NE. 2 . OR. DI . EQ . DJ)  GO  TO  350 
WRITE (6,2003) 

GO  TO  150 
350  1 0= I S TRT 
JO=JSTRT 
WRITE(6, 2004) 

INTERPOLATE 

NPT=IABS (12-11) +IABS ( J2- Jl) -1 
DO  380  M=1 , ITER 

359  IF ( LAMDA. EQ. 2)  GO  TO  360 
1=10+1 INC 

J= JO+JINC 
CALL  MNIMX ( I , J) 

NCODE ( I , J) =1 

CALL  CIRCLE (ANGl, DELPHI, RSTRT, ZSTRT, RC, ZC, I , J) 
WRITE(6, 2005)  I,J,AR(I,J) ,AZ(I,J) 

GO  TO  370 

360  1=10+1 INC 
J=JO 

NCODE ( I, J) =1 
CALL  MNIMX ( I, J) 

CALL  Cl RCLE ( ANG 1 , DE LPH I , RS TRT , ZSTRT , RC , ZC , I , J ) 
WRITE(6, 2005)  I,J,AR(I,J) ,AZ(I,J) 
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J= JO+JINC 
NCODE ( I ,  J) =1 
CALL  MNIMX (I, J) 

CALL  C I RCLE ( ANG1 , DELPH I , RSTRT , ZSTRT ,  RC , ZC , I , J ) 
WRITE (6,2005)  I , J, AR( I , J) , AZ ( I , J) 

370  IO=I 
380  JO=J 

I F ( L AMDA . NE . 2 )  GO  TO  390 
1=10+1 INC 
NCODE ( I, J) =1 
CALL  MNIMX ( I, J) 

CALL  CIRCLE ( ANG1 , DELPH I , RSTRT , Z STRT ,  RC ,  ZC ,  I ,  J ) 
WRITE (6,2005)  I, J, AR( I, J) ,AZ ( I, J) 

390  IF ( KAPPA. EQ . 2 )  GO  TO  150 
ISTRT=I2 
ISTP=I3 
JSTRT=J2 
JSTP=J3 
RSTRT=R2 
RS TP=R3 
ZSTRT=Z2 
ZSTP=Z3 
KAPPA=2 

399  GO  TO  340 

C*  ************************* 

C  CALCULATE  COORDINATES  OF  INTERIOR  POINTS 

C*************************^ 

400  IF (MAXJ. LE. 2)  GO  TO  430 
J2=MAXJ-1 

DO  420  N=1 ,500 
RESID=0 . 

DO  410  J=2 , J2 
Il=IMIN ( J) +1 
I2=IMAX ( J) -1 
DO  410  1=11,12 

IF(NCODE(I,J) .EQ.l)  GO  TO  410 

DR= ( AR ( 1+1 , J) +AR( 1-1 , J) +AR( I, J+l) +AR( I, J-l) ) /4 . 
DZ= ( AZ ( 1+1 , J) +AZ ( 1-1 , J) +AZ ( I , J+l ) +AZ ( I , J-l) ) /4 . 
RESID=RESID+ABS  (DR)  +ABS(DZ) 

AR(I, J) =AR ( I , J) +1 . 8*DR 
AZ (I, J) =AZ ( I , J) +1 . 8*DZ 
410  CONTINUE 

IF(N.EQ.l)  RES1=RESID 
IF(N. EQ.l. AND. RESID. EQ. 0. ) GO  TO  430 
IF ( RESID/RES1 . LT. 1 . E-5)  GO  TO  430 
420  CONTINUE 
430  WRITE (6,2009)  N 

C*  ************************* 
CALL  POINTS 

C*  ************************* 

1000  FORMAT  (515) 

1001  FORMAT  (3(213, 2F8.3) ,15) 


********** 

********** 


-AR(I, J) 
-AZ(I, J) 


********** 

********** 
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2000  FORMAT  (30H1  MESH  GENERATION  INFORMATION// 

1  41H0  MAXIMUM  VALUE  OF  I  IN  THE  MESH - 13/ 

2  41H0  MAXIMUM  VALUE  OF  J  IN  THE  MESH - 13/ 

3  41H0  NUMBER  OF  LINE  SEGMENT  CARDS - 13/ 

4  41H0  NUMBER  OF  BOUNDARY  CONDITION  CARDS - 13/ 

5  41H0  NUMBER  OF  MATERIAL  BLOCK  CARDS - 13///) 

2001  FORMAT  (//88H  INPUT  II  J1  Rl  Z1  12  J2  R2  Z 

12  13  J3  R3  Z3  IPTION/8X , 3 ( 214 , 2F8 . 4) ,16) 

2002  FORMAT  ( 5H  DI=F4.0,5H  DJ=F4.0,7H  DIFF=F4.0,7H  RINC=F8.3,7H  ZI 

1NC=F8 . 3 , 7H  ITER=I3,7H  IINC=I3,7H  JINC=I3,8H  KAPPA=I1) 

2003  FORMAT ( IX , 38H* *BAD  INPUT — THIS  LINE  IS  NOT  DIAGONAL) 

2004  FORMAT  (30H  I  J  AR  AZ) 

2005  FORMAT  (2l5,2F11.6) 

2006  FORMAT  (51H  **  BAD  INPUT  -  THESE  POINTS  DO  NOT  DEFINE  A  CIRCLE,/, 
13X,6F12.4,10X, 2E2  0.8) 

2007  FORMAT ( 19H  CENTER  COORDINATE, (Fll . 6 , IX , Fll . 6 , IX) ) 

2008  FORMAT  (7H  ANG1=F9.6,7H  ANG2=F9.6,7H  DIFF=F3.0,9H  DELPHI=F9.6) 

2009  FORMAT  (//30H  COORDINATES  CALCULATED  AFTER  13 , llH  ITERATIONS) 
RETURN 

END 

SUBROUTINE  MNIMX(I,J) 

COMMON/TD/IMIN ( 100) , IMAX (100) ,JMIN(25) ,JMAX(25) ,MAXI , MAXJ , NMTL, NBC 

IF(J.LT. JMIN(I) )  JMIN ( I) =J 

IF( J.GT. JMAX (I) )  JMAX ( I ) = J 

IF ( I . LT. IMIN ( J) )  IMIN ( J) =1 

IF(I.GT. IMAX(J) )  IMAX ( J) =1 

RETURN 

END 

SUBROUTINE  ADINA ( NFACE , FTHETA) 

INTEGER  CODE 

COMMON/NPDATA/R (1000)  , CODE (1000)  ,XR(1000)  ,Z(1000)  ,XZ(1000)  , 

1NPNUM (25,80) ,T(1000) ,XT(1000) 

COMMON/ARG/RRR ( 5 ) ,ZZZ(5) ,RR(4) ,ZZ(4) ,S(15,15) ,P(15) ,TT(6) , 

1H ( 6 , 15) ,CRZ (6,6) ,XI(10) , ANGLE (4) ,SIG(18) ,EPS(18) ,N 
COMMON/ELDATA/BETA ( 1000) ,EPR(1000) ,PR(200) ,SH(200) ,IX(1000,5) , 

IIP ( 200)  , JP ( 200)  , IS  ( 200)  , JS (200)  ,ALPHA(1000)  ,IT(200)  ,JT(200)  , 

2ST (200) , INP (1000) ,JNP(1000) 

COMMON/BASIC/ACELZ , ANGVEL , ANGACC , TREF , VOL , NUMNP , NUMEL , NUMPC , NUMSC , 
lNUMST 

COMMON/NXMESH/THETAN ( 4 ) ,NST(4) ,NUMS(4,5) ,NPC(8,8) 

DIMENSION  NPG (25,80,9) , FTHETA ( 9) ,NCOD(8) ,NEL(8) ,ID(6) 
COMMON/ANS1/NUMELS ( 4) ,NUMNPS(4) 

COMMON/NXDATA/NTP , NTYPS, NTS , NTOTS 
COMMON/NON AXI/S1 (30,30) ,P1(30) , THETA, BS1 ( 6 , 30) 

COMMON/SOL VE/X (4428) ,Y(4428) ,TEM(4428) , NUMTC , MB AND 

COMMON/TD/IMIN (100) ,IMAX(100) ,JMIN(25) ,JMAX(25) , MAX I , MAXJ, NMTL, NBC 
COMMON / CONVRG/ I DONE 
COMMON/P LANE/NPP 

COMMON/RE SUL T/BS( 6, 15) ,D(6,6) ,C(6,6) ,AR,BB(6,9) ,CNS(6,6) 
COMMON/MATP/RO (6) ,E(12,16,6) ,EE(16) ,AOFTS(6) 

COMMON/TOT AL/NELN 
DATA  NPG/18000*0/ 
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NI=NST ( NTP) 

C  LOOP  ON  #  OF  SEGMENTS  OF  TYPE  NTP 
DO  500  1=1, NI 
DO  200  J=1 ,NUMNP 
NFL=NUMS (NTP, I) 

NFR=NFL+1 
IPl=INP( J) 

JPl=JNP ( J) 

NF=NFL 
DO  200  K=1 , 2 
NC=CODE ( J) +1 

IF (NPG ( IPl , JP1 ,NF) .NE.0)GO  TO  200 
NPG ( IPl , JPl , NF) =IP1+(JP1-1) *25+NF*200 
DO  10  11=1,8 
10  ID  ( 1 1)  =0 

IF (NC . EQ . 3 . OR. NC . EQ. 4 . OR. NC . EQ. 7 . OR. NC . EQ. 8)  ID(1) =1 

IF(NF.EQ.l)  ID ( 3) =1 
IF ( NF . EQ. NTOTS+1)  ID(2)=1 
IF ( R ( J) . EQ .0.0)  ID (2) =1 
IF(R(J) .EQ.0.0)  ID ( 3) =1 

WRITE  NODAL  DATA  ON  FILE  2 

WRITE ( 2 , 100) NPG (IPl , JPl ,NF) , (ID(KK) ,KK=1,6) ,Z(J) ,R(J) 
100  FORMAT ("X" ,I4,1X,I4,5I5, 3F10 . 4 , 4X , "0 " ) 

200  NF=NFR 

DO  350  J=1 , NUMEL 
NF=NUMS (NTP, I) 

I1=IX(J,1) 

I2=IX ( J , 2) 

I3=IX ( J, 3) 

I4=IX ( J, 4) 

IPl=INP ( II) 

JP1=JNP (II) 

NEL(l) =NPG(IP1, JPl, NFL) 

IF ( IPl . EQ. l)NEL(l) =NPG ( IPl , JPl , 1) 

NEL  (  4)  =NPG  ( IPl ,  JPl ,  NFR) 

IPl=INP ( 12 ) 

JP1=JNP ( 12 ) 

NEL (2) =NPG ( IPl , JPl , NFL) 

NEL  (  3)  =NPG  ( IPl ,  JPl , NFR) 

IP1=INP ( 14) 

JPl=JNP ( 14 ) 

NEL ( 5 ) =NPG ( IPl , JPl , NFL) 

IF ( IPl . EQ. 1) NEL ( 5) =NPG ( IPl , JPl , 1) 

NEL ( 8 ) =NPG ( IPl , JPl ,NFR) 

IPl=INP ( 13 ) 

JP1=JNP ( 13 ) 

NEL (6) =NPG ( IPl , JPl ,NFL) 

NEL  ( 7 )  =NPG  ( IPl ,  JPl ,  NFR) 

NELN=NELN+1 
NZ=0 


, FTHETA(NF) 
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N8=8 

WRITE (3,300) NELN, N8  ,  N8  ,  NZ  , IX ( J , 5)  ,NELN,NZ ,NZ 

300  FORMAT (8 15 , F10 . 0) 

WRITE (3,301)  (NEL(K) ,K=1, 8) 

301  FORMAT (815) 

WRITE  MATERIAL  AXES  DATA  ON  FILE  4 

NPI=NEL ( 1) +7000 
NP J=NEL ( 1) +8000 
NPK=NEL ( 1) +9000 
WRITE (4,302)  NELN, NPI , NPJ , NPK 

302  FORMAT ( 4 15 ) 

WRITE  NODAL  POINTS  FOR  MATERIAL  AXES 

FTRAD=FTHETA ( NF) /57. 29578 
CA=COS (ALPHA (J) ) 

CB=COS (BETA ( J) ) 

SA=SIN (ALPHA( J)  ) 

SB=SIN (BETA ( J) ) 

CF=COS (FTRAD) 

SF=SIN ( FTRAD) 

XPI=0 . 0 
YPI=0 . 0 
ZPI=0.0 
XPJ=CA*CB 

YPJ=CA*SB*CF+SA*S? 

ZPJ=CA*SB*SF-SA*CF 

XPK=SA*CB 

YPK=SA*SB*CF-CA*SF 

ZPK=SA*SB*SF+CA*CF 

WRITE ( 7 , 30  3 ) NPI , NZ , NZ , NZ , NZ , NZ , NZ , XP I , YP I , Z  P I 

WRITE (7,303) NP J, NZ , NZ , NZ , NZ , NZ , NZ , XP J, YPJ, ZPJ 
WRITE ( 7 , 3  0  3 ) NPK , NZ , NZ , NZ , NZ , NZ , NZ , XPK , YPK , Z  PK 

303  FORMAT (1X,I4,1X,I4,5I5,3F10.4,4X, "0") 

350  CONTINUE 

DO  500  L=1,NUMPC 
IP1=INP (IP (L) ) 

JP1=JNP ( IP ( L)  ) 

IP2  =  INP ( JP (L)  )  I 

JP2=JNP ( JP (L)  ) 

NEL(l) =NPG ( IPl , JPl , NFL) 

NEL (2) =NPG ( IP2 , JP2 ,NFL) 

NEL ( 3) =NPG(IP2, JP2,NFR) 

NEL (4) =NPG ( IPl , JPl ,NFR) 

NCUR=1 

WRITE (1, 400) NCUR, (NEL ( K) ,K=1,4) 

WRITE (1,401)  PR ( L) ,PR(L) ,PR(L) ,PR(L) 

400  FORMAT ( 9 15 ) 

401  FORMAT ( 5F10 . 0 , 15) 

500  CONTINUE 


38 


o  o 


END 

FUNCTION  NODE ( I ,  J) 

COMMON/TD/IMIN ( 100) , IMAX (100) ,JMIN(25) ,JMAX(25) ,MAXI , MAXJ, NMTL, NBC 
NODE=0 

DO  100  JJ=1,J 
NSTART=IMIN ( JJ) 

NSTOP=IMAX (JJ) 

DO  100  II=NSTART, NSTOP 
NODE=NODE+l 

IF(JJ.EQ. J. AND. II.EQ. I)  RETURN 
100  CONTINUE 
RETURN 
END 

SUBROUTINE  POINTS 
INTEGER  CODE 

COMMON /BASIC/ACELZ , ANGVEL, ANGACC , TREF, VOL , NUMNP , NUMEL, NUMPC , NUMSC , 
lNUMST 

COMMON/MATP/RO ( 6 ) ,E(12,16,6) ,EE(16) ,AOFTS(6) 

COMMON/NPDATA/R (1000) ,CODE(1000) ,XR(1000) ,Z(1000) ,XZ(1000) , 

1NPNUM ( 25 , 80) ,T(1000) ,XT(1000) 

COMMON/E LDATA/BETA (1000) ,EPR(1000) ,PR(200) ,SH(200) ,IX(1000,5) 

IIP ( 200)  , JP (200)  , IS  ( 200)  , JS (200)  ,ALPHA(1000)  ,IT(200)  ,JT(200)  , 
2ST(200) ,INP(1000) ,JNP(1000) 

COMMON/SOL VE/X (4428) ,Y(4428) ,TEM(4428) , NUMTC , MBAND 
COMMON/TD/IMIN ( 100) ,IMAX(100) ,JMIN(25) ,JMAX(25) , MAXI , MAXJ , NMTL, NBC 
COMMON/P LANE/NP P 

iJIMENSLON  AR (25,80) ,AZ(25,80) ,MATRIL ( 100 , 5 ) , BLKANG (100) , BLKALF ( 1 

DIMENSION  IBNG(100) ,NBNG(100) 

EQUIVALENCE  ( R  ( 1)  ,  AR)  ,  ( Z  ( 1)  ,  AZ) 

ESTABLISH  NODAL  POINT  INFORMATION 

************************************ 
NEL=0 
NODSUM=0 
DO  100  J=1,MAXJ 
NSTART=IMIN ( J) 

NSTOP=IMAX ( J) 

DO  100  I=NSTART, NSTOP 
100  NODSUM=NODSUM+l 
NELSUM=0 
JJMAX=MAXJ-1 
DO  110  JJ=1 , JJMAX 

NSTOP=MIN0 ( IMAX ( JJ) , IMAX ( JJ+1) ) -1 
NSTART=MAX0 ( IMIN ( JJ) , IMIN ( JJ+1) ) 

DO  110  I I=NSTART , NSTOP 
110  NELSUM=NELSUM+1 
NUMNP =NODSUM 
NUMEL=NELSUM 
DO  120  J=1 ,MAXJ 
NSTART=IMIN ( J) 

NSTOP=IMAX (J) 

DO  120  I=NSTART, NSTOP 
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NPNUM ( I , J) =NODE ( I , J) 

NP=NPNUM ( I , J) 

INP (NP) =1 
JNP(NP) =J 
R (NP) =AR ( I / J) 

120  Z(NP) =AZ (I, J) 

************************************ 
READ  AND  ASSIGN  BOUNDARY  CONDITIONS 

************************************ 
INITIALIZE 

************************************ 
DO  130  1=1 ,NUMNP 
CODE (I) =0 

IF  (R  ( I)  .  EQ .  0  .  .  AND .  NPP .  EQ .  0 )  CODE(I)=l. 

XR ( I ) =0. 

XZ(I) =0. 

XT ( I) =0.0 
130  T ( I) =0. 

IF (NBC. EQ. 0)  GO  TO  210 
DO  200  IBCON=l,NBC 

READ (5 , 1002)  II , 12 , Jl , J2 , ICN, RCON, ZCON,TCON 
DO  200  1=11,12 
DO  200  J=J1,J2 
NP=NPNUM ( I , J) 

CODE (NP) =ICN 
XR (NP) =RCON 
XT (NP) =TCON 
200  XZ (NP) =ZCON 
210  MPRINT=0 

DO  230  J=1 ,MAXJ 
NSTART=IMIN ( J) 

NSTOP=IMAX ( J) 

DO  230  I=NSTART, NSTOP 
NP=NPNUM ( I , J) 

IF (MPRINT. NE. 0)  GO  TO  220 
WRITE(6, 2000) 

MPRINT=59 

220  MPRINT=MPRINT-1 

230  WRITE(6, 2001)  I , J,NP,CODE (NP) ,R(NP) ,Z(NP) ,XR(NP) ,XZ(NP) ,XT(NP) 

c*  *********************************** 

C  ASSIGN  MATERIALS  IN  BLOCKS 

C*  *********************************** 
DO  300  Ml=l ,NUMEL 
300  IX  (Ml ,  5)  =0 

DO  310  IMTL=1 ,NMTL 

READ  (5,1000)  MTL, (MATRIL ( IMTL, IM) ,IM=2,5) , BLKANG ( IMTL) , BLKALF ( IMT 
1L)  ,  IBNG  ( IMTL)  ,  NBNG  ( IMTL) 

310  MATRIL ( IMTL, 1) =MTL 

c*  *********************************** 

C  ESTABLISH  ELEMENT  INFORMATION 

C*  *********************************** 
J JMAX=MAXJ-1 
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N=0 

MTL=1 

KTL=1 

DO  440  JJ-1 , JJMAX 

NSTOP=MINO (IMAX(JJ)  ,IMAX(JJ+1)  ) -1 
NSTART=MAX0 ( IMIN ( J J) , IMIN(JJ+1) ) 

DO  440  II=NSTART, N3TOP 
NEL=NEL+1 

DO  400  IMTL=1 , NMTL 

IF ( II . LT.MATRIL ( IMTL, 2) )  GO  TO  400 
IF ( I I . GE . MATRIL ( IMTL, 3) )  GO  TO  400 
IF ( JJ. LT.MATRIL ( IMTL, 4) )  GO  TO  400 
IF (J J. GE . MATRIL ( IMTL, 5) )  GO  TO  400 
KAT=IMTL 

MAT=MATRIL ( IMTL, 1) 

400  CONTINUE 

IF(KAT.EQ.KTL)  GO  TO  410 

KTL=KAT 

MTL=MAT 

CO  TO  420 

410  IF ( II . EQ . NSTART)  GO  TO  420 

IF( JJ.NE. JJMAX.OR. II.NE.NSTOP)  GO  TO  440 
M=NEL+1 
IANG=ICNG 
NANG=NCNG 
GO  TO  421 
420  I=NPNUM (II, JJ) 

J=I+1 

K=NPNUM( II+l , JJ+1) 

L=K-1 

M=NEL 

IX (M, 1) =1 

IX (M, 2) =J 

IX (M, 3) =K 

IX (M, 4) =L 

IX  (M,  5)  =MTL 

BETA (M) =BLKANG (KTL) 

ALPHA (M) =BLKALF(KTL) 

IANG=ICNG 
NANG=NCNG 
ICNG=IBNG (KTL) 

NCNG=NBNG (KTL) 

421  NC=2 
430  N=N+1 

IF(M.LE.N)  GO  TO  440 
IX (N, 1) =IX (N-l , 1) +1 
IX(N,2)=IX  (N-l ,  2)  +1 
IX (N, 3) =IX (N-l , 3) +1 
IX (N , 4) =IX (N-l , 4) +1 
IX (N , 5) =IX  (N-l , 5) 

BETA (N) =BETA (N-l ) 

IF ( IANG. EQ. 1)  GO  TO  442 
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ALPHA (N) =ALPHA (N-l ) 

GO  TO  443 

442  CONTINUE 

IF (NC.GT.NANG)  GO  TO  444 
ALPHA (N) =ALPHA (N-l) 

GO  TO  443 
444  NC=1 

ALPHA ( N) =-ALPHA ( N-l ) 

443  CONTINUE 
NC=NC+1 

IF(M.GT.N)  GO  TO  430 
440  CONTINUE 

IF(NUMNP.GT. 2000)  WRITE ( 6 , 2002) 

A**************************,^,^^^^^^^^ 

SET  NODAL  POINT  TEMPERATURE  TO  REFERENCE  TEMPERATURE 
a***************************^^*^^^^* 

I F ( NUMTC . NE . 0 )  RETURN 
DO  500  N=1 , NUMNP 
500  T (N) =TREF 

1000  FORMAT  (5l5,2F10 .0,215) 

1002  FORMAT (415, 110, 3F10 .0) 

2000  FORMAT  (128H1  I  J  NP  TYPE  R-ORDINATE  Z-ORDINA 

1TE  R  LOAD  OR  DISPLACEMENT  Z  LOAD  OR  DISPLACEMENT  T  LOAD  OR  DISP 
2LACEMENT) 

2001  FORMAT  ( 2 15 , 16 , 112 , F13 . 6 , F14 . 6 , E2 6 . 7 , E2 4 . 7 , E24 . 7) 

2002  FORMAT  (35H  BAD  INPUT  -  TOO  MANY  NODAL  POINTS) 

RETURN 
END 

SUBROUTINE  TEMP(R,Z,T) 

COMMON/SOL VE/X (4428) ,Y(4428) ,TEM(4428) , NUMTC, MBAND 
DIMENSION  SMALL ( 20 ), ISM (20) 

*  *********************************** 
INITIALIZE 

C*  *********************************** 
J=1 

JMAX=16 

IF (NUMTC. LT.JMAX)  JMAX=NUMTC 
DO  10  1=1 , JMAX 
SMALL (I) =0 . 

10  ISM ( I) =0 

C*  *********************************** 

C  FIND  THE  JMAX  CLOSEST  POINTS 

C*  *********************************** 
DO  50  1=1, NUMTC 
DSQ=  (X  ( I)  -R)  **2+(Y(I)-Z)  **2 
IF(DSQ.GT. . IE-4)  GO  TO  20 
T=TEM ( I) 

RETURN 

20  IF(I.EQ.l)  SMALL ( 1 ) =DSQ 
IF(I.EQ.l)  ISM (1) =1 
IF(I.EQ.l)  GO  TO  50 

IF ( SMALL ( J) .LE.DSQ.AND.J. LT.JMAX)  SMALL (J+l) =DSQ 
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IF ( SMALL ( J) .LE.DSQ.AND. J.LT. JMAX)  ISM(J+1)=I 
IF (SMALL (J) .LE.DSQ)  GO  TO  40 
DO  30  K=1,J 
JB=J-K  +1 

IF(JB.EQ.O)  GO  TO  40 
SMALL (JB+1) =SMALL ( JB) 

ISM ( JB+1) =ISM ( JB) 

SMALL (JB)  =DSQ 
ISM( JB) =1 

IF(JB.EQ.l)  GO  TO  40 
IF (SMALL (JB-1) .LE.DSQ)  GO  TO  40 
30  CONTINUE 

40  IF( J. LT. JMAX)  J=J+1 
50  CONTINUE 

*********************** 

FIND  THE  THIRD  TEMPERATURE  POINT  BY  AREA 
*********************** 

JCHK= JMAX-2 
J=0 

Il=ISM ( 1) 

I2=ISM (2) 

60  I3=ISM ( J+3 ) 

ARE A= . 5  0  * ( Y ( 1 1 ) *X(I3)-Y(I3) *X(Il) +Y(I3) *X(I2) -Y(I2) *X(I3)  + 

1  Y(I2) *X(I1) -Y(I1) *X(I2) ) 

Dl= (X  (12) -X(I1) ) *  *2+ ( Y ( 12 ) -Y(I1) ) **2 

IF  Dl  IS  APPROXIMATELY  0.  IT  IS  ASSUMED  THAT  THERE  EXISTS  A 
DUPLICATION  OF  INPUT 
IF (Dl .GT. . IE-3)  GO  TO  70 
12  =  13 


************* 

TEST 

************* 


J=J+1 


GO  TO  60 

70  IF ( AREA**2 .GT. . l*Dl *SMALL ( 1) )  GO  TO  80 
J=J+1 


IF(J.LT. JCHK)  GO  TO  60 
WRITE (6,2000)  II, 12, 13, J 
T=TEM ( 11 ) 

RETURN 


c*  *********************************** 
C  FIND  TEMPERATURE  INTERCEPT 

C*  *********************************** 
80  DETA=Y ( II) * (TEM ( 13 ) -TEM ( 12 ) ) +Y ( 12 ) * (TEM ( II) -TEM ( 13 ) ) 

1  +Y (13) * (TEM (12) -TEM( II) ) 

DETB=X ( II ) * (TEM (12) -TEM(I3) ) +X(I2) * (TEM ( 13 ) -TEM ( II) ) 

1  +X  (13)  *  (TEM  (II)  -TEM  (12)  ) 


DETC=TEM ( 11) * (X ( 12 ) *Y ( 13) -X ( 13) *Y ( 12) )+TEM(l2) 

1Y  (13)  )  +TEM  ( 13 )  *(X(II)  *Y(I2)-X(I2)  *Y(II)  ) 

T= (DETA*R+DETB*Z+DETC) / ( 2 . *AREA) 

2000  FORMAT  (28H  ERROR  IN  TEMPERATURE  INPUT, 5H  11=14. 5H 
15H  13=14, 4H  J=I4) 


(X ( 13 ) *Y ( II) -X (II) 


12  =  14 


RETURN 

END 


SUBROUTINE  TEM2 (NUMNP) 
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INTEGER  CODE 

COMMON/NPDATA/R( 1000) ,CODE(1000) ,XR(1000) ,Z(1000) ,XZ(1000) , 
1NPNUM (25,80) ,T(1000) ,XT(1000) 

READ (5, 1000)  TCONST 
DO  100  N=1,NUMNP 
100  T (N) =TCONST 
1000  FORMAT ( F10 . 0) 

RETURN 

END 

END 

N 
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Listing  of  the  NASGEN  Program 
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PROGRAM  NASGEN ( INPUT , OUTPUT, NAS, TAPE5=INPUT , TAPE6=OUTPUT, TAPE1= 

1  NAS) 

INTEGER  CODE 

COMMON/NPDATA/R ( 1000) ,CODE(1000) ,XR(1000) ,Z(1000) ,XZ(1000) , 
1NPNUM(25,80) ,T(1000) ,XT(1000) 

COMMON/ARG/RRR  ( 5 )  ,ZZZ(5)  ,RR(4)  ,  ZZ  (  4)  ,S(15,.15)  ,P(15)  ,TT(6)  , 

1H ( 6 , 15) ,CRZ (6 , 6) ,XI(L0) , ANGLE (4) ,SIG(18) ,EPS(18) ,N 
COMMON/ELDATA/BETA( 1000) ,EPR(1000) ,PR(200) ,SH(200) ,IX(1000,5) , 

IIP ( 200)  , JP ( 200)  , IS ( 200)  , JS ( 200)  ,ALPHA(1000)  ,IT(200)  ,JT(200)  , 
2ST(200) ,INP(1000) ,JNP(1000) 

COMMON/BASIC/ACELZ , ANGVEL, ANGACC, TREF, VOL, NUMNP , NUMEL, NUMPC , NUMSC , 
1NUMST 

COMMON/NXMESH/THETAN (4) ,NST(4) ,NUMS(4,5) ,NPC(8,8) 
COMMON/ANS1/NUMELS (4) ,NUMNPS(4) 

COMMON/NXDATA/NTP , NTYPS , NTS , NTOTS 
COMMON/NONAXI/S1 (30, 30) ,P1(30) , THETA, BS1 ( 6 , 30 ) 

COMMON/SOL VE/X (4428) ,Y(4428) ,TEM(4428) , NUMTC , MB AND 

COMMON/TD/IMIN ( 100) ,IMAX(100) ,JMIN(25) ,JMAX(25) ,MAXI , MAXJ , NMTL, NBC 
COMMON/ CONVRG / I DONE 
C  OMMON/P  L ANE /NP  P 

COMMON/RE SULT/BS( 6, 15)  ,D(6,6)  ,C(6,6)  ,AR,BB(6,9)  ,CNS(6,6) 
COMMON/MATP/RO ( 6 )  ,E(12,16,6)  ,EE(16)  ,AOFTS(6) 

DIMENSION  TITLE ( 20) , FTHETA( 8) 

C*  *********************************** 
C  READ  AND  WRITE  CONTROL  INFORMATION 

C*  *********************************** 

READ (5, 3000)  NTYPS, NTOTS 
DO  150  1=1, NTYPS 

150  READ (5,3001)  THETAN ( I) , NST ( I) 

DO  151  1=1, NTYPS 
J2  =NST ( I) 

151  READ (5,3000)  (NUMS (I, J) , J=1 , J2) 

DO  152  1=  1, NTOTS 

152  READ (5 , 3002)  (NPC ( I , J) , J=1 , 8 ) 

3000  FORMAT(8I10) 

3001  FORMAT ( FI 0 . 5 , 1 1 0 ) 

3002  FORMAT (8110) 

NFACE=NTOTS+l 
FTHETA(l) =0.0 

DO  4100  1=1, NTOTS 
DO  4000  J=l, NTYPS 
K2=NST ( J) 

DO  4000  K=1,K2 

4000  I F ( NUMS ( J , K ) . EQ . I ) GO  TO  4100 
4100  FTHETA (1+1 ) =FTHETA( I) +THETAN ( J) 

WRITE (6,4200) ( (FTHETA ( I ) ,1) ,1=1, NTOTS) 

4200  FORMAT ( "  ANGLE  IN  DEGREES  ",5X,"  BEGINNING  SEGMENT  #"/ 

1  (E20 .10,110) ) 

WRITE(S, 3010) 

3010  FORMAT ( "1 ", "SEGMENT  DATA  FOR  NONAX I SYMMETRIC  PROBLEM") 

WRITE (6 , 3011)  NTYPS, NTOTS 

3011  FORMAT  ("  NUMBER  OF  TYPES  OF  SEGMENTS  =  ",I5,//, 
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"  NUMBER  OF  TOTAL  SEGMENTS  =",I5) 

DO  153  1=1 , NTYPS 
WRITE(6, 3012)  I , THETAN ( I) ,NST( I) 

3012  FORMAT ( "  ",///,"  SEGMENT  TYPE  =  " ,15/,”  THETA  =  ",F10.5/, 

1  "  NUMBER  OF  SEGMENTS  OF  THIS  TYPE  =  ",I5) 

J2  =  NST(I) 

WRITE  (6,3013)  (NUMS(I,.J)  ,J=1,J2) 

3013  FORMAT ("  "  SEGMENT  NUMBERS  IN  GLOBAL  SYSTEM  ARE  ".515) 

153  CONTINUE 

DO  154  1=1 , NTOTS 

154  WRITE (6,3014) I, (NPC(I,J) ,J=1,8) 

3014  FORMAT ( "  " , "CONNECTING  NODES  FOR  SEGMENT" , 15 , "  ARE",8l5) 

DO  950  NTP  =  1, NTYPS 

THETA=  THETAN (NTP)  /5 7. 295780 

50  READ ( 5 , 1000  ) TITLE, NNLA, NUMTC , NUMMAT, NUMPC ,NUMSC , NUMST, TREF 
1, INERT, NLINC, INCI, INCF,IPLOT 

WRITE ( 6 , 2  0  0  0 ) TITLE , NNLA, NUMTC , NUMMAT , NUMPC , NUMSC , NUMST , TREF , INERT , 
1NLINC 
NPP=0 

0*  *********************************** 
C  GENERATE  FINITE  ELEMENT  MESH 

C*  *********************************** 
100  CALL  MESH 

NUMELS (NTP)  =  NUMEL 
NUMNPS (NTP)  =  NUMNP 

C*  *********************************** 
C  READ  AND  WRITE  T5MPERATURE  DATA 

C*  *********************************** 
103  IF (NUMTC. EQ.0)  GO  TO  440 

IF (NUMTC. GT. 0)  READ(5,1001)  (X  ( I)  , Y ( I )  , TEM ( I)  , 1=1 , NUMTC) 

IF (NUMTC. EQ. -2)  CALL  TEM2 (NUMNP) 

IF (NUMTC. EQ. -2)  GO  TO  440 
MPRINT=0 

DO  210  1=1, NUMTC 
IF (MPRINT.NE. 0)  GO  TO  200 
WRITE(6, 2001) 

MPRINT=5  9 

200  MP R I NT=MP R I NT- 1 

210  WRITE  (6 , 2002)  X  ( I)  ,  Y  ( I)  , TEM  ( I) 

MPRINT=0 

DO  230  N=l, NUMNP 
IF (MPRINT. NE . 0)  GO  TO  220 
WRITE(6, 2003) 

MPRINT=5  9 

220  MPRINT=MPRINT-1 

CALL  TEMP (R(N) ,Z(N) ,T(N)) 

230  WRITE (6,2004)  N, R (N) , Z (N) , T (N) 

440  MPRINT=0 

DO  460  N=l, NUMEL 
IF (MPRINT.NE. 0)  GO  TO  450 
WRITE (6,2008) 

MPRINT=59 
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450  MPRINT=MPRINT-1 
II=IX(N,1) 

JJ=IX (N, 2) 

KK=IX (N , 3) 

LL=IX (N, 4) 

TEM  IS  TEMPORARY  STORAGE  FOR  ELEMENT  TEMPERATURES 

TEM  (N)  =  ( T  ( 1 1 )  +T  ( J J)  +T  ( KK)  +T  ( LL)  )  /4 . 0  0 
460  WRITE(6, 2009)  N,  (IX (N, I)  , 1  =  1 , 5)  ,BETA(N)  , ALPHA(N)  ,TEM(N) 

DO  470  K=1 , NUMEL 
470  T (K) =TEM (K) 

************************************ 
READ  AND  WRITE  MATERIAL  PROPERTIES 

(2*  **********************************  * 
500  CONTINUE 

DO  510  M=1 , NUMMAT 

READ (5, 1004)  MTYPE, (NT, RO (MTYPE) , AOFTS ( MTYPE) ) 

WRITE(6, 2010)  MTYPE, NT, RO (MTYPE) 

READ (5, 1005) ( (E ( I , J, MTYPE) ,J=1,14) ,1=1, NT) 

IF (AOFTS (MTYPE) .NE.l.)  WRITE ( 6 , 2011) ( (E ( I , J, MTYPE) ,J=1,13) ,1=1, NT) 
IF (AOFTS (MTYPE) .EQ.l.)  WRITE ( 6 , 2012) ( (E ( I , J, MTYPE) ,J=1,13) ,1=1, NT) 
DO  510  I=NT, 12 
DO  510  J=1 , 16 

510  E (I, J, MTYPE) =E(NT,J, MTYPE) 

DO  900  NL=1 ,NLINC 

WRITE (6 , 2030)  NL 
ACELZ=0 .00 
ANGVEL=0 .00 
ANGACC=0 .00 

IF ( INERT  .EQ.  0)  GO  TO  511 

IF (NL  .NE.  1  .AND.  INCI  . EQ .  0)  GO  TO  511 

Q**  ***************************************  ****************************** 

C  READ  AND  WRITE  DYNAMIC  FORCES 

(2*  *********************************************************************  * 
READ (5, 10 30)  ACELZ,  ANGVEL,  ANGACC 
WRITE (6, 2031)  ACELZ,  ANGVEL,  ANGACC 

511  CONTINUE 

Q*  *********************************** 

C  READ  AND  WRITE  PRESSURE  AND  SHEAR  BOUNDARY  CONDITIONS 

C*  *********************************** 
IF (NL  .NE.  1  .AND.  INCF  .EQ.  0)  GO  TO  700 
600  IF ( NUMPC . EQ . 0 )  GO  TO  630 
MPRINT=0 

DO  620  L=l, NUMPC 
IF (MPRINT. NE . 0)  GO  TO  610 
WRITE (6,2013) 

MPRINT=58 

610  MPRINT=MPRINT-1 

READ (5,1006)  IP (L) , JP (L) ,PR(L) 

620  WRITE(6, 2014)  IP ( L) , JP ( L) , PR ( L) 

630  I F ( NUMSC . EQ . 0 )  GO  TO  701 
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MPRINT=0 

DO  650  L=1 , NUMSC 
IF (MPRINT. NE. 0)  GO  TO  640 
WRITE (6 , 2015) 

MPRINT=58 

640  MPRINT=MPRINT-1 

READ  (5,1006)  IS  ( L)  ,  JS  ( L)  ,  SH  ( L) 

650  WRITE (6,2014)  IS (L) , JS (L) , SH (L) 

701  IF (NUMST. EQ . 0 )  GO  TO  700 
MPRINT=0 

DO  680  L=l, NUMST 
IF (MPRINT. NE. 0)  GO  TO  670 
WRITE(6, 2025) 

MPRINT=58 

670  MPRINT=MPRINT-1 

READ  (5,10  06)  IT  ( L)  ,  JT  ( L)  ,  ST  ( L) 

680  WRITE (6, 2014) IT ( L) ,JT(L) , ST ( L) 
************'k'k'k'k**ie**iei(**'fC'kje£ici(** 
DETERMINE  BANDWIDTH,  INITIALIZE  ELASTIC-PLASTIC  RATIO, 

AND  CONVERT  BETA  FROM  DEGREES  TO  RADIANS 
********************************* 
700  J=0 


DO  710  N=1 , NUMEL 
IX (N, 5) =IABS ( IX (N, 5) ) 

DO  710  1=1,4 
DO  710  L=1 , 4 
KK=IABS (IX(N,I)-IX(N,L) ) 
IF(KK.GE.J)  J=KK 
710  CONTINUE 


720 

721 


MBAND=3  * J+3 
IF (NL. GT. 1)  GO  TO  721 
DO  720  N=l, NUMEL 
EPR(N) =1 . 

ALPHA(N) =ALPHA(N) /57. 295780 
BETA (N) =BETA(N) /57. 295780 
CONTINUE 
900  CONTINUE 

CALL  NASTRN ( NFACE , FTHETA) 

950  CONTINUE 

WRITE (1,4500) ( (MTYPE, E (1,4, MTYPE) , E ( 1 , 10 , MTYPE) , E ( 1 , 7 , MTYPE) 
1  , RO (MTYPE) ) ,MTYPE=1 ,NUMMAT) 

4500  FORMAT ("MAT1",4X, 18, 4E8. 2) 

WRITE (1,4600) ACELZ 

4600  FORMAT ( "CORD1C", 9X, " 1", 5X, "101", 5X, "102", 5X,"103"/"GRID",9X, 

1  "1 01 ", 13X, "0. 0  ", 5X, "0. 0",5X,"0.0",10X, "123456 "/"GRID", 9X, 

2  "102" , 13X, "0 . 0" , 5X, "0 . 0" , 5X, "1 . 0 ", 10X , "123456 "/"GRID" , 9X , 

3  "103" , 13X , " 1 . 0 " , 5X , "0 . 0 " , 5X , " 1 . 0 " , 10X , "123456 "/"GRAV" , 11X 

4  , "2 " , 7X, "1 " ,E8 . 2 , 5X, "0 . 0" , 5X, "0 . 0 " , 4X, "-1 . 0"/"PSOLID" , 

5  9X , " 1 " , 7X , " 1 " , 7X , " 1 "/"ENDDATA" ) 

910  CONTINUE 


1000  FORMAT ( 20A4/6 15 , F5 . 0 , 515 ) 

1001  FORMAT ( 3F10 . 0) 


* 


* 


* 


* 


* 


* 
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1004  FORMAT  (2I5,2f10.0) 

1005  FORMAT (7F10 . 0) 

1006  FORMAT  (2I5,F10.0) 

1030  FORMAT ( 3F10 . 0) 

2000  FORMAT  (2H1  ,20A4/ 

1  33H0  NUMBER  OF  APPROXIMATIONS - 14/ 

2  33H0  NUMBER  OF  TEMPERATURE  CARDS - 14/ 

3  3  3H0  NUMBER  OF  MATERIALS - 14/ 

4  3  3H0  NUMBER  OF  PRESSURE  CARDS - 14/ 

5  33H0  NUMBER  OF  SHEAR  CARDS - 14/ 

6  33H0  NUMBER  OF  TORSION  CARDS - 14/ 

7  33H0  REFERENCE  TEMPERATURE - E12.4/ 

8  3 3H0  NUMBER  OF  INERTIA  CARDS - 14/ 

9  33H0  NUMBER  OF  LOAD  INCREMENTS - 14/) 

2001  FORMAT  ( lHl , 13X , lHR, 14X , 1HZ , 14X , 1HT) 

2002  FORMAT  (3F15.3) 

2003  FORMAT  (35H1  N  R  Z  T) 

2004  FORMAT  ( 15 , 2F10 . 4 ,F10 . 0) 


2008  FORMAT  (74H1  EL  I  J  K  L  MATERIAL  ANGLE  BETA 
1LPHA  TEMPERATURE) 

2009  FORMAT  ( 15 , 414 , 18 , Fll . 1, 2F13 . 3) 

2010  FORMAT  ( lHl , "MATERIAL  IDENTIFICATION  NUMBER  =",12/ 
llH  , "NO.  OF  MATERIAL  TEMPERATURE  CARDS  =",12/ 

2 1H  , "MASS  DENSITY  =",E15.7) 

2011  FORMAT  ( 1H  , "TEMPERATURE  =",E15.7/ 

llH  , "MODULUS  OF  ELASTICITY-EN  =",E15.7/ 

2 lH  , "MODULUS  OF  ELASTICITY-ES  =",E15.7/ 

3 1H  , "MODULUS  OF  ELASTICITY-ET  =",E15.7/ 

41H  , "POISSON  RATIO-NUNS  =",E15.7/ 

51H  , "POISSON  RATIO-NUNT  =" ,E15.7/ 

6 lH  , "POISSON  RATIO-NUST  =",E15.7/ 

71H  , "SHEAR  MODULUS-GNS  =",E15.7/ 

81H  , "SHEAR  MODULUS-GST  =",E15.7/ 

9lH  , "SHEAR  MODULUS-GTN  =",E15.7/ 

llH  , "COEFFICIENT  OF  THERMAL  EXPANSION-AN  =",E15.7/ 

2 lH  , "COEFFICIENT  OF  THERMAL  EXPANSION-AS  =",E15.7/ 

3 lH  , "COEFFICIENT  OF  THERMAL  EXPANSION-AT  =",E15.7/) 

2012  FORMAT  (lH  , "TEMPERATURE  =",E15.7/ 

llH  , "MODULUS  OF  ELASTICITY-EN  =",E15.7/ 

21H  , "MODULUS  OF  ELASTICITY-ES  =" ,E15.7/ 

3 1H  , "MODULUS  OF  ELASTICITY-ET  =",El5.7/ 

41H  , "POISSON  RATIO-NUNS  =",E15.7/ 

51H  , "POISSON  RATIO-NUNT  =",E15.7/ 

6 lH  , "POISSON  RATIO-NUST  =",E15.7/ 

7 1H  , "SHEAR  MODULUS-GNS  =",E15.7/ 

81H  , "SHEAR  MODULUS-GST  =",E15.7/ 

91H  , "SHEAR  MODULUS-GTN  =",El5.7/ 
llH  , "FREE  THERMAL  STRAIN-FN  =",E15.7/ 

21H  , "FREE  THERMAL  STRAIN-FS  =",E15.7/ 

31H  , "FREE  THERMAL  STRAIN-FT  =",El5.7/) 

2013  FORMAT  (30H1  PRESSURE  BOUNDARY  CONDITIONS/20H  I  J 

2014  FORMAT  (2I5,F10.1) 


ANGLE  A 


PRESSURE) 
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015  FORMAT  (27H1  SHEAR  BOUNDARY  CONDITIONS/17H  I  J  SHEAR) 

016  FORMAT  (26H  THE  SYSTEM  CONVERGED  IN  12  ,  llH  ITERATIONS) 

017  FORMAT  (33H  THE  SYSTEM  DID  NOT  CONVERGE  IN  I2,llH  ITERATIONS) 

024  FORMAT  (43H0  THE  AXISYMMETRIC  OPTION  NAS  BEEN  SELECTED) 

025  FORMAT (30H1  TORSION  BOUNDARY  CONDITIONS/17H  I  J  SHEAR) 

030  FORMAT ( lHl , " LOAD  STEP=",I4) 

031  FORMAT ( 1H0  , "AXIAL  ACCELERATION  =",El2.4/ 

11H0  , "ANGULAR  VELOCITY  =",E12.4/ 

2 1H0  , "ANGULAR  ACCELERATION^' ,E12 . 4) 

920  STOP 
END 

SUBROUTINE  ANGLE  (R, Z , RC, ZC , ANG) 

C  FIND  ANGLE  OF  INCLINATION  BETWEEN  0  AND  2*PI 

q*  *********************************** 

PI=3. 1415927 
Dl= (Z-ZC) 

D2= (R-RC) 

IF(ABS(R-RC) .GT.l.E-8)  GO  TO  100 
ANG=PI/2 . 

IF(D1.GT. l.E-8)  RETURN 

ANG=-ANG 

RETURN 

C*  *********************************** 

C  ALLOW  CIRCLE  TO  CROSS  AXIS 

C*  *********************************** 
100  ANG=ATAN2 (Dl ,D2) 

RETURN 

END 

SUBROUTINE  CIRCLE (ANGl , DELPHI , RSTRT , ZSTRT, RC ,  ZC  ,  I ,  J) 

INTEGER  CODE 

COMMON/TD/IMIN (100) ,IMAX(100) ,JMIN(25) ,JMAX(25) , MAX I , MAX J , NMTL , NBC 
COMMON/NPDATA/R( 1000)  ,CODE(1000)  ,XR(1000)  ,Z(1000)  ,XZ(1000)  , 

1NPNUM (25,80) ,T(1000) ,XT(1000) 

DIMENSION  AR(25, 80) ,AZ (25,80) 

EQUIVALENCE  ( R ( 1 )  , AR)  ,  ( Z  ( 1)  ,  AZ ) 

Q*  *********************************** 

C  FIND  INTERSECTION  OF  LINE  AND  CIRCLE  =  NEW  R  AND  Z 

C*  *********************************** 

ANG1=ANG1+DELPHI 

RR=SQRT ( (RSTRT-RC) * *2+ ( ZSTRT-ZC)  **2) 

AR ( I , J) =RC+RR*COS ( ANGl ) 

AZ (I, J) =ZC+RR*SIN (ANGl) 

RETURN 

END 


SUBROUTINE  MESH 
INTEGER  CODE 

DIMENSION  AR( 25,80) ,AZ(25,80) ,NCODE ( 25 , 80) 

COMMON/TD/IMIN (100) ,IMAX(100) ,JMIN(25) ,JMAX(25) , MAX I , MAX J , NMTL, NBC 
COMMON /NPDATA/R( 10  00)  , CODE (1000)  ,XR(1000)  ,Z(1000)  ,XZ(1000)  , 
lNPNUM (25,80) ,T(1000) ,XT(1000) 

COMMON/ELDATA/BETA ( 1000 ) ,EPR(1000) ,PR(200) ,SH(200) ,IX(1000,5) , 

IIP ( 200)  , JP ( 200 )  , IS  (  200)  , JS ( 200)  ,ALPHA(1000)  ,IT(200)  ,JT(200)  , 
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2ST (200) , INP (1000) ,JNP(1000) 

EQUIVALENCE  (R ( 1)  , AR)  ,  (Z(l)  ,AZ)  ,  (IX(1,1)  , NCODE) 

C*  *********************************** 
C  MESH  CONTROL  INFORMATION 

0*  *********************************** 
READ  (5,1000)  MAXI , MAX J, NS EG , NBC , NMTL 
WRITE(6, 2000)  MAXI , MAX J, NS  EG , NBC , NMTL 
0*  *********************************** 
C  INITIALIZE 

0*  *********************************** 
ISEG=-1 
PI=3 .1415927 
DO  110  J=1 ,80 
DO  100  1=1,25 
NCODE ( I , J) =0 
AR(I, J) =0. 

AZ ( I , J) =0  . 

JMAX (I) =0 
100  JMIN ( I ) =MAXI 
IMIN ( J) =MAXJ 
110  IMAX ( J) =0 

0*  *********************************** 
C  LINE  SEGMENT  CARDS 

0*  *********************************** 
150  ISEG=ISEG+1 

159  IF( ISEG. EQ.NSEG)  GO  TO  400 

READ (5,1001)  II , J1 , R1 , Z1 , 12 , J2 , R2 , Z2 , 13 , J3 , R3 , Z3 , IPTION 

WRITE (6, 2001) Il,Jl,Rl,Zl,I2, J2,R2,Z2, 13, J3,R3,Z3, IPTION 

IPTION=IPTION+l 

AR ( 1 1 , J 1 )  =Rl 

AZ (II, Jl) =Z1 

NCODE ( II, Jl) =1 

CALL  MNIMX ( II , Jl) 

GOTO  (150,200,200,300,300,200,200),  IPTION 
c*  *********************************** 
C  GENERATE  STRAIGHT  LINES  ON  BOUNDARY 

0*  *********************************** 
200  DI=  ABS ( FLOAT ( I2-I1) ) 

DJ=  ABS (FLOAT (J2-J1) ) 

AR(I2,J2) =R2 
AZ(I2,J2)=Z2 
NCODE ( 12, J2)  =1 
CALL  MNIMX (12 , J2 ) 

ISTRT=Il 

ISTP=I2 

JSTRT=Jl 

JSTP=J2 

DIFF=MAXl (DI,DJ) 

ITER=DIFF-1 . 

I INC=0 
JINC=0 

IF (I2.NE.il)  IINC=( 12-11) /IABS( 12-11) 
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IF( J2.NE. Jl)  JINC=( J2-J1) /IABS( J2-J1) 

KAPPA=1 

IF (12 .NE. II. AND. J2 .NE. Jl .AND. IPTION.NE. 3)  KAPPA=2 
IF (KAPPA. EQ. 2)  DIFF=2 . *DIFF 
RINC= (R2-R1) /DIFF 
ZINC= ( Z2-Z1) /DIFF 

WRITE (6, 2002)  DI , DJ, DIFF, RINC , ZINC , ITER, I INC , JINC , KAPPA 

CHECK  FOR  INPUT  ERROR 

IF (KAPPA. NE. 2 .OR. DI . EQ. DJ)  GO  TO  210 
WRITE(6, 2003) 

GO  TO  150 

INTERPOLATE 


210  1=11 
J=J1 

WRITE(6, 2004) 

DC  230  M=1 , ITER 

IF ( ITER. EQ. 0 .AND. IPTION. EQ. 2)  GO  TO  230 

IF ( ITER. EQ. 0 . AND. IPTION. EQ. 6)  GO  TO  230 

IF ( ITER. EQ. 0. AND. IPTION. EQ. 7)  GO  TO  230 

IF (KAPPA. EQ. 2)  GO  TO  220 

IOLD=I 

1=1+1 INC 

JOLD= J 

J=J+JINC 

AR( I , J) =AR ( IOLD , JOLD) +RINC 
AZ (I, J) =AZ (IOLD, JOLD) +ZINC 
WRITE(6, 2005)  I,J,AR(I,J) ,AZ(I,J) 

CALL  MN IMX ( I , J ) 

NCODE ( I , J) =1 
GO  TO  230 
220  CONTINUE 

IF ( II . GT . 12 . AND. IPTION. EQ. 7)  GO  TO  221 
IF ( II . LT . 12 . AND. IPTION. EQ . 6)  GO  TO  221 
IOLD=I 
1=1+1 INC 

AR( I , J) =AR ( IOLD , J) +RINC 
AZ (I, J) =AZ (IOLD, J) +ZINC 
WRITE (6 , 2005)  I,J,AR(I,J) ,AZ(I,J) 

NCODE ( I, J) =1 
CALL  MNIMX ( I , J) 

JOLD= J 
J=J+JINC 

AR ( I , J) =AR ( I , JOLD) +RINC 
AZ (I, J) =AZ (I, JOLD) +ZINC 
NCODE ( I, J) =1 

WRITE(6, 2005)  I,J,AR(I,J) ,AZ(I,J) 

CALL  MNIMX ( I, J) 

GO  TO  230 
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221  JOLD= J 
J=J+JINC 

AR ( I / J) =AR ( I , JOLD) +RINC 
AZ (I, J) =AZ (I, JOLD) +ZINC 
NCODE ( I / J) =1 

WRITE (6 , 2005)  I,J,AR(I,J) ,AZ(I,J) 

CALL  MNIMX(I,J) 

IOLD=I 
1=1+1 INC 

AR( I , J) =AR ( IOLD , J) +RINC 
AZ (I, J)  =AZ (IOLD, J) +ZINC 
NCODE ( I, J) =1 

WRITE(6, 2005)  I,J,AR(I,J) ,AZ(I,J) 

CALL  MNIMX ( I , J) 

230  CONTINUE 

IF (KAPPA. EQ. 1)  GO  TO  150 

IF(Il.GT. 12. AND. IPTION.EQ.7)  GO  TO  231 

IF(I1.LT.I2. AND . IPTION. EQ . 6)  GO  TO  231 

IOLD=I 

I=I+IINC 

AR ( I , J) =AR ( IOLD , J) +RINC 
AZ ( I, J) =AZ (IOLD, J) +ZINC 
GO  TO  232 

231  CONTINUE 
JOLD= J 
J=J+JINC 

AR ( I , J) =AR ( I , JOLD) +RINC 
AZ (I, J) =AZ (I, JOLD) +ZINC 

232  CONTINUE 
NCODE ( I, J) =1 

WRITE(6, 2005)  I,J,AR(I,J) ,AZ(I,J) 

CALL  MNIMX ( I , J) 

GO  TO  150 

Q*  *********************************** 

C  GENERATE  CIRCULAR  ARCS  ON  BOUNDARY 

C*  *********************************** 
300  AR( 12 , J2) =R2 
AZ (12, J2) =Z2 
NCODE ( 12, J2)  =  1 
CALL  MNIMX ( 12, J2) 

IF( IPTION. EQ. 5)  GO  TO  320 

FIND  CENTER  OF  CIRCLE 

AR ( 13 , J3 ) =R3 
AZ (13, J3) =Z3 
NCODE ( 13, J3) =1 
CALL  MNIMX ( 13, J3) 

SLAC= (Z2-Z1) /(R2-R1) 

SLBF=-1 . /SLAC 
SLCE= ( Z  3  —  Z  2 ) /  ( R3-R2 ) 

SLD  AR ( I , J) =AR ( IOLD , J) +RINC 

AZ (I, J)  =AZ (IOLD, J) +ZINC 
GO  TO  232 
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231  CONTINUE 
JOLD= J 
J= J+JINC 

AR(I, J) =AR( I , JOLD) +RINC 
AZ (I, J) =AZ (I, JOLD) +ZINC 

232  CONTINUE 
NCODE ( I ,  J) =1 

WRITE(6, 2005)  I ,  J, AR ( I ,  J) ,AZ(I,J) 

CALL  MNIMX ( I, J) 

GO  TO  150 

c*  *********************************** 
C  GENERATE  CIRCULAR  ARCS  ON  BOUNDARY 

Q*  *********************************** 

300  AR(I2, J2) =R2 
AZ (12, J2) =Z2 
NCODE ( 12, J2)  =  1 
CALL  MNIMX ( 1 2 , J2 ) 

IF ( IPTION. EQ . 5)  GO  TO  320 

FIND  CENTER  OF  CIRCLE 

AR (13 , J  3 ) =R3 
AZ (13, J3) =Z3 
NCODE ( 13 , J3 ) =1 
CALL  MNIMX ( 13, J3) 

SLAC= (Z2-Z1) / (R2-R1) 

SLBF=-1 . /SLAC 
SLCE= ( Z3-Z2) / (R3-R2) 

3LDF=-1 . /SLCE 

CHECK  FOR  INPUT  ERROR 

IF(ABS(SLAC-SLCE) .GT. .001)  GO  TO  310 
WRITE (6, 2006)  R1,Z1,R2,Z2,R3,Z3, SLAC, SLCE 
GO  TO  150 

310  R4=Rl+(R2-Rl)/2. 

Z4=Zl+(Z2-Zl)/2. 

R5=R2+ (R3-R2) /2 . 

Z5=Z2+ ( Z3-Z2 ) /2 . 

BBF=Z4-SLBF*R4 
BDF=Z5-SLDF*R5 
RC= (BBF-BDF) / ( SLDF-SLBF) 

ZC=SLBF*RC+BBF 
WRITE(6, 2007)  RC,ZC 
KAPPA=1 
GO  TO  330 
320  KAPPA=2 
RC=R3 
ZC=Z3 

330  ISTRT=I1 
ISTP=I2 
JSTRT= J1 


55 


UUU  U  UU  CJCJU 


JSTP=J2 

RSTRT=Rl 

RSTP=R2 

ZSTRT=Zl 

ZSTP=Z2 

340  CALL  ANGLE (RSTRT, ZSTRT, RC, ZC,ANG1) 

CALL  ANGLE ( RSTP , ZSTP ,  RC , ZC , ANG2) 

IF ( ANG2 . LE. ANG1)  ANG2=2 . 0*PI+ANG2 

FIND  ANGULAR  INCREMENT 

DI=  ABS ( FLOAT ( ISTP-ISTRT) ) 

DJ=  ABS (FLOAT (JSTP-JSTRT) ) 

I INC=0 
JINC=0 

IF ( ISTRT . NE. ISTP)  IINC= ( ISTP-ISTRT) /IABS (ISTP-ISTRT) 
I F  ( JSTRT .  NE .  vJSTP)  JINC=  ( JSTP-JSTRT)  /IABS  ( JSTP-JSTRT) 
LAMDA=1 

IF ( I INC . NE . O.AND. JINC.NE. 0)  LAMDA=2 
DIFF=MAX1 (DI,DJ) 

ITER=DIFF-1 . 

I F ( LAMDA . EQ . 2 )  DIFF=2.*DIFF 

DELPHI= (ANG2-ANG1) /DIFF 

WRITE ( 6 , 2008)  ANG1 ,ANG2 , DIFF, DELPHI 

CHECK  FOR  INPUT  ERROR 

IF ( LAMDA. NE . 2 . OR. DI . EQ . DJ)  GO  TO  350 
WRITE (6 , 2003) 

GO  TO  150 
350  IO= ISTRT 
JO=JSTRT 
WRITE(6, 2004) 

INTERPOLATE 

NPT=IABS( 12-11) +  IABS ( J2-J1)  -1 
DO  380  M=1 , ITER 

359  IF ( LAMDA. EQ. 2)  GO  TO  360 
1=10+1 INC 

J= JO+JINC 
CALL  MNIMX ( I , J) 

NCODE ( I , J) =1 

CALL  Cl RCLE(ANG1, DELPHI, RSTRT, ZSTRT, RC,ZC, I, J) 
WRITE(6, 2005)  I,J,AR(I,J) ,AZ(I,J) 

GO  TO  370 

360  1=10+1 INC 
J=JO 

NCODE ( I, J) =1 
CALL  MNIMX ( I, J) 

CALL  Cl RCLE(ANG1, DELPHI, RSTRT, ZSTRT, RC,ZC, I, J) 

WRITE (6,2005)  I,J,AR(I,J) ,AZ(I,J) 
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J=JO+JINC 
NCODE ( I / J)  =1 
CALL  MNIMX(I,J) 

CALL  CIRCLE ( ANG1 , DELPH I , RS TRT , ZSTRT , RC , ZC , I , J) 
WRITE(6, 2005)  I,J,AR(I,J) ,AZ(I,J) 

370  IO=I 
380  JO=J 

IF (LAMDA. NE. 2)  GO  TO  390 
1=10+1 INC 
NCODE ( I, J) =1 
CALL  MNIMX(I,J) 

CALL  CIRCLE ( ANG1 , DELPH I, RSTRT, ZSTRT, RC, ZC, I, J) 
WRITE (6,2005)  I,J,AR(I,J) ,AZ(I,J) 

390  IF (KAPPA. EQ. 2)  GO  TO  150 
ISTRT=I2 
ISTP=I3 
JSTRT= J2 
JSTP=J3 
RSTRT=R2 
RSTP=R3 
ZSTRT=Z2 
ZSTP=Z3 
KAPPA=2 

399  GO  TO  340 

C*  ************************* 

C  CALCULATE  COORDINATES  OF  INTERIOR  POINTS 

C*  ************************* 

400  IF (MAXJ. LE. 2)  GO  TO  430 
J2=MAXJ-1 

DO  420  N=1 ,500 
RESID=0 . 

DO  410  J=2,J2 
Il=IMIN(J) +1 
I2=IMAX ( J) -1 
DO  410  1=11,12 

IF(NCODE(I, J) .EQ. 1)  GO  TO  410 

DR= ( AR ( 1+1 , J) +AR( 1-1 , J) +AR ( I , J+l)  +AR(I, J-l)  )  /4 . 
DZ= ( AZ ( 1+1 , J) +AZ (1-1 ,J) +AZ (I,J+1) +AZ ( I , J-l) ) /A. 
RESID=RESID+ABS (DR)  +ABS(DZ) 

AR ( I , J) =AR ( I , J) +1 . 8*DR 
AZ(I, J) =AZ (I, J) +1.8*DZ 
410  CONTINUE 

IF(N.EQ.l)  RES1=RESID 
IF (N. EQ . 1 . AND. RE SID. EQ . 0 . ) GO  TO  430 
IF ( RESID/RES1 . LT. 1 . E-5)  GO  TO  430 
420  CONTINUE 
430  WRITE (6,2009)  N 

C*  ************************* 

CALL  POINTS 

C*  ************************* 

1C00  FORMAT  (515) 

1C01  FORMAT  (3(213, 2F8. 3) ,15) 


********** 

********** 


-AR  ( I ,  J) 
-AZ (I, J) 


********** 

********** 
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2000  FORMAT  (30H1  MESH  GENERATION  INFORMATION// 

1  41H0  MAXIMUM  VALUE  OF  I  IN  THE  MESH - 13/ 

2  41H0  MAXIMUM  VALUE  OF  J  IN  THE  MESH - 13/ 

3  41H0  NUMBER  OF  LINE  SEGMENT  CARDS - 13/ 

4  41H0  NUMBER  OF  BOUNDARY  CONDITION  CARDS - 13/ 

5  41H0  NUMBER  OF  MATERIAL  BLOCK  CARDS - 13///) 

2001  FORMAT  (//88H  INPUT  11  J1  Rl  Z1  12  J2  R2  Z 

12  13  J3  R3  Z3  IPTION/8X, 3 (214  , 2F8 .4)  ,16) 

2002  FORMAT  (5H  DI=F4.0,5H  DJ=F4.0,7H  DIFF=F4.0,7H  RINC=F8.3,7H  ZI 

1NC=F8.3,7H  ITER=I3 , 7H  IINC=I3,7H  JINC=I3,8H  KAPPA=Il) 

2003  FORMAT ( IX ,  3 8H* *BAD  INPUT — THIS  LINE  IS  NOT  DIAGONAL) 

2004  FORMAT  (30H  I  J  AR  AZ) 

2005  FORMAT  (2I5,2F11.6) 

2006  FORMAT  (51H  **  BAD  INPUT  -  THESE  POINTS  DO  NOT  DEFINE  A  CIRCLE,/, 
13X,6F12.4,10X,2E20.8) 

2007  FORMAT (19H  CENTER  COORDINATE, (Fll . 6 , IX , Fll . 6 , IX) ) 

2008  FORMAT  (7H  ANGl=F9.6,7H  ANG2=F9.6,7H  DIFF=F3.0,9H  DELPHI=F9.6) 

2009  FORMAT  (//30H  COORDINATES  CALCULATED  AFTER  13 , 11H  ITERATIONS) 
RETURN 

END 

SUBROUTINE  MNIMX(I,J) 

COMMON/TD/IMIN (100) , IMAX (100) , JMIN (25) ,JMAX(25) , MAX I , MAX J , NMTL , NBC 

IF ( J. LT. JMIN ( I) )  JMIN ( I) =J 

IF(J.GT. JMAX (I) )  JMAX ( I ) = J 

IF ( I . LT. IMIN ( J) )  IMIN ( J) =1 

IF ( I . GT . IMAX ( J) )  IMAX ( J) =1 

RETURN 

END 

SUBROUTINE  ADINA (NFACE , FTHETA) 

INTEGER  CODE 

COMMON/NPDATA/R (1000)  ,CODE(1000)  ,XR(1000)  ,Z  (1000)  ,XZ  (1000)  , 

1NPNUM (25,80) ,T(1000) ,XT(1000) 

COMMON/ARG/RRR ( 5 ) , ZZZ ( 5) ,RR(4) ,ZZ(4) ,S(15,15) ,P(15) ,TT(6) , 

1H ( 6 , 15) ,CRZ (6,6) ,XI(10) , ANGLE (4) ,SIG(18) ,EPS(18) ,N 
COMMON/E  LDATA/BETA (1000)  ,EPR(1000)  ,PR(200)  ,SH(200)  ,IX(1000,5)  , 

IIP ( 200) , JP ( 200) , IS ( 200) , JS ( 200) ,ALPHA(1000) ,IT(200) ,JT(200) , 

2ST (200) , INP (1000) ,JNP(1000) 

COMMON/BAS IC/ACELZ , ANGVEL , ANGACC , TRE  F , VOL , NUMNP , NUMEL , NUMPC , NUMSC , 
1NUMST 

COMMON/NXMESH/THETAN(4) , NST ( 4) ,NUMS(4,5) ,NPC(8,8) 

DIMENSION  NPG (25,80,9) , FTHETA ( 9) ,NCOD(8) ,NEL(8) ,ID(6) 

COMMON/ANSI /NUMELS (4) ,NUMNPS(4) 

COMMON/NXDATA/NTP , NTYPS, NTS , NTOTS 
COMMON/NONAX I /S 1 (30, 30) ,Pl(30) , THETA, BS1 (6 , 30) 

COMMON/SOLVE/X (4428) ,Y(4428) ,TEM(4428) , NUMTC , MB AND 

COMMON/TD/IMIN (100) ,IMAX(100) ,JMIN(25) ,JMAX(25) , MAX I , MAX J, NMTL, NBC 

COMMON/CONVRG/IDONE 

COMMON/P LANE/NPP 

COMMON/RESULT/BS (6,15) ,D(6,6) ,C(6,6) ,AR,BB(6,9) ,CNS(6,6) 
COMMON/MATP/RO ( 6 ) ,E ( 12 , 16 , 6) ,EE(16) ,AOFTS(6) 

COMMON/TOTAL/NELN 
DATA  NPG/18000*0/ 
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NI=NST  (NTP) 

C  LOOP  ON  #  OF  SEGMENTS  OF  TYPE  NTP 
DO  500  1=1, NI 
DO  200  J=1 , NUMNP 
NFL=NUMS (NTP, I) 

NFR=NFL+1 
IP1=INP ( J) 

JP1=JNP (J) 

NF=NFL 
DO  200  K=1 , 2 
NC=CODE ( J) +1 

IF (NPG ( IP1 , JP1 ,NF) .NE.0)GO  TO  200 
NPG(IP1,JP1,NF) =IP1+ ( JP1-1) *25+NF*200 
DO  10  11=1,8 
10  ID (II) =0 

IF(NC.EQ. 3.0R.NC.EQ. 4.0R.NC.EQ. 7 . OR. NC . EQ . 8 )  ID(1) =1 
IF(NF.EQ.l)  ID (3) =1 
IF (NF . EQ . NTOTS+1 )  ID(2)=1 

IF  (R  ( J)  .  EQ  .0.0)  ID  ( 2)  =1 
IF ( R ( J) .EQ.0.0)  ID ( 3 ) =1 

WRITE  NODAL  DATA  ON  FILE  2 


WRITE (2,100) NPG ( IP1 , JP1 ,NF) , (ID(KK) ,KK=1,6) ,Z(J) 
100  FORMAT ( "X" , 14 , IX, 14 , 515 , 3F10 . 4 , 4X , "0 " ) 

200  NF=NFR 


DO  350  J=1 ,NUMEL 
NF=NUMS (NTP, I) 

I1=IX(J,1) 

I2=IX ( J, 2) 

I3=IX ( J, 3) 

I4  =  IX  ( J, 4) 

IP1=INP ( II ) 

JPl=JNP ( II ) 

NEL ( 1 ) =NPG ( IP1 , JP1 , NFL) 

IF ( IP1 . EQ. 1) NEL ( 1) =NPG ( IP1 , JPl , 1) 
NEL ( 4 ) =NPG ( IP1 , JPl , NFR) 

IP1=INP (12) 

JP1=JNP ( 12 ) 

NEL (2) =NPG ( IP1 , JPl , NFL) 

NEL ( 3 ) =NPG ( IP1 , JPl , NFR) 

IP1=INP (14) 

JP1=JNP (14) 

NEL ( 5 ) =NPG ( IP1 , JPl , NFL) 

IF ( IP1 . EQ . 1) NEL (5) =NPG ( IP1 , JPl , 1) 
NEL ( 8 ) =NPG ( IP1 , JPl , NFR) 

IP1=INP (13 ) 

JP1=JNP ( 13 ) 

NEL ( 6 ) =NPG ( IP1 , JPl , NFL) 

NEL ( 7 ) =NPG ( IP1 , JPl , NFR) 

NELN=NELN+1 

NZ=0 


,R( J)  , FTHETA (NF) 
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N8=8 

WRITE (3,300) NELN,N8 , N8 ,NZ , IX ( J, 5) ,NELN, NZ , NZ 

300  FORMAT (8 15 , F10 . 0) 

WRITE (3,301)  (NEL (K) , K=1 , 8 ) 

301  FORMAT ( 8 15 ) 

WRITE  MATERIAL  AXES  DATA  ON  FILE  4 

NPI=NEL  ( .1)  +7000 
NPJ=NEL ( 1) +8000 
NPK=NEL ( 1) +9000 
WRITE (4,302)  NELN, NPI , NPJ , NPK 

302  FORMAT (4 15) 

WRITE  NODAL  POINTS  FOR  MATERIAL  AXES 

FTRAD=FTHETA (NF)  /57. 29578 
CA=COS (ALPHA (J) ) 

CB=COS ( BETA ( J) ) 

SA=SIN (ALPHA (J) ) 

SB=SIN (BETA ( J) ) 

CF=COS (FTRAD) 

SF=SIN (FTRAD) 

XPI  =  0 .0 
YPI=0 . 0 
ZPI=0.0 
XP J=CA*CB 

YPJ=CA*SB*CF+SA*SF 

ZPJ=CA*SB*SF-SA*CF 

XPK=SA*CB 

YPK=SA*SB*CF-CA*SF 

ZPK=SA*SB*SF+CA*CF 

WRITE (7 , 303) NPI ,NZ,NZ,NZ,NZ, NZ , NZ , XPI , YPI ,  ZPI 
WRITE (7, 303)NPJ,NZ,NZ,NZ,NZ,NZ,NZ,XPJ, YPJ, ZPJ 
WRITE (7, 303) NPK , NZ , NZ , NZ , NZ , NZ , N Z , XPK , YPK , ZPK 

303  FORMAT (IX, 14, IX, 14, 515, 3F10 . 4 , 4X , "0 " ) 

350  CONTINUE 

DO  500  L=1 , NUMPC 
IPl=INP ( IP ( L) ) 

JPl=JNP ( IP ( L) ) 

IP2=INP ( JP (L) ) 

JP2=JNP(JP(L) ) 

NEL  ( 1)  =NPG  ( IPl ,  JPl ,  NFL) 

NEL ( 2) =NPG ( IP2 , JP2 , NFL) 

NEL  (  3)  =NPG  ( IP2  ,  JP2  , NFR) 

NEL(4) =NPG ( IPl , JPl , NFR) 

NCUR=1 

WRITE (1,400) NCUR, (NEL ( K) ,K=1,4) 

WRITE (1,401)  PR (L) ,PR(L) ,PR(L) ,PR(L) 

400  FORMAT (915) 

401  FORMAT (5F10. 0,1 5) 

500  CONTINUE 
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END 

FUNCTION  NODE ( I , J) 

COMMON/TD/IMIN (100) ,IMAX(100) ,JMIN(25) ,JMAX(25) ,MAXI , MAXJ, NMTL, NBC 
NODE=0 

DO  100  JJ=1,J 
NSTART=IMIN ( JJ) 

NSTOP=IMAX (JJ) 

DO  100  I I=NSTART , NSTOP 
N  ODE  =NODE + 1 

IF(JJ.EQ. J. AND. II. EQ. I)  RETURN 
100  CONTINUE 
RETURN 
END 

SUBROUTINE  POINTS 
INTEGER  CODE 

COMMON/BAS IC/ACELZ , ANGVEL, ANGACC , TREF, VOL , NUMNP , NUMEL, NUMPC , NUMSC , 
1NUMST 

COMMON/MATP/RO ( 6 ) ,  E  ( 12 , 16 , 6)  ,  EE  ( 16) ,AOFTS(6) 

COMMON/NPDATA/R( 1000)  ,CODE(1000)  ,XR(1000)  ,Z(1000)  ,XZ(1000)  , 

1NPNUM (25,80) , T ( 1000) ,XT(1000) 

COMMON/ELDATA/BETA ( 1000 )  ,EPR(1000)  ,PR(200)  ,SH(200)  ,IX(1000,5)  , 

IIP ( 200)  , JP (200)  , IS  (  200)  , JS ( 200)  ,ALPHA(1000)  ,IT(200)  ,JT(200)  , 
2ST(200) ,INP(1000) ,JNP(1000) 

COMMON/SOLVE/X ( 4428) ,Y(4428) ,TEM(4428) , NUMTC , MB AND 

COMMON/TD/IMIN (100) ,IMAX(100) ,JMIN(25) ,JMAX(25) , MAX I , MAXJ, NMTL, NBC 
C  OMMON/P L ANE/NP  P 

lDIMENSI0N  AR (25,80)  , AZ ( 25 , 8 0)  ,MATRIL ( 100 , 5 )  , BLKANG (100)  , BLKALF ( 1 

DIMENSION  IBNG(100) ,NBNG(100) 

EQUIVALENCE  (R( 1)  , AR)  ,  (Z  (1)  ,  AZ) 

ESTABLISH  NODAL  POINT  INFORMATION 

************************************ 
NEL=0 
NODSUM=0 
DO  100  J=1 ,MAXJ 
NSTART=IMIN ( J) 

NS  TOP  =  IMAX ( J ) 

DO  100  I =NS TART, NSTOP 
100  NODSUM=NODSUM+l 
NELSUM=0 
J  JMAX  =MAX J - 1 
DO  110  JJ=1 , JJMAX 

NSTOP=MIN0 (IMAX (JJ) , IMAX ( JJ+1) ) -1 
NSTART=MAX0 ( IMIN ( JJ) ,IMIN(JJ+1) ) 

DO  110  II=NSTART, NSTOP 
110  NELSUM=NELSUM+1 
NUMNP =NODSUM 
NUMEL=NELSUM 
DO  120  J=1 ,MAXJ 
NSTART=IMIN ( J) 

NSTOP=IMAX ( J)  > 

DO  120  I=NSTART, NSTOP 
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NPNUM ( I , J) =NODE ( I , J) 

NP=NPNUM(I, J) 

INP  (NP)  =1 
JNP (NP) =J 
R  (NP) =AR(I, J) 

120  Z (NP)  =AZ (I, J) 

************************************ 
READ  AND  ASSIGN  BOUNDARY  CONDITIONS 

************************************ 
INITIALIZE 

************************************ 
DO  130  I=1,NUMNP 
CODE (I) =0 

IF ( R ( I) .EQ.O. . AND. NPP. EQ. 0)  CODE(I)=l. 

XR  ( I)  =0. 

XZ(I)=0. 

XT  (I) =0.0 
130  T(I) =0. 

IF(NBC.EQ.O)  GO  TO  210 
DO  200  IBCON=l ,  NBC 

READ (5,1002)  II , 12 , J1 , J2 , ICN, RCON, ZCON, TCON 
DO  200  1=11,12 
DO  200  J=J1,J2 
NP=NPNUM( I, J) 

CODE (NP) =ICN 
XR (NP) =RCON 
XT ( NP) =TCON 
200  XZ (NP) =ZCON 
210  MPRINT=0 

DO  230  J=1,MAXJ 
NSTART=IMIN( J) 

NS  TOP = IMAX ( J ) 

DO  230  I=NSTART, NSTOP 
NP=NPNUM ( I , J) 

IF (MPRINT. NE. 0)  GO  TO  220 
WRITE (6 , 2000) 

MPRINT=59 

220  MPRINT=MPRINT-1 

230  WRITE (6,2001)  I , J, NP ,CODE (NP) , R (NP) , Z (NP) , XR (NP) , XZ (NP) , XT (NP) 

c*  *********************************** 

C  ASSIGN  MATERIALS  IN  BLOCKS 

C*  *********************************** 
DO  300  Ml=l ,NUMEL 
300  IX (Ml, 5) =0 

DO  310  IMTL=1 ,NMTL 

READ  (5,1000)  MTL, (MATRIL ( IMTL, IM) ,IM=2,5) , BLKANG ( IMTL) , BLKALF ( IMT 
1L)  ,  IBNG  ( IMTL)  ,  NBNG  ( IMTL) 

310  MATRIL ( IMTL, 1) =MTL 

C*  *********************************** 
C  ESTABLISH  ELEMENT  INFORMATION 

C*  *********************************** 
J JMAX=MAXJ-1 
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N=0 

MTL=1 

KTL=1 

DO  440  JJ=1 , JJMAX 

NSTOP=MINO (IMAX ( JJ) , IMAX (JJ+1) ) -1 
NSTART=MAX0 (IMIN(JJ) ,IMIN(JJ+1) ) 

DO  440  II=NSTART,NSTOP 
NEL=NEL+1 

DO  400  IMTL=1,NMTL 
IF ( I I . LT.MATRIL ( IMTL, 2) )  GO  TO  400 
IF ( I I . GE . MATRIL ( IMTL, 3) )  GO  TO  400 
IF (JJ. LT.MATRIL (IMTL, 4) )  GO  TO  400 
IF( JJ .MR (T) GOO  KT  MATM(L)  4CIE  KQ.KTL)  GOO 

G4  K=KT 

0 (II.Q.NSTAXI.NE.NSTOP)  GO  TO  440 
M=NEL+1 
IANG=ICNG 
NANG=NCNG 
GO  TO  421 
420  I=NPNUM (II, JJ) 

=K  K= ( , + 

IX (M , )  12 

I)  =K 

X,L  I) 

AHgFIBE (KT 

N  * INTERRUPTED* 

get,diab 
* TERMINATED* 

/diab,p=yes , s=41 
#PAGES=  53 


TA 
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N=0 

MTL=1 

KTL=1 

DO  440  J J=1 , J JMAX 

NSTOP=MINO ( IMAX ( JJ) , IMAX ( JJ+1) ) -1 
NSTART=MAX0 (IMIN(JJ) ,IMIN(JJ+1) ) 

DO  440  I I=NSTART, NSTOP 
NEL=NEL+1 

DO  400  IMTL=1 , NMTL 

IF ( II. LT.MATRIL ( IMTL, 2) )  GO  TO  400 
IF ( I I . GE . MATRIL ( IMTL, 3) )  GO  TO  400 
IF (JJ. LT.MATRIL (IMTL, 4) )  GO  TO  400 
IF( JJ.GE.MATRIL( IMTL, 5) )  GO  TO  400 
KAT=IMTL 

MAT=MATRIL( IMTL, 1) 

400  CONTINUE 

IF (KAT. EQ. KTL)  GO  TO  410 

KTL=KAT 

MTL=MAT 

GO  TO  420 

410  IF ( I I . EQ. NSTART)  GO  TO  420 

IF(JJ.NE. JJMAX.OR. II.NE.NSTOP)  GO  TO  440 
M=NEL+1 
IANG=ICNG 
NANG=NCNG 
GO  TO  421 
420  I=NPNUM(I I, JJ) 

J=I  +  1 

K=NPNUM( II+l , JJ+1) 

L=K-1 

M=NEL 

IX (M, 1) =1 

IX (M, 2) =J 

IX (M, 3) =K 

IX (M, 4) =L 

IX (M, 5) =MTL 

BETA (M) =BLKANG (KTL) 

ALPHA (M) =BLKALF (KTL) 

I ANG=ICNG 
NANG=NCNG 
ICNG=I3NG (KTL) 

NCNG=NBNG (KTL) 

421  NC=2 
430  N=N+1 

IF(M.LE.N)  GO  TO  440 
IX (N , 1) =IX (N-l , 1) +1 
IX (N, 2) =IX (N-l , 2) +1 
IX (N, 3) =IX (N-l , 3) +1 
IX (N, 4) =IX (N-l , 4) +1 
IX (N, 5) =IX (N-l , 5) 

BETA (N) =BETA (N-l) 

IF ( IANG . EQ. 1)  GO  TO  442 


64 


non 


ALPHA ( N) = ALPHA ( N-l ) 

GO  TO  443 

442  CONTINUE 

IF(NC.GT.NANG)  GO  TO  444 
ALPHA (N) =ALPHA (N-l) 

GO  TO  443 
444  NC=1 

ALPHA ( N) =-ALPHA ( N-l ) 

443  CONTINUE 
NC=NC+1 

IF(M.GT.N)  GO  TO  430 
440  CONTINUE 

IF(NUMNP.GT. 2000)  WRITE(6 , 2002) 

******************************** 
SET  NODAL  POINT  TEMPERATURE  TO  REFERENCE  TEMPERATURE 
******************************** 
I F ( NUMTC . NE . 0 )  RETURN 
DO  500  N=1 , NUMNP 
500  T (N) =TREF 
1000  FORMAT  (515 , 2F10 . 0 , 215 ) 

1002  FORMAT (415,110, 3F10 . 0) 

2000  FORMAT  (128H1  I  J  NP  TYPE  R-ORDINATE 

1TE  R  LOAD  OR  DISPLACEMENT  Z  LOAD  OR  DISPLACEMENT  T  LOAD 
2LACEMENT) 

2001  FORMAT  ( 215 , 16 , 112 , F13 . 6 , F14 . 6 , E26 . 7 , E2 4 . 7 , E24 . 7) 

2002  FORMAT  (35H  BAD  INPUT  -  TOO  MANY  NODAL  POINTS) 

RETURN 

END 

SUBROUTINE  TEMP(R,Z,T) 

COMMON/SOLVE/X (4428) ,Y(4428) ,TEM(4428) , NUMTC, MB AND 
DIMENSION  SMALL (20) , ISM (20) 

fcfcfcifr**************  &&*********** 

C  INITIALIZE 

c*  ******************************* 
J=1 

JMAX=16 

IF ( NUMTC. LT. UMAX)  JMAX=NUMTC 
DO  10  1=1 , JMAX 
SMALL (I) =0. 

10  ISM (I) =0 

C*  ******************************* 
C  FIND  THE  JMAX  CLOSEST  POINTS 

0*  ******************************* 
DO  50  1=1, NUMTC 
DSQ=(X(I)  -R)  *  *2+  ( Y  ( I )  -Z)  **2 
IF(DSQ.GT. .IE-4)  GO  TO  20 
T=TEM ( I) 

RETURN 

20  IF(I.EQ.l)  SMALL (1)=DSQ 
IF(I.EQ.l)  ISM (1) =1 
IF(I.EQ.l)  GO  TO  50 

IF ( SMALL (J) . LE.DSQ.AND. J.LT. JMAX)  SMALL ( J+l ) =DSQ 


*  *  *  * 

*  *  *  * 

-ORDINA 
OR  DISP 


*  *  *  * 

*  *  *  * 


*  *  *  * 

*  *  *  * 
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IF ( SMALL (J) .LE.DSQ. AND. J. LT. JMAX)  ISM(J+1)=I 
IF ( SMALL (J) .LE.DSQ)  GO  TO  40 
DO  30  K=1 , J 
JB=J-K  +1 

IF(JB.EQ.O)  GO  TO  40 
SMALL (JB+1) =SMALL ( JB) 

ISM ( JB+1) =ISM ( JB) 

SMALL (JB) =DSQ 
ISM( JB) =1 

IF(JB.EQ.l)  GO  TO  40 
IF (SMALL (JB-1) .LE.DSQ)  GO  TO  40 
30  CONTINUE 

40  IF(J.LT. JMAX)  J=J+1 
50  CONTINUE 

C*  ********************** 

C  FIND  THE  THIRD  TEMPERATURE  POINT  BY  AREA 

0*  ********************** 

JCHK=JMAX-2 
J=0 

I 1 = I SM ( 1 ) 

12= ISM ( 2) 

60  I3=ISM ( J+3 ) 

ARE A= . 5  0  * ( Y (II) *X ( 1 3 ) - Y (13)  *X ( 1 1 ) +Y (13)  *X(I2)-Y(I2) *X(I3)  + 

1  Y(I2) *X(I1) -Y (II) *X(I2) ) 

Dl= (X ( 12) -X ( 11) ) **2+ ( Y (12) -Y(I1) ) **2 
C  IF  Dl  IS  APPROXIMATELY  0.  IT  IS  ASSUMED  THAT  THERE  EXISTS  A 

C  DUPLICATION  OF  INPUT 

IF (Dl .GT. . IE-3)  GO  TO  70 
12  =  13 
J=J+1 
GO  TO  60 

70  IF (AREA* *2 . GT. . l*Dl * SMALL ( 1) )  GO  TO  80 
J=J+1 

IF ( J. LT. JCHK)  GO  TO  60 
WRITE(6, 2000)  Il,I2fl3,J 
T=TEM ( II ) 

RETURN 

q*  *********************************** 
C  FIND  TEMPERATURE  INTERCEPT 

0*  *********************************** 
80  DETA=Y (II) * (TEM( 13) -TEM( 12) ) +Y (12) * (TEM( 11) -TEM( 13) ) 

1  +Y(I3) * (TEM(I2) -TEM(Il) ) 

DETB-X(Il) * ( TEM ( 12 ) -TEM ( 13 ) )+X(I2) * ( TEM ( 13 ) -TEM ( 11 ) ) 

1  +X  ( 1 3 )  *  ( TEM  (II)  -TEM  (12)) 

DETC=TEM ( 11 ) *(X(I2) *Y(I3)-X(I3) *Y(I2) )+TEM(I2) * (X ( 13) *Y( II) -X (II) * 
1Y  ( 13 )  )+TEM(I3)  *(X(II)  *Y  ( 12 )  -X  ( 12 )  *Y  ( 1 1)  ) 

T= ( DETA*R+DETB  *  Z  +DETC ) / ( 2 . *AREA) 

2000  FORMAT  (28H  ERROR  IN  TEMPERATURE  INPUT, 5H  Il=l4,5H  12=14, 

15H  13=14, 4H  J=I4) 

RETURN 

END 

SUBROUTINE  TEM2 (NUMNP) 


************* 

TEST 

************* 
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100 

1000 


N 


INTEGER  CODE 

COMMON /NPDATA/R ( 1000) ,CODE(1000) ,XR(1000) 
1NPNUM (25,80) ,T(1000) ,XT(1000) 

READ (5, 1000)  TCONST 
DO  100  N=1 , NUMNP 
T (N) =TCONST 
FORMAT (FI 0.0) 

RETURN 

END 

END 


,Z  (1000)  ,XZ (1000) 


/ 
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APPENDIX  C 


Listing  of  File  ADI, 
the  Input  Cards  for  ADINA  Program 
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X  20i  i 

X  401  1 

X  202  1 

X  402  1 

X  203  1 

X  403  1 

X  204  1 

X  404  1 

X  205  1 

X  405  1 

X  206  1 

X  406  1 

X  207  0 

X  407  0 

X  208  0 

X  408  0 

X  226  0 

X  426  0 

X  227  0 

X  427  0 

X  228  0 

X  428  0 

X  229  0 

X  429  0 

X  230  0 

X  430  0 

X  231  0 

X  431  0 

X  232  0 

X  432  0 

X  233  0 

X  433  0 

X  251  0 

X  451  0 

X  252  0 

X  452  0 

X  253  0 

X  453  0 

X  254  0 

X  454  0 

X  255  0 

X  455  0 

X  256  0 

X  456  0 

X  257  0 

X  457  0 

X  258  0 

X  458  0 

X  276  0 

X  476  0 

X  277  0 

X  477  0 

X  278  0 

X  478  0 

X  279  0 

X  479  0 

X  280  0 

X  480  0 

X  281  0 


110  0 
110  0 
110  0 
110  0 
110  0 
110  0 
110  0 
110  0 
110  0 
110  0 
110  0 
110  0 
0  10  0 
0  10  0 
0  10  0 
0  10  0 
110  0 
110  0 
0  10  0 
0  10  0 
0  10  0 
0  10  0 
0  10  0 
0  10  0 
0  10  0 
0  10  0 
0  10  0 
0  10  0 
0  10  0 
0  10  0 
0  10  0 
0  10  0 
110  0 
110  0 
0  10  0 
0  10  0 
0  10  0 
0  10  0 
0  10  0 
0  10  0 
0  10  0 
0  10  0 
0  10  0 
0  10  0 
0  10  0 
0  10  0 
0  10  0 
0  10  0 
110  0 
110  0 
0  10  0 
0  10  0 
0  10  0 
0  10  0 
0  10  0 
0  10  0 
0  10  0 
0  10  0 
0  10  0 


0 

7 ♦ 9950 

0.0000 

0.0000 

0  ADNPD 

0 

7*9950 

0.0000 

18.8000 

0 

0 

7*9950 

.2363 

0.0000 

0 

0 

7*9950 

.2363 

18.3000 

0 

0 

7*9950 

.4725 

0.0000 

0 

0 

7*9950 

.4725 

18.8000 

0 

0 

7*9950 

.5050 

0,0000 

0 

0 

7.9950 

.5050 

18.3000 

0 

0 

7.9950 

.5375 

0.0000 

0 

0 

7.9950 

.5375 

18.8000 

0 

0 

8*0218 

1.0467 

0.0000 

0 

0 

8.0218 

1.0467 

18.8000 

0 

0 

8.0486 

1.5558 

0.0000 

0 

0 

8*0486 

1.5558 

18,8000 

0 

0 

8*0754 

2.0650 

0.0000 

0 

0 

8*0754 

2.0650 

18.8000 

0 

0 

9.3024 

0.0000 

0 . 0000 

0 

0 

9.3024 

0.0000 

18.8000 

0 

0 

9.3234 

.2362 

0.0000 

0 

0 

9.3234 

.2362 

18.3000 

0 

0 

9*3024 

.4725 

0.0000 

0 

0 

9.3024 

.4725 

13,8000 

0 

0 

9.3230 

.5091 

0.0000 

0 

0 

9.3230 

.5091 

18.8000 

0 

0 

9.3024 

.5375 

0.0000 

0 

0 

9.3024 

.5375 

18.8000 

0 

0 

9.2469 

1.0393 

0.0000 

0 

0 

9.2469 

1.0393 

18.8000 

0 

0 

9.0533 

1.5509 

0.0000 

0 

0 

9.0533 

1.5509 

18.3000 

0 

0 

8.8254 

2.0650 

0 . 0000 

0 

0 

8.8254 

2.0650 

18.3000 

0 

0 

10.6098 

0.0000 

0.0000 

0 

0 

10.6098 

0.0000 

18.8000 

0 

0 

10.6944 

.2363 

0.0000 

0 

0 

10.6944 

.2363 

13.3000 

0 

0 

10.6098 

.4725 

0.0000 

0 

0 

10.6098 

.4725 

18.8000 

0 

0 

10.6919 

.5214 

0 . 0000 

0 

0 

10.6919 

.5214 

18.8000 

0 

0 

10.6098 

.5375 

0 . 0000 

0 

0 

10.6098 

.5375 

18.8000 

0 

0 

10.6098 

1  .0222 

0.0000 

0 

0 

10.6098 

1.0222 

18.3000 

0 

0 

10.0926 

1.5436 

0.0000 

0 

0 

10.0926 

1 .5436 

18.8000 

0 

0 

9.5754 

2.0650 

0.0000 

0 

0 

9.5754 

2.0650 

18.3000 

0 

0 

12.2166 

0.0000 

0.0000 

0 

0 

12.2166 

0.0000 

13.8000 

0 

0 

12.2338 

.2362 

0.0000 

0 

0 

12.2338 

.2362 

18.8000 

0 

0 

12.2300 

.4725 

0.0000 

0 

0 

12.2300 

.4725 

18.8000 

0 

0 

12.2255 

.5666 

0.0000 

0 

0 

12.2255 

.5666 

18.8000 

0 

0 

12.2211 

.6607 

0.0000 

0 

0 

12.2211 

♦  6607 

18.8000 

0 

0 

12.2166 

.7548 

0.0000 

0 
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X  481  0 

X  301  0 

X  501  0 

X  302  0 

X  502  0 

X  303  0 

X  503  0 

X  304  0 

X  504  0 

X  305  0 

X  505  0 

X  306  0 

X  506  0 

X  326  0 

X  526  0 

X  327  0 

X  527  0 

X  328  0 

X  528  0 

X  351  0 

X  551  0 

X  352  0 

X  552  0 

X  353  0 

X  553  0 

X  376  0 

X  576  0 

X  377  0 

X  577  0 

X  378  0 

X  578  0 

X  601  0 

X  602  0 

X  603  0 

X  604  0 

X  605  0 

X  606  0 

X  607  0 

X  608  0 

X  626  0 

X  627  0 

X  628  0 

X  629  0 

X  630  0 

X  631  0 

X  632  0 

X  633  0 

X  651  0 

X  652  0 

X  653  0 

X  654  0 

X  655  0 

X  656  0 

X  657  0 

X  658  0 

X  676  0 

X  677  0 

X  678  0 

X  679  0 

X  680  0 


0  1  0 

1  1  0 

1  1  0 

0  10 

0  10 

0  10 

0  10 

0  10 

0  1  o 

0  10 

0  10 

0  1  o 

0  1  o 

1  1  0 

1  1  0 

0  10 

0  10 

0  10 

0  1  o 

110 
1  1  0 

0  10 

0  10 

0  1  o 

0  1  o 

1  1  0 

1  1  0 

0  10 

0  1  o 

0  10 

0  1  o 

0  10 

0  10 

Oio 

0  10 

0  10 

Oio 
0  1  o 

0  1  o 

0  1  o 

Oio 
0  1  o 

0  10 

0  10 

0  10 

0  10 

0  10 

0  10 

0  10 

0  10 

0  1  o 

0  10 

0  10 

0  10 

Oio 

0  10 

0  1  o 

0  10 

0  1  o 

0  10 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


12.2166 

.7548 

18.8000 

0 

13.8234 

0.0000 

0.0000 

0 

13.8234 

0.0000 

18.8000 

0 

13.7950 

.2363 

0.0000 

0 

13.7950 

.2363 

18.3000 

0 

13.8234 

.4725 

0.0000 

0 

13.8234 

.4725 

18.8000 

0 

13.8234 

.4775 

0.0000 

0 

13.8234 

.4775 

18.8000 

0 

13.8234 

.4825 

0.0000 

0 

13.8234 

.4825 

18.8000 

0 

13.8234 

.4875 

0 . 0000 

0 

13.8234 

.4875 

13.8000 

0 

15.3068 

0.0000 

0.0000 

0 

15.3068 

0.0000 

18.8000 

0 

15.2992 

.2362 

0.0000 

0 

15.2992 

.2362 

18.8000 

0 

15.3068 

.4725 

0.0000 

0 

15.3068 

.4725 

18.8000 

0 

16.7903 

0.0000 

0.0000 

0 

16.7903 

0.0000 

18.8000 

0 

16.7883 

.2362 

0.0000 

0 

16.7883 

.2362 

18.3000 

0 

16.7903 

.4725 

0.0000 

0 

16.7903 

.4725 

18.8000 

0 

18.2737 

0.0000 

0.0000 

0 

18.2737 

0.0000 

18.3000 

0 

18.2737 

.2363 

0.0000 

0 

18.2737 

.2363 

18.8000 

0 

18.2737 

.4725 

0.0000 

0 

18.2737 

♦  4725 

18.8000 

0 

7.9950 

0.0000 

37.6000 

0 

7.9950 

.2363 

37.6000 

0 

7.9950 

.4725 

37.6000 

0 

7.9950 

.5050 

37.6000 

0 

7.9950 

.5375 

37.6000 

0 

8.0218 

1.0467 

37,6000 

0 

8.0486 

1.5558 

37,6000 

0 

8.0754 

2.0650 

37.6000 

0 

9.3024 

0.0000 

37.6000 

0 

9.3234 

.2362 

37.6000 

0 

9.3024 

.4725 

37.6000 

0 

9.3230 

.5091 

37.6000 

0 

9.3024 

.5375 

37,6000 

0 

9.2469 

1.0393 

37,6000 

0 

9.0533 

1.5509 

37.6000 

0 

8.8254 

2.0650 

37.6000 

0 

10 . 6098 
10  *  6944 
10.6098 
10.6919 
10.6098 
10.6098 
10.0926 
9.5754 
12.2166 
12.2338 
12.2300 
12.2255 
12.2211 


0.0000 

.2363 

.4725 

.5214 

.5375 

1.0222 

1.5436 

2.0650 

0.0000 

.2362 

.4725 

.5666 

.6607 


37.6000 

37.6000 

37.6000 

37.6000 

37.6000 

37.6000 

37.6000 

37.6000 

37.6000 

37.6000 

37.6000 

37.6000 

37.6000 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


71 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


0  10  0 
0  10  0 
0  10  0 
0  10  0 
0  10  0 
0  10  0 
0  10  0 
0  10  0 
0  10  0 
0  10  0 
0  10  0 
0  10  0 
0  10  0 
0  10  0 
0  10  0 
0  10  0 
0  10  0 
0  10  0 
0  10  0 
0  10  0 
0  10  0 
0  10  0 
0  10  0 
0  10  0 
0  10  0 
0  10  0 
0  10  0 
0  10  0 
0  10  0 
0  10  0 
0  10  0 
0  10  0 
0  10  0 
0  10  0 
0  10  0 
0  10  0 
0  10  0 
0  10  0 
0  10  0 
0  10  0 
0  10  0 
0  10  0 
0  10  0 
0  10  0 
0  10  0 
0  10  0 
0  10  0 
0  10  0 
0  10  0 
0  10  0 
0  10  0 
0  10  0 
0  10  0 
0  10  0 
0  10  0 
0  10  0 
0  10  0 
0  10  0 
0  10  0 
0  10  0 


0 

12.2166 

.7548 

37.6000 

0 

0 

13.8234 

0.0000 

37,6000 

0 

0 

13.7950 

.2363 

37  *  gOO 0 

0 

0 

13.8234 

.4725 

37.6000 

0 

0 

13.8234 

.4775 

37.6000 

0 

0 

13.8234 

.4825 

37,6000 

0 

0 

13.8234 

.4875 

37.6000 

0 

0 

15.3068 

0,0000 

37.6000 

0 

0 

15.2992 

.2362 

37.6000 

0 

0 

15.3068 

.4725 

37,6000 

0 

0 

16.7903 

0,0000 

37.6000 

0 

0 

16.7883 

.2362 

37.6000 

0 

0 

16.7903 

.4725 

37,6000 

0 

0 

18.2737 

0.0000 

37.6000 

0 

0 

18.2737 

.2363 

37,6000 

0 

0 

18.2737 

.4725 

37.6000 

0 

0 

7.9950 

0.0000 

56,4000 

0 

0 

7.9950 

.2363 

56.4000 

0 

0 

7.9950 

.4725 

56.4000 

0 

0 

7.9950 

.5050 

56.4000 

0 

0 

7.9950 

.5375 

56.4000 

0 

0 

8.0218 

1.0467 

56.4000 

0 

0 

8.0486 

1.5558 

56.4000 

0 

0 

8.0754 

2.0650 

56,4000 

0 

0 

9.3024 

0.0000 

56,4000 

0 

0 

9.3234 

.2362 

56.4000 

0 

0 

9.3024 

.4725 

56.4000 

0 

0 

9.3230 

.5091 

56.4000 

0 

0 

9.3024 

.5375 

56.4000 

0 

0 

9.2469 

1.0393 

56.4000 

0 

0 

9.0533 

1.5509 

56.4000 

0 

0 

8.8254 

2.0650 

56.4000 

0 

0 

10.6098 

0.0000 

56.4000 

0 

0 

10.6944 

.2363 

56.4000 

0 

0 

10.6098 

.4725 

56.4000 

0 

0 

10.6919 

.5214 

56.4000 

0 

0 

10.6098 

.5375 

56.4000 

0 

0 

10.6098 

1.0222 

56,4000 

0 

0 

10.0926 

1.5436 

56.4000 

0 

0 

9.5754 

2.0650 

56.4000 

0 

0 

12.2166 

0.0000 

56.4000 

0 

0 

12.2338 

.2362 

56.4000 

0 

0 

12.2300 

.4725 

56.4000 

0 

0 

12.2255 

.5  666 

56.4000 

0 

0 

12.2211 

.6607 

56.4000 

0 

0 

12.2166 

.7548 

56.4000 

0 

0 

13.8234 

0.0000 

56.4000 

0 

0 

13.7950 

.2363 

56.4000 

0 

0 

13.8234 

.4725 

56.4000 

0 

0 

13.8234 

.4775 

56.4000 

0 

0 

13.8234 

.4825 

56.4000 

0 

0 

13.8234 

.4875 

56.4000 

0 

0 

15.3068 

0.0000 

56.4000 

0 

0 

15.2992 

.2362 

56.4000 

0 

0 

15.3068 

.4725 

56.4000 

0 

0 

16.7903 

0.0000 

56.4000 

0 

0 

16,7883 

.2362 

56.4000 

0 

0 

16.7903 

.4725 

56.4000 

0 

0 

18.2737 

0.0000 

56.4000 

0 

0 

18.2737 

.2363 

56.4000 

0 

72 


X  978  0  0 

X1001  1  1 

X1002  1  1 

X1003  1  1 

X1004  1  1 

XI 005  1  i 

X1006  1  1 

X1007  0  0 

X1008  0  0 

X1026  0  1 

X1027  0  0 

X1028  0  0 

X1029  0  0 

X1030  0  0 

X1031  0  0 

X1032  0  0 

X1033  0  0 

X1051  0  1 

X1052  0  0 

X1053  0  0 

X1054  0  0 

X1055  0  0 

X1056  1  1 

X1057  1  1 

X1058  0  0 

X1076  0  1 

X1077  0  0 

X1078  0  0 

X1079  0  0 

X1080  0  0 

X1081  1  1 

X  882  1  1 

X1082  0  0 

X1101  0  1 

XI 102  0  0 

XI 103  0  0 

X1104  0  0 

X1105  0  0 

XI 106  1  1 

X  907  1  1 

XI 107  0  0 

XI 126  0  1 

X1127  0  0 

XI 128  0  0 

X1151  0  1 

XI 152  0  0 

XI 153  0  0 

XI 176  0  1 

XI 177  0  0 

X1178  0  0 

7201  0  0 

8201  0  0 

9201  0  0 

7202  0  0 

8202  0  0 

9202  0  0 

7203  0  0 

8203  0  0 

9203  0  0 

7204  0  0 


10  0 
10  0 
10  0 
10  0 
10  0 
10  0 
10  0 
10  0 
10  0 
10  0 
10  0 
10  0 
10  0 
10  0 
10  0 
10  0 
10  0 
10  0 
10  0 
10  0 
10  0 
1  0  0 
10  0 
10  0 
10  0 
10  0 
10  0 
10  0 
10  0 
10  0 
10  0 
10  0 
10  0 
10  0 
10  0 
10  0 
10  0 
10  0 
10  0 
10  0 
10  0 
10  0 
10  0 
10  0 
10  0 
10  0 
10  0 
10  0 
10  0 
10  0 
0  0  0 
0  0  0 
0  0  0 
0  0  0 
0  0  0 
0  0  0 
0  0  0 
0  0  0 
0  0  0 
0  0  0 


0 

18.2737 

.4725 

56.4000 

0 

0 

7.9950 

0.0000 

59.9000 

0 

0 

7.9950 

.2363 

59.9000 

0 

0 

7.9950 

.4725 

59.9000 

0 

0 

7.9950 

.5050 

59.9000 

0 

0 

7.9950 

.5375 

59.9000 

0 

0 

8.0218 

1.0467 

59.9000 

0 

0 

8.0486 

1.5558 

59.9000 

0 

0 

8.0754 

2.0650 

59.9000 

0 

0 

9.3024 

0.0000 

59,9000 

0 

0 

9.3237 

.2363 

59.9000 

0 

0 

9.3024 

.4725 

59.9000 

0 

0 

9.3229 

.5091 

59,9000 

0 

0 

9.3024 

.5375 

59.9000 

0 

0 

9.2273 

1.0241 

59.9000 

0 

0 

9.0485 

1.5471 

59.9000 

0 

0 

8.8254 

2.0650 

59.9000 

0 

0 

10.6098 

0.0000 

59.9000 

0 

0 

10.6942 

.2362 

59.9000 

0 

0 

10.6098 

.4725 

59.9000 

0 

0 

10.6920 

.5214 

59.9000 

0 

0 

10.6098 

.5375 

59.9000 

0 

0 

10.5366 

.9650 

59 . 9000 

0 

0 

10.0926 

1.5436 

59.9000 

0 

0 

9.5754 

2.0650 

59.9000 

0 

0 

12.2166 

0.0000 

59.9000 

0 

0 

12.2341 

.2363 

59,9000 

0 

0 

12.2300 

.4725 

59.9000 

0 

0 

12.2255 

.5666 

59.9000 

0 

0 

12.2211 

.6607 

59.9000 

0 

0 

12.2166 

.7548 

59,9000 

0 

0 

11.9580 

1.5436 

56.4000 

0 

0 

11.9580 

1.5436 

59.9000 

0 

0 

13.8234 

0.0000 

59.9000 

0 

0 

13.7950 

.2363 

59.9000 

0 

0 

13.8234 

.4725 

59.9000 

0 

0 

13.8234 

.4775 

59.9000 

0 

0 

13.8234 

.4825 

59.9000 

0 

0 

13.8234 

.4875 

59.9000 

0 

0 

13.8234 

1.5436 

56.4000 

0 

0 

13.8234 

1.5436 

59.9000 

0 

0 

15.3068 

0.0000 

59.9000 

0 

0 

15.2993 

.2363 

59.9000 

0 

0 

15.3068 

.4725 

59.9000 

0 

0 

16.7903 

0.0000 

59.9000 

0 

0 

16.7884 

.2363 

59.9000 

0 

0 

16.7903 

.4725 

59.9000 

0 

0 

18.2737 

0.0000 

59.9000 

0 

0 

18.2737 

.2363 

59.9000 

0 

0 

18.2737 

.4725 

59.9000 

0 

0 

0.0000 

0.0000 

0.0000 

0  ADNPX 

0 

1.0000 

0.0000 

0.0000 

0 

0 

0.0000 

0.0000 

1 . 0000 

0 

0 

0,0000 

0.0000 

0.0000 

0 

0 

1.0000 

0.0000 

0.0000 

0 

0 

0.0000 

0.0000 

1.0000 

0 

0 

0.0000 

0.0000 

0 . 0000 

0 

0 

1.0000 

0.0000 

0.0000 

0 

0 

0.0000 

0.0000 

1.0000 

0 

0 

0.0000 

0.0000 

0.0000 

0 

73 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


0 

0 

0 

0 

0 

1.0000 

0.0000 

0 . 0000 

0 

0 

0 

0 

0 

0 

0.0000 

0.0000 

1.0000 

0 

0 

0 

0 

0 

0 

0.0000 

0.0000 

0.0000 

0 

0 

0 

0 

0 

0 

1.0000 

0.0000 

0.0000 

0 

0 

0 

0 

0 

0 

0.0000 

0.0000 

1.0000 

0 

0 

0 

0 

0 

0 

0.0000 

0.0000 

0.0000 

0 

0 

0 

0 

0 

0 

1.0000 

0.0000 

0,0000 

0 

0 

0 

0 

0 

0 

0.0000 

0.0000 

1.0000 

0 

0 

0 

0 

0 

0 

0.0000 

0.0000 

0,0000 

0 

0 

0 

0 

0 

0 

1.0000 

0.0000 

0.0000 

0 

0 

0 

0 

0 

0 

0.0000 

0.0000 

1.0000 

0 

0 

0 

0 

0 

0 

0.0000 

0.0000 

0.0000 

0 

0 

0 

0 

0 

0 

1.0000 

0.0000 

0.0000 

0 

0 

0 

0 

0 

0 

0.0000 

0.0000 

1.0000 

0 

0 

0 

0 

0 

0 

0.0000 

0.0000 

0.0000 

0 

0 

0 

0 

0 

0 

1.0000 

0.0000 

0.0000 

0 

0 

0 

0 

0 

0 

0.0000 

0.0000 

1.0000 

0 

0 

0 

0 

0 

0 

0.0000 

0.0000 

0.0000 

0 

0 

0 

0 

0 

0 

1.0000 

0.0000 

0.0000 

0 

0 

0 

0 

0 

0 

0.0000 

0.0000 

1.0000 

0 

0 

0 

0 

0 

0 

0 . 0000 

0.0000 

0,0000 

0 

0 

0 

0 

0 

0 

1.0000 

0.0000 

0.0000 

0 

0 

0 

0 

0 

0 

0.0000 

0.0000 

1  .0000 

0 

0 

0 

0 

0 

0 

0.0000 

0.0000 

0.0000 

0 

0 

0 

0 

0 

0 

1.0000 

0.0000 

0.0000 

0 

0 

0 

0 

0 

0 

0.0000 

0.0000 

1.0000 

0 

0 

0 

0 

0 

0 

0.0000 

0.0000 

0.0000 

0 

0 

0 

0 

0 

0 

1.0000 

0.0000 

0.0000 

0 

0 

0 

0 

0 

0 

0.0000 

0.0000 

1 .0000 

0 

0 

0 

0 

0 

0 

0.0000 

0.0000 

0.0000 

0 

0 

0 

0 

0 

0 

1.0000 

0.0000 

0.0000 

0 

0 

0 

0 

0 

0 

0.0000 

0.0000 

1.0000 

0 

0 

0 

0 

0 

0 

0.0000 

0.0000 

0.0000 

0 

0 

0 

0 

0 

0 

1 . 0000 

0,0000 

0.0000 

0 

0 

0 

0 

0 

0 

0.0000 

0.0000 

1.0000 

0 

0 

0 

0 

0 

0 

0.0000 

0.0000 

0.0000 

0 

0 

0 

0 

0 

0 

1.0000 

0.0000 

0.0000 

0 

0 

0 

0 

0 

0 

0.0000 

0.0000 

1.0000 

0 

0 

0 

0 

0 

0 

0.0000 

0.0000 

0.0000 

0 

0 

0 

0 

0 

0 

1.0000 

0.0000 

0.0000 

0 

0 

0 

0 

0 

0 

0.0000 

0.0000 

1.0000 

0 

0 

0 

0 

0 

0 

0.0000 

0.0000 

0.0000 

0 

0 

0 

0 

0 

0 

1.0000 

0.0000 

0.0000 

0 

0 

0 

0 

0 

0 

0.0000 

0.0000 

1.0000 

0 

0 

0 

0 

0 

0 

0.0000 

0.0000 

0.0000 

0 

0 

0 

0 

0 

0 

1.0000 

0.0000 

0.0000 

0 

0 

0 

0 

0 

0 

0.0000 

0.0000 

1 . 0000 

0 

0 

0 

0 

0 

0 

0.0000 

0.0000 

0.0000 

0 

0 

0 

0 

0 

0 

1.0000 

0.0000 

0.0000 

0 

0 

0 

0 

0 

0 

0.0000 

0.0000 

1.0000 

0 

0 

0 

0 

0 

0 

0 . 0000 

0.0000 

0.0000 

0 

0 

0 

0 

0 

0 

1.0000 

0.0000 

0.0000 

0 

0 

0 

0 

0 

0 

0.0000 

0.0000 

1.0000 

0 

0 

0 

0 

0 

0 

0.0000 

0.0000 

0.0000 

0 

0 

0 

0 

0 

0 

1.0000 

0.0000 

0.0000 

0 

0 

0 

0 

0 

0 

0.0000 

0.0000 

1.0000 

0 

0 

0 

0 

0 

0 

0.0000 

0.0000 

0.0000 

0 

0 

0 

0 

0 

0 

1.0000 

0.0000 

0.0000 

0 

0 

0 

0 

0 

0 

0.0000 

0.0000 

1 .0000 

0 

0 

0 

0 

0 

0 

0.0000 

0.0000 

0.0000 

0 

74 


8280 

9280 

7301 

8301 

9301 

7302 

8302 

9302 

7326 

8326 

9326 

7327 

8327 

9327 

7351 

8351 

9351 

7352 

8352 

9352 

7201 

8201 

9201 

7402 

8402 

9402 

7403 

8403 

9403 

7404 

8404 

9404 

7405 

8405 

9405 

7406 

8406 

9406 

7407 

8407 

9407 

7226 

8226 

9226 

7427 

8427 

9427 

7428 

8428 

9428 

7429 

8429 

9429 

7430 

8430 

9430 

7431 

8431 

9431 

7432 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


1.0000 

0,0000 

0.0000 

0.0000 

0,0000 

1.0000 

0.0000 

0.0000 

0.0000 

1.0000 

0,0000 

0.0000 

0.0000 

0.0000 

1.0000 

0.0000 

0.0000 

0.0000 

1,0000 

0.0000 

0.0000 

0.0000 

0.0000' 

1.0000 

0.0000 

0,0000 

0.0000 

1.0000 

0.0000 

0,0000 

0.0000 

0,0000 

1.0000 

0.0000 

0.0000 

0.0000 

1.0000 

0.0000 

0.0000 

0.0000 

0.0000 

1,0000 

0.0000 

0.0000 

0.0000 

1.0000 

0.0000 

0.0000 

0.0000 

0.0000 

1.0000 

0.0000 

0.0000 

0.0000 

1,0000 

0.0000 

0.0000 

0.0000 

0.0000 

1.0000 

0.0000 

0.0000 

0.0000 

1.0000 

0.0000 

0.0000 

0.0000 

-.3223 

.9466 

0.0000 

0.0000 

0.0000 

1.0000 

0.0000 

0.0000 

0.0000 

-.3223 

.9466 

0.0000 

0.0000 

0,0000 

1.0000 

0.0000 

0.0000 

0.0000 

-.3223 

.9466 

0.0000 

0.0000 

0.0000 

1.0000 

0.0000 

0.0000 

0.0000 

-.3223 

.9466 

0.0000 

0.0000 

0,0000 

1.0000 

0.0000 

0.0000 

0,0000 

-.3223 

.9466 

0.0000 

0.0000 

0.0000 

1.0000 

0.0000 

0,0000 

0.0000 

-.3223 

.9466 

0.0000 

0,0000 

0.0000 

1.0000 

0.0000 

0.0000 

0.0000 

-.3223 

.9466 

0.0000 

0.0000 

0.0000 

1.0000 

0.0000 

0.0000 

0.0000 

-.3223 

.9466 

0.0000 

0.0000 

0.0000 

1.0000 

0.0000 

0.0000 

0.0000 

-.3223 

.9466 

0.0000 

0.0000 

0.0000 

1.0000 

0.0000 

0.0000 

0.0000 

-.3223 

.9466 

0.0000 

0.0000 

0.0000 

1.0000 

0.0000 

0.0000 

0.0000 

-.3223 

.9466 

0.0000 

0.0000 

0.0000 

1.0000 

0.0000 

0.0000 

0.0000 

-.3223 

.9466 

0.0000 

0.0000 

0.0000 

1.0000 

0.0000 

0.0000 

0.0000 

-.3223 

.9466 

0.0000 

0.0000 

0.0000 

75 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


0 

0 

0 

0 

0 

1.0000 

0.0000 

0.0000 

0 

0 

0 

0 

0 

0 

0.0000 

-.3223 

.9466 

0 

0 

0 

0 

0 

0 

0.0000 

0.0000 

0.0000 

0 

0 

0 

0 

0 

0 

1.0000 

0.0000 

0.0000 

0 

0 

0 

0 

0 

0 

0.0000 

-.3223 

.9466 

0 

0 

0 

0 

0 

0 

0.0000 

0.0000 

0.0000 

0 

0 

0 

0 

0 

0 

1.0000 

0.0000 

0.0000 

0 

0 

0 

0 

0 

0 

0.0000 

-.3223 

.9466 

0 

0 

0 

0 

0 

0 

0,0000 

0.0000 

0.0000 

0 

0 

0 

0 

0 

0 

1.0000 

0.0000 

0.0000 

0 

0 

0 

0 

0 

0 

0.0000 

-.3223 

.9466 

0 

0 

0 

0 

0 

0 

0.0000 

0.0000 

0.0000 

0 

0 

0 

0 

0 

0 

1,0000 

0,0000 

0.0000 

0 

0 

0 

0 

0 

0 

0.0000 

-.3223 

.9466 

0 

0 

0 

0 

0 

0 

0.0000 

0.0000 

0.0000 

0 

0 

0 

0 

0 

0 

1.0000 

0.0000 

0.0000 

0 

0 

0 

0 

0 

0 

0.0000 

-.3223 

.9466 

0 

0 

0 

0 

0 

0 

0,0000 

0,0000 

0.0000 

0 

0 

0 

0 

0 

0 

1.0000 

0,0000 

0.0000 

0 

0 

0 

0 

0 

0 

0.0000 

-.3223 

.9466 

0 

0 

0 

0 

0 

0 

0.0000 

0.0000 

0.0000 

0 

0 

0 

0 

0 

0 

1.0000 

0.0000 

0.0000 

0 

0 

0 

0 

0 

0 

0,0000 

-.3223 

.9466 

0 

0 

0 

0 

0 

0 

0.0000 

0,0000 

0.0000 

0 

0 

0 

0 

0 

0 

1,0000 

0,0000 

0.0000 

0 

0 

0 

0 

0 

0 

0.0000 

-.3223 

.9466 

0 

0 

0 

0 

0 

0 

0.0000 

0.0000 

0.0000 

0 

0 

0 

0 

0 

0 

1.0000 

0.0000 

0.0000 

0 

0 

0 

0 

0 

0 

0.0000 

-.3223 

.9466 

0 

0 

0 

0 

0 

0 

0.0000 

0.0000 

0.0000 

0 

0 

0 

0 

0 

0 

1.0000 

0.0000 

0.0000 

0 

0 

0 

0 

0 

0 

0.0000 

-.3223 

» 9466 

0 

0 

0 

0 

0 

0 

0.0000 

0.0000 

0.0000 

0 

0 

0 

0 

0 

0 

1.0000 

0.0000 

0.0000 

0 

0 

0 

0 

0 

0 

0.0000 

-.3223 

.9466 

0 

0 

0 

0 

0 

0 

0.0000 

0.0000 

0,0000 

0 

0 

0 

0 

0 

0 

1.0000 

0.0000 

0.0000 

0 

0 

0 

0 

0 

0 

0.0000 

-.3223 

.9466 

0 

0 

0 

0 

0 

0 

0.0000 

0.0000 

0.0000 

0 

0 

0 

0 

0 

0 

1.0000 

0.0000 

0.0000 

0 

0 

0 

0 

0 

0 

0.0000 

-.3223 

♦  9466 

0 

0 

0 

0 

0 

0 

0.0000 

0.0000 

0,0000 

0 

0 

0 

0 

0 

0 

1.0000 

0.0000 

0.0000 

0 

0 

0 

0 

0 

0 

0.0000 

-.3223 

.9466 

0 

0 

0 

0 

0 

0 

0.0000 

0.0000 

0,0000 

0 

0 

0 

0 

0 

0 

1.0000 

0.0000 

0.0000 

0 

0 

0 

0 

0 

0 

0.0000 

-.3223 

.9466 

0 

0 

0 

0 

0 

0 

0.0000 

0.0000 

0.0000 

0 

0 

0 

0 

0 

0 

1.0000 

0.0000 

0.0000 

0 

0 

0 

0 

0 

0 

0.0000 

-.3223 

.9466 

0 

0 

0 

0 

0 

0 

0.0000 

0.0000 

0.0000 

0 

0 

0 

0 

0 

0 

1.0000 

0.0000 

0.0000 

0 

0 

0 

0 

0 

0 

0.0000 

-.6101 

.7923 

0 

0 

0 

0 

0 

0 

0.0000 

0.0000 

0.0000 

0 

0 

0 

0 

0 

0 

1.0000 

0.0000 

0.0000 

0 

0 

0 

0 

0 

0 

0.0000 

-.6101 

.7923 

0 

0 

0 

0 

0 

0 

0.0000 

0.0000 

0.0000 

0 

0 

0 

0 

0 

0 

1 .0000 

0.0000 

0.0000 

0 

0 

0 

0 

0 

0 

0.0000 

-.6101 

.7923 

0 

0 

0 

0 

0 

0 

0.0000 

0.0000 

0.0000 

0 

76 


8604 

9604 

7605 

8605 

9605 

7606 

8606 

9606 

7607 

8607 

9607 

7226 

8226 

9226 

7627 

8627 

9627 

7628 

8628 

9628 

7629 

8629 

9629 

7630 

8630 

9630 

7631 

8631 

9631 

7632 

8632 

9632 

7251 

8251 

9251 

7652 

8652 

9652 

7653 

8653 

9653 

7654 

8654 

9654 

7655 

8655 

9655 

7276 

8276 

9276 

7677 

8677 

9677 

7678 

8678 

9678 

7679 

8679 

9679 

7680 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


1.0000 

0.0000 

0 . 0000 

0.0000 

-.6101 

.7923 

0.0000 

0.0000 

0.0000 

1.0000 

0.0000 

0.0000 

0 . 0000 

-.6101 

.7923 

0.0000 

0.0000 

0.0000 

1.0000 

0.0000 

0.0000 

0.0000 

-.6101 

.7923 

0.0000 

0.0000 

0.0000 

1.0000 

0.0000 

0.0000 

0.0000 

-.6101 

.7923 

0.0000 

0.0000 

0 . 0000 

1.0000 

0.0000 

0.0000 

0.0000 

-.6101 

.7923 

0.0000 

0.0000 

0.0000 

1.0000 

0.0000 

0.0000 

0.0000 

-.6101 

.7923 

0.0000 

0.0000 

0.0000 

1.0000 

0.0000 

0.0000 

0.0000 

-.6101 

.7923 

0.0000 

0.0000 

0.0000 

1.0000 

0.0000 

0.0000 

0.0000 

-.6101 

.7923 

0.0000 

0.0000 

0.0000 

1.0000 

0.0000 

0.0000 

0.0000 

-.6101 

.7923 

0.0000 

0.0000 

0.0000 

1.0000 

0.0000 

0.0000 

0.0000 

-.6101 

.7923 

0.0000 

0.0000 

0.0000 

1.0000 

0.0000 

0.0000 

0.0000 

-.6101 

.7923 

0.0000 

0.0000 

0.0000 

1.0000 

0.0000 

0.0000 

0.0000 

-.6101 

.7923 

0.0000 

0.0000 

0.0000 

1.0000 

0.0000 

0.0000 

0.0000 

-.6101 

.7923 

0.0000 

0.0000 

0.0000 

1.0000 

0.0000 

0.0000 

0.0000 

-.6101 

.7923 

0 . 0000 

0.0000 

0.0000 

1.0000 

0.0000 

0.0000 

0.0000 

-.6101 

.7923 

0.0000 

0.0000 

0.0000 

1.0000 

0.0000 

0.0000 

0.0000 

-.6101 

.7923 

0.0000 

0.0000 

0 . 0000 

1.0000 

0.0000 

0.0000 

0.0000 

-.6101 

.7923 

0.0000 

0.0000 

0.0000 

1.0000 

0.0000 

0.0000 

0.0000 

-.6101 

.7923 

0.0000 

0.0000 

0.0000 

1.0000 

0.0000 

0.0000 

0.0000 

-.6101 

.7923 

0.0000 

0.0000 

0.0000 

1.0000 

0.0000 

0.0000 

0,0000 

-.6101 

.7923 

0.0000 

0.0000 

0.0000 

77 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


0  0  0 
0  0  0 
0  0  0 
0  0  0 
0  0  0 
0  0  0 
0  0  0 
0  0  0 
0  0  0 
0  0  0 
0  0  0 
0  0  0 
0  0  0 
0  0  0 
0  0  0 
0  0  0 
0  0  0 
0  0  0 
0  0  0 
0  0  0 
0  0  0 
0  0  0 
0  0  0 
0  0  0 
0  0  0 
0  0  0 
0  0  0 
0  0  0 
0  0  0 
0  0  0 
0  0  0 
0  0  0 
0  0  0 
0  0  0 
0  0  0 
0  0  0 
0  0  0 
0  0  0 
0  0  0 
0  0  0 
0  0  0 
0  0  0 
0  0  0 
0  0  0 
0  0  0 
0  0  0 
0  0  0 
0  0  0 
0  0  0 
0  0  0 
0  0  0 
0  0  0 
0  0  0 
0  0  0 
0  0  0 
0  0  0 
0  0  0 
0  0  0 
0  0  0 
0  0  0 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


0 

1.0000 

0.0000 

0.0000 

0 

0 

0 . 0000 

-.6101 

.7923 

0 

0 

0.0000 

0.0000 

0.0000 

0 

0 

1.0000 

0.0000 

0.0000 

0 

0 

0.0000 

-.6101 

.7923 

0 

0 

0.0000 

0.0000 

0.0000 

0 

0 

1.0000 

0.0000 

0.0000 

0 

0 

0.0000 

-.6101 

.7923 

0 

0 

0.0000 

0.0000 

0.0000 

0 

0 

1.0000 

0.0000 

0.0000 

0 

0 

0.0000 

-.6101 

.7923 

0 

0 

0.0000 

0.0000 

0.0000 

0 

0 

1.0000 

0.0000 

0.0000 

0 

0 

0.0000 

-.6101 

.7923 

0 

0 

0.0000 

0.0000 

0.0000 

0 

0 

1.0000 

0.0000 

0.0000 

0 

0 

0.0000 

-.6101 

.7923 

0 

0 

0.0000 

0.0000 

0.0000 

0 

0 

1.0000 

0.0000 

0.0000 

0 

0 

0.0000 

-.6101 

.7923 

0 

0 

0.0000 

0.0000 

0 . 0000 

0 

0 

1.0000 

0.0000 

0.0000 

0 

0 

0.0000 

-.8329 

.5534 

0 

0 

0.0000 

0.0000 

0.0000 

0 

0 

1.0000 

0.0000 

0.0000 

0 

0 

0.0000 

-.8329 

.5534 

0 

0 

0.0000 

0.0000 

0.0000 

0 

0 

1.0000 

0.0000 

0.0000 

0 

0 

0.0000 

-.8329 

.5534 

0 

0 

0.0000 

0.0000 

0.0000 

0 

0 

1 . 0000 

0.0000 

0 . 0000 

0 

0 

0.0000 

-.8329 

.5534 

0 

0 

0.0000 

0.0000 

0.0000 

0 

0 

1.0000 

0.0000 

0.0000 

0 

0 

0.0000 

-.8329 

.5534 

0 

0 

0.0000 

0.0000 

0.0000 

0 

0 

1.0000 

0.0000 

0.0000 

0 

0 

0 » 0000 

-.8329 

.5534 

0 

0 

0.0000 

0.0000 

0.0000 

0 

0 

1 . 0000 

0.0000 

0.0000 

0 

0 

0 . 0000 

-.8329 

.5534 

0 

0 

0.0000 

0.0000 

0.0000 

0 

0 

1 . 0000 

0.0000 

0.0000 

0 

0 

0.0000 

-.8329 

.5534 

0 

0 

0.0000 

0.0000 

0.0000 

0 

0 

1.0000 

0.0000 

0.0000 

0 

0 

0 . 0000 

-.8329 

.5534 

0 

0 

0.0000 

0.0000 

0 , OOOO 

0 

0 

1.0000 

0.0000 

0.0000 

0 

0 

0.0000 

-.8329 

.5534 

0 

0 

0.0000 

0.0000 

0.0000 

0 

0 

1.0000 

0.0000 

0.0000 

0 

0 

0.0000 

-.8329 

.5534 

0 

0 

0.0000 

0.0000 

0.0000 

0 

0 

1.0000 

0.0000 

0.0000 

0 

0 

0.0000 

-.8329 

.5534 

0 

0 

0.0000 

0.0000 

0.0000 

0 

0 

1.0000 

0.0000 

0.0000 

0 

0 

0.0000 

-.8329 

.5534 

0 

0 

0.0000 

0.0000 

0.0000 

0 

78 


8832 

9832 

7251 

8251 

9251 

7852 

8852 

9852 

7853 

8853 

9853 

7854 

8854 

9854 

7855 

8855 

9855 

7856 

8856 

9856 
7276 
8276 
9276 

7877 

8877 

9877 

7878 

8878 

9878 

7879 

8879 

9879 

7880 

8880 

9880 

7881 
8881 

9881 
7301 
8301 
9301 
7902 
8902 
9902 
7326 
8326 
9326 
7927 
8927 
9927 
7351 
8351 
9351 
7952 
8952 
9952 

3  122 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


0 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


0 


10000000 


10000 

10000000* 


10000000. 


0 

0 

1.0000 

0.0000 

0.0000 

0 

0 

0 

0.0000 

-.8329 

.5534 

0 

0 

0 

0.0000 

0.0000 

0.0000 

0 

0 

0 

1.0000 

0.0000 

0.0000 

0 

0 

0 

0.0000 

-.8329 

.5534 

0 

0 

0 

0.0000 

0.0000 

0,0000 

0 

0 

0 

1.0000 

0.0000 

0.0000 

0 

0 

0 

0.0000 

-.8329 

.5534 

0 

0 

0 

0.0000 

0.0000 

0.0000 

0 

0 

0 

1.0000 

0.0000 

0.0000 

0 

0 

0 

0.0000 

-.8329 

.5534 

0 

0 

0 

0.0000 

0.0000 

0.0000 

0 

0 

0 

1.0000 

0.0000 

0.0000 

0 

0 

0 

0.0000 

-.8329 

.5534 

0 

0 

0 

0.0000 

0.0000 

0.0000 

0 

0 

0 

1.0000 

0.0000 

0.0000 

0 

0 

0 

0.0000 

-.8329 

.5534 

0 

0 

0 

0.0000 

0.0000 

0.0000 

0 

0 

0 

1.0000 

0.0000 

0.0000 

0 

0 

0 

0.0000 

-.8329 

.5534 

0 

0 

0 

0.0000 

0.0000 

0.0000 

0 

0 

0 

1.0000 

0.0000 

0.0000 

0 

0 

0 

0.0000 

-.8329 

.5534 

0 

0 

0 

0.0000 

0.0000 

0.0000 

0 

0 

0 

1.0000 

0.0000 

0.0000 

0 

0 

0 

0.0000 

-.8329 

.5534 

0 

0 

0 

0.0000 

0.0000 

0.0000 

0 

0 

0 

1.0000 

0.0000 

0.0000 

0 

0 

0 

0.0000 

-.8329 

.5534 

0 

0 

0 

0.0000 

0.0000 

0.0000 

0 

0 

0 

1.0000 

0.0000 

0.0000 

0 

0 

0 

0.0000 

-.8329 

.5534 

0 

0 

0 

0.0000 

0.0000 

0.0000 

0 

0 

0 

1.0000 

0.0000 

0.0000 

0 

0 

0 

0.0000 

-.8329 

.5534 

0 

0 

0 

0.0000 

0.0000 

0.0000 

0 

0 

0 

1.0000 

0.0000 

0.0000 

0 

0 

0 

0.0000 

-.8329 

.5534 

0 

0 

0 

0.0000 

0.0000 

0.0000 

0 

0 

0 

1.0000 

0.0000 

0.0000 

0 

0 

0 

0.0000 

-.8329 

.5534 

0 

0 

0 

0.0000 

0.0000 

0.0000 

0 

0 

0 

1.0000 

0.0000 

0.0000 

0 

0 

0 

0.0000 

-.8329 

.5534 

0 

0 

0 

0.0000 

0.0000 

0.0000 

0 

0 

0 

1.0000 

0.0000 

0.0000 

0 

0 

0 

0.0000 

-.8329 

.5534 

0 

0 

0 

0.0000 

0.0000 

0.0000 

0 

0 

0 

1.0000 

0.0000 

0.0000 

0 

0 

0 

0.0000 

-.8329 

.5534 

0 

0 

0 

0.0000 

0.0000 

0.0000 

0 

0 

0 

1.0000 

0.0000 

0.0000 

0 

0 

0 

0.0000 

-.8329 

.5534 

0 

0 

0 

0.0000 

0.0000 

0.0000 

0 

0 

0 

1.0000 

0.0000 

0.0000 

0 

0 

0 

0.0000 

-.8329 

.5534 

0 

0 

2  3  0 

122ADECC 

320 

ADMPC 

♦ 

♦ 

320 

320  3790000 

♦ 

377 — 

79 


0000 


0000. 


0000. 


3790000. 

2  .67000 

20000000.  20000000.  20000000.  .290  .290  .290  7750000.  775 


7750000. 


3 

♦ : 

10000 

10000000. 

10000000* 

3790000. 

1 

7201 

8201 

9201 

2 

7202 

8202 

9202 

3 

7203 

8203 

9203 

4 

7204 

8204 

9204 

5 

7205 

8205 

9205 

6 

7206 

8206 

9206 

7 

7207 

8207 

9207 

8 

7226 

8226 

9226 

9 

7227 

8227 

9227 

10 

7228 

8228 

9228 

11 

7229 

8229 

9229 

12 

7230 

8230 

9230 

13 

7231 

8231 

9231 

14 

7232 

8232 

9232 

15 

7251 

8251 

9251 

16 

7252 

8252 

9252 

17 

7253 

8253 

9253 

18 

7254 

8254 

9254 

19 

7255 

8255 

9255 

20 

7276 

8276 

9276 

21 

7277 

8277 

9277 

22 

7278 

8278 

9278 

23 

7279 

8279 

9279 

24 

7280 

8280 

9280 

25 

7301 

8301 

9301 

26 

7302 

8302 

9302 

27 

7326 

8326 

9326 

28 

7327 

8327 

9327 

29 

7351 

8351 

9351 

30 

7352 

8352 

9352 

31 

7201 

8201 

9201 

32 

7402 

8402 

9402 

33 

7403 

8403 

9403 

34 

7404 

8404 

9404 

35 

7405 

8405 

9405 

36 

7406 

8406 

9406 

37 

7407 

8407 

9407 

38 

7226 

8226 

9226 

39 

7427 

8427 

9427 

40 

7428 

8428 

9428 

41 

7429 

8429 

9429 

42 

7430 

8430 

9430 

43 

7431 

8431 

9431 

44 

7432 

8432 

9432 

45 

7251 

8251 

9251 

46 

7452 

8452 

9452 

47 

7453 

8453 

9453 

48 

7454 

8454 

9454 

49 

7455 

8455 

9455 

50 

7276 

8276 

9276 

10000000. 


.320  .320 


.320  3790000.  379 


ADMAX 


80 


51  7477  8477  9477 

52  7478  8478  9478 

53  7479  8479  9479 

54  7480  8480  9480 

55  7301  8301  9301 

56  7502  8502  9502 

57  7326  8326  9326 

58  7527  8527  9527 

59  7351  8351  9351 

60  7552  8552  9552 

61  7201  8201  9201 

62  7602  8602  9602 

63  7603  8603  9603 

64  7604  8604  9604 

65  7605  8605  9605 

66  7606  8606  9606 

67  7607  8607  9607 

68  7226  8226  9226 

69  7627  8627  9627 

70  7628  8628  9628 

71  7629  8629  9629 

72  7630  8630  9630 

73  7631  8631  9631 

74  7632  8632  9632 

75  7251  8251  9251 

76  7652  8652  9652 

77  7653  8653  9653 

78  7654  8654  9654 

79  7655  8655  9655 

80  7276  8276  9276 

81  7677  8677  9677 

82  7678  8678  9678 

83  7679  8679  9679 

84  7680  8680  9680 

85  7301  8301  9301 

86  7702  8702  9702 

87  7326  8326  9326 

88  7727  8727  9727 

89  7351  8351  9351 

90  7752  8752  9752 

91  7201  8201  9201 

92  7802  8802  9802 

93  7803  8803  9803 

94  7804  8804  9804 

95  7805  8805  9805 

96  7806  8806  9806 

97  7807  8807  9807 

98  7226  8226  9226 

99  7827  8827  9827 

100  7828  8828  9828 

101  7829  8829  9829 

102  7830  8830  9830 

103  7831  8831  9831 

104  7832  8832  9832 

105  7251  8251  9251 

106  7852  8852  9852 

107  7853  8853  9853 

108  7854  8854  9854 

109  7855  8855  9855 

110  7856  8856  9856 


81 


8276 

8877 

8878 

8879 

8880 

8881 

8301 

8902 

8326 

8927 

8351 

8952 

8 

402 

8 

403 

8 

404 

8 

405 

8 

406 

8 

407 

3 

408 

8 

427 

8 

428 

8 

429 

8 

430 

8 

431 

8 

432 

8 

433 

8 

452 

8 

453 

8 

454 

8 

455 

8 

456 

8 

477 

8 

478 

8 

479 

8 

480 

8 

481 


9276 

9877 

9878 

9879 

9880 

9881 
9301 
9902 
9326 
9927 
9351 
9952 

02  I-  0  0  ADELD 


401 

226 

227 

427 

426 

0 

O 

2 

0 

0 

402 

227 

228 

428 

427 

0 

1 

3 

0 

0 

403 

228 

229 

429 

428 

0 

1 

4 

0 

0 

404 

229 

230 

430 

429 

0 

1 

5 

0 

0 

405 

230 

231 

431 

430 

0 

1 

6 

0 

0 

406 

231 

232 

432 

431 

0 

1 

7 

0 

0 

407 

232 

233 

433 

432 

0 

2 

8 

0 

0 

426 

251 

252 

452 

451 

0 

2 

9 

0 

0 

427 

252 

253 

453 

452 

0 

1 

10 

0 

0 

428 

253 

254 

454 

453 

0 

1 

11 

0 

0 

429 

254 

255 

455 

454 

0 

1 

12 

0 

0 

430 

255 

256 

456 

455 

0 

1 

13 

0 

0 

431 

256 

257 

457 

456 

0 

1 

14 

0 

0 

432 

257 

258 

458 

457 

0 

2 

15 

0 

0 

451 

276 

277 

477 

476 

0 

2 

16 

0 

0 

452 

277 

278 

478 

477 

0 

1 

17 

0 

0 

453 

278 

279 

479 

478 

0 

1 

18 

0 

0 

454 

279 

280 

480 

479 

0 

1 

19 

0 

0 

455 

280 

281 

481 

480 

0 

2 

20 

0 

0 

476 

301 

302 

502 

501 

0 

2 

21 

0 

0 

477 

302 

303 

503 

502 

0 

1 

22 

0 

0 

478 

303 

304 

504 

503 

0 

1 

23 

0 

0 

479 

304 

305 

505 

504 

0 

1 

24 

0 

0 

480 

305 

306 

506 

505 

82 


25 

8 

8 

0 

2 

25 

0 

0 

301 

302 

502 

501 

326 

327 

527 

526 

26 

8 

8 

0 

2 

26 

0 

0 

302 

303 

503 

502 

327 

328 

528 

527 

27 

8 

8 

0 

2 

27 

0 

0 

326 

327 

527 

526 

351 

352 

552 

551 

28 

8 

8 

0 

2 

28 

0 

0 

327 

328 

528 

527 

352 

353 

553 

552 

29 

8 

8 

0 

1 

29 

0 

0 

351 

352 

552 

551 

376 

377 

577 

576 

30 

8 

8 

0 

1 

30 

0 

0 

352 

353 

553 

552 

377 

378 

578 

577 

31 

8 

8 

0 

2 

31 

0 

0 

201 

402 

602 

601 

226 

427 

627 

626 

32 

8 

8 

0 

2 

32 

0 

0 

402 

403 

603 

602 

427 

428 

628 

627 

33 

8 

8 

0 

1 

33 

0 

0 

403 

404 

604 

603 

428 

429 

629 

628 

34 

8 

8 

0 

1 

34 

0 

0 

404 

405 

605 

604 

429 

430 

630 

629 

35 

8 

8 

0 

1 

35 

0 

0 

405 

406 

606 

605 

430 

431 

631 

630 

36 

8 

8 

0 

1 

36 

0 

0 

406 

407 

607 

606 

431 

432 

632 

631 

37 

8 

8 

0 

1 

37 

0 

0 

407 

408 

608 

607 

432 

433 

633 

632 

38 

8 

8 

0 

2 

38 

0 

0 

226 

427 

627 

626 

251 

452 

652 

651 

39 

8 

8 

0 

2 

39 

0 

0 

427 

428 

628 

627 

452 

453 

653 

652 

40 

8 

8 

0 

1 

40 

0 

0 

428 

429 

629 

628 

453 

454 

654 

653 

41 

8 

8 

0 

1 

41 

0 

0 

429 

430 

630 

629 

454 

455 

655 

654 

42 

8 

8 

0 

1 

42 

0 

0 

430 

431 

631 

630 

455 

456 

656 

655 

43 

8 

8 

0 

1 

43 

0 

0 

431 

432 

632 

631 

456 

457 

657 

656 

44 

8 

8 

0 

1 

44 

0 

0 

432 

433 

633 

632 

457 

458 

658 

657 

45 

8 

8 

0 

2 

45 

0 

0 

251 

452 

652 

651 

276 

477 

677 

676 

46 

8 

8 

0 

2 

46 

0 

0 

452 

453 

653 

652 

477 

478 

678 

677 

47 

8 

8 

0 

1 

47 

0 

0 

453 

454 

654 

653 

478 

479 

679 

678 

48 

8 

8 

0 

1 

48 

0 

0 

454 

455 

655 

654 

479 

480 

680 

679 

49 

8 

8 

0 

1 

49 

0 

0 

455 

456 

656 

655 

480 

481 

681 

680 

50 

8 

8 

0 

2 

5Q 

0 

0 

276 

477 

677 

676 

301 

5Q2 

702 

701 

51 

8 

8 

0 

2 

51 

0 

0 

477 

478 

678 

677 

502 

503 

703 

702 

52 

8 

8 

0 

1 

52 

0 

0 

478 

479 

679 

678 

503 

504 

704 

703 

53 

8 

8 

0 

1 

53 

0 

0 

479 

480 

680 

679 

504 

505 

705 

704 

54 

8 

8 

0 

1 

54 

0 

0 

480 

481 

681 

680 

505 

5Q6 

706 

705 

83 


55 

301 

56 

502 

57 

326 

58 

527 

59 

351 

60 

552 

61 

201 

62 

602 

63 

603 

64 

604 

65 

605 

66 

606 

67 

607 

68 

226 

69 

627 

70 

628 

71 

629 

72 

630 

73 

631 

74 

632 

75 

251 

76 

652 

77 

653 

78 

654 

79 

655 

80 

276 

81 

677 

82 

678 

83 

679 

84 

680 


8 

8 

0 

2 

55 

0 

0 

502 

702 

701 

326 

527 

727 

726 

8 

8 

0 

2 

56 

0 

0 

503 

703 

702 

527 

528 

728 

727 

8 

8 

0 

2 

57 

0 

0 

527 

727 

726 

351 

552 

752 

751 

8 

8 

0 

2 

58 

0 

0 

528 

728 

727 

552 

553 

753 

752 

8 

8 

0 

1 

59 

0 

0 

552 

752 

751 

376 

577 

777 

776 

8 

8 

0 

1 

60 

0 

0 

553 

753 

752 

577 

578 

778 

777 

8 

8 

0 

2 

61 

0 

0 

602 

802 

801 

226 

627 

827 

826 

8 

8 

0 

2 

62 

0 

0 

603 

803 

802 

627 

628 

828 

827 

8 

8 

0 

1 

63 

0 

0 

604 

804 

803 

628 

629 

829 

828 

8 

8 

0 

1 

64 

0 

0 

605 

805 

804 

629 

630 

830 

829 

8 

8 

0 

1 

65 

0 

0 

606 

806 

805 

630 

631 

831 

830 

8 

8 

0 

1 

66 

0 

0 

607 

807 

806 

631 

632 

832 

831 

8 

8 

0 

1 

67 

0 

0 

608 

808 

807 

632 

633 

833 

832 

8 

8 

0 

2 

68 

0 

0 

627 

827 

826 

251 

652 

852 

851 

8 

8 

0 

2 

69 

0 

0 

628 

828 

827 

652 

653 

853 

852 

8 

8 

0 

1 

70 

0 

0 

629 

829 

828 

653 

654 

354 

853 

8 

8 

0 

1 

71 

0 

0 

630 

830 

829 

654 

655 

855 

854 

8 

8 

0 

1 

72 

0 

0 

631 

831 

830 

655 

656 

856 

855 

8 

8 

0 

1 

73 

0 

0 

632 

832 

831 

656 

657 

857 

856 

8 

8 

0 

1 

74 

0 

0 

633 

833 

832 

657 

658 

858 

857 

8 

8 

0 

2 

75 

0 

0 

652 

852 

851 

276 

677 

877 

876 

8 

8 

0 

2 

76 

0 

0 

653 

853 

852 

677 

678 

878 

877 

8 

8 

0 

1 

77 

0 

0 

654 

854 

853 

678 

679 

879 

878 

8 

8 

0 

1 

78 

0 

0 

655 

855 

854 

679 

680 

880 

879 

8 

8 

0 

1 

79 

0 

0 

656 

856 

855 

680 

681 

881 

880 

8 

8 

0 

2 

80 

0 

0 

677 

877 

876 

301 

702 

902 

901 

8 

8 

0 

2 

81 

0 

0 

678 

878 

877 

702 

703 

903 

902 

8 

8 

0 

1 

82 

0 

0 

679 

879 

878 

703 

704 

904 

903 

8 

8 

0 

1 

83 

0 

0 

680 

880 

879 

704 

705 

905 

904 

8 

8 

0 

1 

84 

0 

0 

681 

881 

880 

705 

706 

906 

905 

84 


85 

8 

8 

0 

301 

702 

902 

901 

86 

8 

8 

0 

702 

703 

903 

902 

87 

8 

8 

0 

326 

727 

927 

926 

88 

8 

8 

0 

727 

728 

928 

927 

89 

8 

8 

0 

351 

752 

952 

951 

90 

8 

8 

0 

752 

753 

953 

952 

91 

8 

8 

0 

201 

802 

1002 

1001 

92 

8 

8 

0 

802 

803 

1003 

1002 

93 

8 

8 

0 

803 

804 

1004 

1003 

94 

8 

8 

0 

804 

805 

1005 

1004 

95 

8 

8 

0 

805 

806 

1006 

1005 

96 

8 

8 

0 

806 

807 

1007 

1006 

97 

8 

8 

0 

807 

808 

1008 

1007 

98 

8 

8 

0 

226 

827 

1027 

1026 

99 

8 

8 

0 

827 

828 

1028 

1027 

100 

8 

8 

0 

828 

829 

1029 

1028 

101 

8 

8 

0 

829 

830 

1030 

1029 

102 

8 

8 

0 

830 

831 

1031 

1030 

103 

8 

8 

0 

831 

832 

1032 

1031 

104 

8 

8 

0 

832 

833 

1033 

1032 

105 

8 

8 

0 

251 

852 

1052 

1051 

106 

8 

8 

0 

852 

853 

1053 

1052 

107 

8 

8 

0 

853 

854 

1054 

1053 

108 

8 

8 

0 

854 

855 

1055 

1054 

109 

8 

8 

0 

855 

856 

1056 

1055 

110 

8 

8 

0 

856 

857 

1057 

1056 

111 

8 

8 

0 

276 

877 

1077 

1076 

112 

8 

8 

0 

877 

878 

1078 

1077 

113 

8 

8 

0 

878 

879 

1079 

1078 

114 

8 

8 

0 

879 

880 

1080 

1079 

85 

0 

0 

727 

927 

926 

86 

0 

0 

728 

928 

927 

87 

0 

0 

752 

952 

951 

88 

0 

0 

753 

953 

952 

89 

0 

0 

777 

977 

976 

90 

0 

0 

778 

978 

977 

91 

0 

0 

827 

1027 

1026 

92 

0 

0 

828 

1028 

1027 

93 

0 

0 

829 

1029 

1028 

94 

0 

0 

830 

1030 

1029 

95 

0 

0 

831 

1031 

1030 

96 

0 

0 

832 

1032 

1031 

97 

0 

0 

833 

1033 

1032 

98 

0 

0 

852 

1052 

1051 

99 

0 

0 

853 

1053 

1052 

100 

0 

0 

854 

1054 

1053 

101 

0 

0 

855 

1055 

1054 

102 

0 

0 

856 

1056 

1055 

103 

0 

0 

857 

1057 

1056 

104 

0 

0 

858 

1058 

1057 

105 

0 

0 

877 

1077 

1076 

106 

0 

0 

878 

1070 

1077 

107 

0 

0 

879 

1079 

1078 

108 

0 

0 

880 

1080 

1079 

109 

0 

0 

881 

1081 

1080 

110 

0 

0 

882 

1082 

1081 

111 

0 

0 

902 

1102 

1101 

112 

0 

0 

903 

1103 

1102 

113 

0 

0 

904 

1104 

1103 

114 

0 

0 

905 

1105 

1104 

2 

326 

2 

727 

2 

351 

2 

752 

1 

376 

1 

777 

2 

226 

2 

827 

1 

828 

1 

829 

1 

830 

1 

831 

1 

832 

2 

251 

2 

852 

1 

853 

1 

854 

1 

855 

1 

856 

1 

857 

2 

276 

9 

877 

1 

878 

1 

879 

1 

880 

3 

881 

2 

301 

2 

902 

1 

903 

1 

904 


85 


115 

880 

116 

881 

117 

301 

118 

902 

119 

326 

120 

927 

121 

351 

122 

952 


8 

8 

0 

1 

115 

0 

0 

881 

1081 

1080 

905 

906 

1106 

1105 

8 

8 

0 

3 

116 

0 

0 

882 

1082 

1081 

906 

907 

1107 

1106 

8 

8 

0 

2 

117 

0 

0 

902 

1102 

1101 

326 

927 

1127 

1126 

8 

8 

0 

2 

118 

0 

0 

903 

1103 

1102 

927 

928 

1128 

1127 

8 

8 

0 

2 

119 

0 

0 

927 

1127 

1126 

351 

952 

1152 

1151 

8 

8 

0 

2 

120 

0 

0 

928 

1128 

1127 

952 

953 

1153 

1152 

8 

8 

0 

1 

121 

0 

0 

952 

1152 

1151 

376 

977 

1177 

1176 

8 

8 

0 

1 

122 

0 

0 

953 

1153 

1152 

977 

978 

1178 

1177 

86 


APPENDIX  D 


Listing  of  File  NAS, 
the  Input  Cards  for  NASTRAN  Program 


87 


GRID 

2003 

1 

1.2500 

GRID 

4003 

1 

1.2500 

GRID 

2004 

1 

1.5000 

GRID 

4004 

1 

1.5000 

GRID 

2005 

1 

1.7500 

GRID 

4005 

1 

1.7500 

GRID 

2006 

1 

2.0000 

GRID 

4006 

1 

2.0000 

GRID 

2028 

1 

1.2500 

GRID 

4028 

1 

1.2500 

GRID 

2029 

1 

1.5000 

GRID 

4029 

1 

1.5000 

GRID 

2030 

1 

1.7500 

GRID 

4030 

1 

1.7500 

GRID 

2031 

1 

2.0000 

GRID 

403.1. 

1 

2.0000 

GRID 

2053 

1 

1.2500 

GRID 

4053 

1 

1.2500 

GRID 

2054 

1 

1.5000 

GRID 

4054 

1 

1.5000 

GRID 

2055 

1 

1 .7500 

GRID 

4055 

1 

1.7500 

GRID 

2056 

1 

2.0000 

GRID 

4056 

1 

2.0000 

GRID 

2077 

1 

1.0000 

GRID 

4077 

1 

1.0000 

GRID 

2078 

1 

1.2500 

GRID 

4078 

1 

1.2500 

GRID 

2079 

1 

1.5000 

GRID 

4079 

1 

1.5000 

GRID 

2080 

1 

1.7500 

GRID 

4080 

1 

1.7500 

GRID 

2081 

1 

2.0000 

GRID 

4081 

1 

2.0000 

GRID 

2102 

1 

1.0000 

GRID 

4102 

1 

1.0000 

GRID 

2103 

1 

1.2500 

GRID 

4103 

1 

1.2500 

GRID 

2104 

1 

1.5000 

GRID 

4104 

1 

1.5000 

GRID 

2.105 

1 

1.7500 

GRID 

4105 

1 

1.7500 

GRID 

2106 

1 

2.0000 

GRID 

4106 

1 

2.0000 

GRID 

2127 

1 

1.0000 

GRID 

4.127 

1 

1.0000 

GRID 

2128 

1 

1.2500 

GRID 

4128 

1 

1.2500 

GRID 

2129 

1 

1.5000 

GRID 

4 1 29 

1 

1.5000 

GRID 

2130 

1 

1.7500 

GRID 

4130 

1 

1.7500 

GRID 

2.131 

1 

2.0000 

GRID 

4131 

1 

2.0000 

GRID 

2152 

1 

1.0000 

GRID 

4152 

1 

1 .0000 

GRID 

2153 

1 

1.2500 

GRID 

4153 

1 

1.2500 

GRID 

2154 

1 

1.5000 

0.0000 

0.0000 

:i. 

2456 

12.0000 

0.0000 

1 

456 

0.0000 

0.0000 

;l. 

2456 

12.0000 

0.0000 

l 

456 

0.0000 

0.0000 

1 

2456 

12.0000 

0.0000 

1 

456 

0.0000 

0.0000 

i 

2456 

12.0000 

0.0000 

i 

456 

0.0000 

.5000 

i 

2456 

12.0000 

.5000 

i 

456 

0.0000 

.  5000 

1. 

2456 

12.0000 

.5000 

:L 

456 

0.0000 

.5000 

1 

2456 

12.0000 

.5000 

1 

456 

0.0000 

.  5000 

1. 

2456 

12.0000 

.5000 

:i. 

456 

0.0000 

1.0000 

:i. 

2456 

12.0000 

1.0000 

:i. 

456 

0.0000 

1.0000 

1 

2456 

12.0000 

1.0000 

i 

456 

0.0000 

1.0000 

i 

2456 

12.0000 

1.0000 

1. 

456 

0.0000 

1.0000 

1. 

2456 

12.0000 

1 . 0000 

:l. 

456 

0.0000 

1.5000 

1 

2456 

12.0000 

1.5000 

1 

456 

0.0000 

1 . 5000 

:1 

2456 

12.0000 

1.5000 

1 

456 

0.0000 

1.5000 

1 

2456 

12.0000 

1.5000 

1 

456 

0.0000 

1  .5000 

i 

2456 

12.0000 

1 .5000 

1 

456 

0.0000 

1.5000 

l 

2456 

12.0000 

1.5000 

1 

456 

0 . 0000 

2.0000 

1 

2456 

12.0000 

2.0000 

1 

456 

0.0000 

2.0000 

1 

2456 

12.0000 

2.0000 

l 

456 

0.0000 

2.0000 

1 

2456 

12.0000 

2.0000 

l 

456 

0.0000 

2.0000 

1 

2456 

12.0000 

2.0000 

1 

456 

0.0000 

2.0000 

1 

2456 

12.0000 

2.0000 

i 

456 

0.0000 

2.5000 

1 

2456 

12.0000 

2.5000 

1 

456 

0.0000 

2.5000 

1 

2456 

12.0000 

2.5000 

1 

456 

0.0000 

2.5000 

1 

2456 

12.0000 

2.5000 

i 

456 

0.0000 

2.5000 

l 

2456 

12.0000 

2.5000 

i 

456 

0.0000 

2.5000 

.1 

2456 

12.0000 

2.5000 

l 

456 

0.0000 

3.0000 

l 

2456 

12.0000 

3.0000 

l 

456 

0.0000 

3.0000 

l 

2456 

12.0000 

3.0000 

:l 

456 

0.0000 

3.0000 

1 

2456 

88 


GRID 

GRID 

GRID 

GRID 

GRID 

GRID 

GRID 

GRID 

GRID 

GRID 

GRID 

GRID 

GRID 

GRID 

GRID 

GRID 

GRID 

GRID 

GRID 

GRID 

GRID 

GRID 

GRID 

GRID 

GRID 

CHEXA 

2001 

+BC  2001 
CHEXA 

2002 

+BC  2002 
CHEXA 

2003 

+BC  2003 
CHEXA 

2004 

+BC  2004 
CHEXA 

2005 

+BC  2005 
CHEXA 

2006 

+BC  2006 
CHEXA 

2007 

+BC  2007 
CHEXA 

2008 

+BC  2008 
CHEXA 

2009 

+BC  2009 
CHEXA 

2010 

+BC  2010 
CHEXA 

2011 

+BC  2011 
CHEXA 

2012 


4154 

2155 

4155 

2156 

4156 

2177 

4177 

2178 

4178 

2179 

4179 

2180 

4180 
2181 

4181 
2202 

4202 

2203 

4203 

2204 

4204 

2205 

4205 

2206 
4  206 
2001 


4029 

2002 


4030 

2003 


4031 

2004 


4054 

2005 

4055 

2006 

4056 

2007 

4079 

2008 


4080 

2009 

4081 

2010 


4103 

2011 


4104 

2012 


1 

1  *  5000 

12.0000 

3.0000 

1 

1  *  7500 

0.0000 

3.0000 

1 

1  *  7500 

12.0000 

3.0000 

1 

2  .0000 

0.0000 

3.0000 

i 

2 .0000 

12.0000 

3.0000 

i 

1*0000 

0.0000 

3.5000 

1 

1*0000 

12.0000 

3.5000 

1 

1*2500 

0.0000 

3 .5000 

1 

1*2500 

12.0000 

3.5000 

1 

1*5000 

0.0000 

3.5000 

1 

1*5000 

12.0000 

3.5000 

i 

1*7500 

0.0000 

3.5000 

l 

1*7500 

12.0000 

3.5000 

l 

2.0000 

0.0000 

3.5000 

l 

2*0000 

12.0000 

3.5000 

l 

1.0000 

0.0000 

4.0000 

l 

1.0000 

12.0000 

4.0000 

l 

1.2500 

0.0000 

4.0000 

l 

1*2500 

12.0000 

4.0000 

l 

1.5000 

0.0000 

4.0000 

1 

1.5000 

12.0000 

4.0000 

l 

1.7500 

0.0000 

4.0000 

l 

1.7500 

12.0000 

4.0000 

l 

2.0000 

0.0000 

4.0000 

1 

2.0000 

12.0000 

4.0000 

l 

2003 

2004 

4004 

4028 

1 

2004 

2005 

4005 

4029 

1 

2005 

2006 

4006 

4030 

1 

2028 

2029 

4029 

4053 

1 

2029 

2030 

4030 

4054 

1 

2030 

2031 

4031 

4055 

1 

2053 

2054 

4054 

4078 

1 

2054 

2055 

4055 

4079 

i 

2055 

2056 

4056 

4080 

1 

2077 

2078 

4078 

4102 

1 

2078 

2079 

4079 

4103 

1 

2079 

2080 

4080 

1  456 

1  2456 

1  456 

1  2456 

1  456 

1  2456 

1  456 

1  2456 

1  456 

1  2456 

1  456 

1  2456 

1  456 

1  2456 

1  456 

1  23456 

1  3456 

1  2456 

1  456 

1  2456 

1  456 

1  2456 

1  456 

1  2456 

1  456 

4003  2028  2029CBC 


4004 

2029 

2030CBC 

4005 

2030 

2031 CBC 

4028 

2053 

2054CBC 

4029 

2054 

2055CBC 

4030 

2055 

2056CBC 

4053 

2078 

2079CBC 

4054 

2079 

2080CBC 

4055 

2080 

2081 CBC 

4077 

2102 

2103CBC 

4078 

2103 

2104CBC 

4079 

2104 

2105CBC 

89 


+BC  203.2 
CHEXA 

2013 

4105 

2013 

4104 

1 

2080 

2081 

4081 

4080 

2105 

2106CBC 

iBC  2013 
CHEXA 

2014 

4106 

2014 

4105 

1 

2102 

2103 

4103 

4102 

2127 

2128CBC 

+  BC  2014 
CHEXA 

2015 

4128 

2015 

4127 

1 

2103 

2104 

4104 

4103 

2128 

2129CBC 

+BC  2015 
CHEXA 

2016 

4129 

2016 

4128 

1 

2104 

2105 

4105 

4104 

2129 

2130CBC 

+BC  2016 
CHEXA 

2017 

4130 

2017 

412? 

1 

2105 

2106 

4106 

4105 

2130 

2131CBC 

+BC  2017 
CHEXA 

2018 

4131 

2018 

4130 

1 

2127 

2128 

4128 

4127 

2152 

2153CBC 

+BC  2018 
CHEXA 

2019 

4153 

2019 

4152 

1 

2128 

2129 

4129 

4128 

2153 

2154CBC 

+BC  2019 
CHEXA 

2020 

4154 

2020 

4153 

1 

2129 

2130 

4130 

4129 

2154 

2155CBC 

•fBC  2020 
CHEXA 

2021 

4155 

2021 

4154 

1 

2130 

2131 

4131 

4130 

2155 

2156CBC 

•fBC  2021 
CHEXA 

2022 

4156 

2022 

4155 

1 

2152 

2153 

4153 

4152 

2177 

2178CBC 

•fBC  2022 
CHEXA 

2023 

4178 

2023 

4177 

1 

2153 

2154 

4154 

4153 

2178 

2179CBC 

+BC  2023 
CHEXA 

2024 

4179 

2024 

4178 

1 

2154 

2155 

4155 

4154 

2179 

2180CBC 

+BC  2024 
CHEXA 

2025 

4180 

2025 

4179 

1 

n  i  kk 

4l •  A  w  w 

2156 

4156 

4155 

2180 

2181CBC 

•fBC  2025 
CHEXA 

2026 

4181 

2026 

4180 

1 

2177 

2178 

4178 

4177 

2202 

2203CBC 

+BC  2026 
CHEXA 

2027 

4203 

2027 

4202 

1 

2178 

2179 

4179 

4178 

2203 

2204CBC 

+BC  2027 
CHEXA 

2028 

4204 

2028 

4203 

1 

2179 

2180 

4180 

4179 

2204 

2205CBC 

+BC  2028 
CHEXA 

2029 

4205 

2029 

4204 

1 

2180 

2181 

4181 

4180 

2205 

2206CBC 

+BC  2029 
PLOAD 

PLOAD 

PLOAD 

PLOAD 

PLOAD 

PLOAD 

PLOAD 

PLOAD 

4206  4205 
130000  *00 

1 30000  *  00 
130000*00 
130000*00 
130000*00 
130000*00 
130000.00 
130000.00 

2077 

2078 
2053 
2028 

2003 

2004 

2005 

2006 

2078 

2053 

2028 

2003 

2004 

2005 

2006 
2031 

4078 

4053 

4028 

4003 

4004 

4005 

4006 
4031 

4077 

4078 
4053 
4028 

4003 

4004 

4005 

4006 

90 


F'LOAD 

130000 

♦  00 

2031 

2056 

4056 

4031 

PLOAD 

130000 

.00 

2056 

2081 

4031 

4056 

GF:  I  Ei 

10003 

1 

1.2500 

78.0000 

0.0000 

1 

456 

GRID 

12003 

1 

1.2500 

90.0000 

0.0000 

1 

2456 

GRID 

10004 

1 

1.5000 

78.0000 

0.0000 

1 

456 

GR I D 

12004 

1 

1.5000 

90.0000 

0.0000 

1 

2456 

GRID 

10005 

i 

1.7500 

78.0000 

0.0000 

1 

456 

GRID 

12005 

1 

1.7500 

90.0000 

0.0000 

1 

2456 

GRID 

10006 

1 

2.0000 

78.0000 

0.0000 

1 

456 

GRID 

12006 

1 

2.0000 

90.0000 

0.0000 

1 

2456 

GRID 

10028 

1 

1.2500 

78.0000 

.5000 

1 

456 

GRID 

12028 

1 

1.2500 

90.0000 

♦  5000 

1 

2456 

GRID 

10029 

1 

1.5000 

78.0000 

♦  5000 

1 

456 

GRID 

12029 

1 

1.5000 

90,0000 

.5000 

1 

2456 

GRID 

10030 

1 

1.7500 

78.0000 

.5000 

1 

456 

GRID 

12030 

1 

1.7500 

90.0000 

.5000 

I 

2456 

GRID 

10031 

1 

2.0000 

78.0000 

.5000 

1 

456 

GRID 

12031 

1 

2.0000 

90,0000 

.5000 

1 

2456 

GRID 

10053 

1 

1.2500 

78.0000 

1.0000 

1 

456 

GRID 

12053 

1 

1.2500 

90,0000 

1.0000 

1 

2456 

GRID 

10054 

1 

1.5000 

78.0000 

1.0000 

1 

456 

GRID 

12054 

1 

1.5000 

90.0000 

1.0000 

1 

2456 

GRID 

10055 

1 

1.7500 

78.0000 

1.0000 

1 

456 

GRID 

12055 

1 

1.7500 

90.0000 

1,0000 

1 

2456 

GRID 

10056 

1 

2.0000 

78.0000 

1 .0000 

1 

456 

GRID 

12056 

1 

2.0000 

90 » 0000 

1 .0000 

1 

2456 

GF;  IEi 

10077 

1 

1.0000 

78.0000 

1 .5000 

1 

456 

GRID 

12077 

1 

1.0000 

90.0000 

1 ,5000 

1 

2456 

GRID 

10078 

1 

1.2500 

78.0000 

1 .5000 

1 

456 

GRID 

12078 

1 

1.2500 

90.0000 

1 .5000 

1 

2456 

GRID 

10079 

1 

1.5000 

78.0000 

1  .5000 

1 

456 

GRID 

12079 

1 

1.5000 

90.0000 

1.5000 

1 

2456 

GRID 

10080 

1 

1.7500 

78.0000 

1.5000 

1 

456 

GRID 

12080 

1 

1.7500 

90.0000 

1 .5000 

1 

2456 

GRID 

10081 

1 

2.0000 

78.0000 

1.5000 

1 

456 

GRID 

12081 

1 

2.0000 

90 . 0000 

1.5000 

1 

2456 

GRID 

10102 

1 

1.0000 

78.0000 

2.0000 

1 

456 

GRID 

12102 

1 

1.0000 

90.0000 

2.0000 

1 

2456 

G  R I  Ei 

10103 

1 

1.2500 

78.0000 

2.0000 

1 

456 

GRID 

12103 

1 

1.2500 

90.0000 

2.0000 

1 

2456 

GRID 

10104 

1 

1.5000 

78.0000 

2.0000 

1 

456 

GRID 

12104 

1 

1.5000 

90.0000 

2.0000 

1 

2456 

GRID 

10105 

1 

1.7500 

78.0000 

2.0000 

1 

456 

GRID 

12105 

1 

1.7500 

90.0000 

2.0000 

1 

2456 

GR  I  Ei 

10106 

1 

2.0000 

78.0000 

2.0000 

1 

456 

GRID 

12106 

1 

2.0000 

90.0000 

2.0000 

I 

2456 

GRID 

10127 

1 

1.0000 

78.0000 

2.5000 

1 

456 

GRID 

12127 

1 

1.0000 

90.0000 

2.5000 

1 

2456 

GRID 

10128 

1 

1.2500 

78.0000 

2.5000 

1 

456 

G  R I  Ei 

12128 

1 

1.2500 

90.0000 

2.5000 

1 

2456 

GRID 

10129 

1 

1.5000 

78.0000 

2.5000 

1 

456 

GRID 

12129 

1 

1.5000 

90.0000 

2,5000 

1 

2456 

GRID 

10130 

1 

1.7500 

78.0000 

2.5000 

1 

456 

GRID 

12130 

1 

1.7500 

90.0000 

2.5000 

1 

2456 

GRID 

10131 

1 

2.0000 

78.0000 

2.5000 

1 

456 

GRID 

12131 

1 

2.0000 

90.0000 

2.5000 

1 

2456 

GR  I  Ei 

10152 

1 

1.0000 

78.0000 

3.0000 

I 

456 

GRID 

12152 

1 

1.0000 

90.0000 

3.0000 

1 

2456 

GRID 

10153 

3. 

1.2500 

78.0000 

3.0000 

1 

456 

GRID 

12153 

1 

1.2500 

90.0000 

3.0000 

1 

2456 
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GRID 

10154 

GRID 

12154 

GRID 

10155 

GRID 

12155 

GRID 

10156 

GRID 

12156 

GRID 

10177 

GRID 

12177 

GRID 

10178 

GRID 

12178 

GRID 

10179 

GRID 

12179 

GRID 

10180 

GRID 

12180 

GRID 

10181 

GRID 

12181 

GRID 

10202 

GRID 

12202 

GRID 

10203 

GRID 

12203 

GRID 

10204 

GRID 

12204 

GRID 

10205 

GRID 

12205 

GRID 

10206 

GRID 

12206 

CHEXA 

1000.1 

10001 


1  1*5000  78*0000  3*0000 

1  1*5000  90*0000  3.0000 

1  1.7500  78*0000  3.0000 

1  1*7500  90.0000  3.0000 

1  2.0000  78.0000  3.0000 

1  2.0000  90.0000  3.0000 

1  1.0000  78.0000  3.5000 

1  1.0000  90.0000  3.5000 

1  1.2500  78.0000  3.5000 

1  1.2500  90.0000  3.5000 

1  1.5000  78.0000  3.5000 

1  1.5000  90.0000  3.5000 

1  1.7500  78.0000  3.5000 

1  1.7500  90.0000  3.5000 

1  2.0000  78.0000  3.5000 

1  2.0000  90.0000  3.5000 

1  1.0000  78.0000  4.0000 

1  1.0000  90.0000  4.0000 

1  1.2500  78.0000  4.0000 

1  1.2500  90.0000  4.0000 

1  1.5000  78.0000  4.0000 

1  1.5000  90.0000  4.0000 

1  1.7500  78.0000  4.0000 

1  1.7500  90.0000  4.0000 

1  2.0000  78.0000  4.0000 

1  2.0000  90.0000  4.0000 

1  10003  10004  12004 


1  456 

1  2456 

1  456 

1  2456 

1  456 

1  2456 

1  456 

1  2456 

1  456 

1  2456 

1  456 

1  2456 

1  456 

1  2456 

1  456 

1  2456 

1  3456 

1  23456 

1  456 

1  2456 

1  456 

1  2456 

1  456 

1  2456 

1  456 

1  2456 

12003  10028  1 0029CBC 


+BC10001 

.12029 

12028 

CHEXA 

10002 

10002 

1 

10004 

10005 

12005 

.12004 

10029 

10030CBC 

+BC10002 

12030 

12029 

CHEXA 

10003 

10003 

1 

10005 

10006 

12006 

12005 

10030 

10031CBC 

+BC10003 

12031 

12030 

CHEXA 

10004 

10004 

1 

10028 

10029 

12029 

12028 

10053 

10054CBC 

+BC10004 

12054 

12053 

CHEXA 

10005 

10005 

1 

10029 

10030 

12030 

12029 

10054 

1 0055CBC 

+BC10005 

12055 

12054 

CHEXA 

10006 

10006 

1 

10030 

10031 

12031 

12030 

10055 

10056CBC 

+BC10006 

12056 

12055 

CHEXA 

10007 

10007 

1 

10053 

10054 

12054 

12053 

10078 

10079CBC 

+BC10007 

12079 

12078 

CHEXA 

10008 

10008 

1 

10054 

10055 

12055 

12054 

10079 

10080CBC 

+BC10008 

12080 

12079 

CHEXA 

10009 

10009 

1 

10055 

10056 

12056 

12055 

10080 

1 0081.  CBC 

+BC10009 

12081 

12080 

CHEXA 

10010 

10010 

1 

10077 

10078 

12078 

12077 

10102 

10103CBC 

+BC10010 

12103 

12102 

CHEXA 

10011 

10011 

1 

10078 

10079 

12079 

12078 

10103 

1 0104CBC 

+BC1001 1 

12104 

12103 

CHEXA 

10012 

1 

10079 

10080 

12080 

12079 

10104 

10105CBC 
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10012 

+BC10012 

CHEXA 

10013 

12105 

10013 

12104 

1 

10080 

10081 

12081 

12080 

10105 

10106CBC 

+BC10013 

CHEXA 

10014 

12106 

10014 

12105 

1 

10102 

10103 

12103 

12102 

10127 

10128CBC 

+BC10014 

CHEXA 

10015 

12128 

10015 

12127 

1 

10103 

10104 

12104 

12103 

10128 

10129CPC 

+BC10015 

CHEXA 

10016 

12129 

10016 

12128 

1 

10104 

10105 

12105 

12104 

10129 

10130CBC 

+BC10016 

CHEXA 

10017 

12130 

10017 

12129 

1 

10105 

10106 

12106 

12105 

10130 

10131CBC 

+BC10017 

CHEXA 

10018 

12131 

10018 

12130 

1 

10127 

10128 

12128 

12127 

10152 

10153CBC 

+BC10018 

CHEXA 

10019 

12153 

10019 

12152 

1 

10128 

10129 

12129 

12128 

10153 

1 0154CBC 

+BC10019 

CHEXA 

10020 

12154 

10020 

12153 

1 

10129 

10130 

12130 

12129 

10154 

10155CBC 

+BC10020 

CHEXA 

10021 

12155 

10021 

12154 

1 

10130 

10131 

12131 

12130 

10155 

1 0156CBC 

+BC10021 

CHEXA 

10022 

12156 

10022 

12155 

1 

10152 

10153 

12153 

12152 

10177 

10178CBC 

+BC10022 

CHEXA 

10023 

12178 

10023 

12177 

1 

10153 

10154 

12154 

12153 

10178 

1 0179CBC 

+BC10023 

CHEXA 

10024 

12179 

10024 

12178 

1 

10154 

10155 

12155 

12154 

10179 

1 0180CBC 

+BC10024 

CHEXA 

10025 

12180 

10025 

12179 

1 

10155 

10156 

12156 

12155 

10180 

10181CBC 

+BC10025 

CHEXA 

10026 

12181 

10026 

12180 

1 

10177 

10178 

12178 

12177 

10202 

10203CBC 

+BC10026 

CHEXA 

10027 

12203 

10027 

12202 

1 

10178 

10179 

12179 

12178 

10203 

10204CBC 

+BC10027 

CHEXA 

10028 

12204 

10028 

12203 

1 

10179 

10180 

12180 

12179 

10204 

10205CBC 

+BC10028 

CHEXA 

10029 

12205 

10029 

12204 

1 

10180 

10181 

12181 

12180 

10205 

10206CBC 

+BC10029 

PLOAD 

PLOAD 

PLOAD 

PLOAD 

PLOAD 

12206  12205 
130000*00 
130000*00 
130000*00 
130000*00 
130000*00 

10077 

10078 
10053 
10028 
10003 

10078 

10053 

10028 

10003 

10004 

12078 

12053 

12028 

12003 

12004 

12077 

12078 
12053 
12028 
12003 

PLOAD 

PLOAD 

130000*00 

130000.00 

10004 

10005 

10005 

lOOOt* 

12005 

12006 

12004 

12005 
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PLOAD 

1 

PLOAD 

1 

PLOAD 

1 

GRID 

4002 

GRID 

6002 

GRID 

6003 

GRID 

6004 

GRID 

6005 

GRID 

6006 

GRID 

4027 

GRID 

6027 

GRID 

6028 

GRID 

6029 

GRID 

6030 

Gi"<  Hi 

6031 

GRID 

4052 

GRID 

6052 

GRID 

6053 

GRID 

6054 

GRID 

6055 

GRID 

6056 

GRID 

6077 

GRID 

6078 

GRID 

6079 

GRID 

6080 

GRID 

6081 

GRID 

6102 

GRID 

6103 

GRID 

6104 

GRID 

6105 

GRID 

6106 

GRID 

6127 

GRID 

6128 

GRID 

6129 

GRID 

6130 

GRID 

6131 

GRID 

6152 

GRID 

6153 

GRID 

6154 

GRID 

6155 

GRID 

6156 

GRID 

6177 

GRID 

6178 

GRID 

6179 

GRID 

6180 

GRID 

6181 

GRID 

6202 

GRID 

6203 

GRID 

6204 

GRID 

6205 

GRID 

6206 

CHEXA 

4001 

4001 

+BC  4001 

6028 

CHEXA 

4002 

4002 

+BC  4002 

6029 

CHEXA 

4003 

4003 

+BC  4003 

6030 

30000  *  00 

10006 

1 0  -j  3 1 

30000  *  00 

10031 

10056 

30000  *  00 

10056 

10081 

1 

1  . 0000 

12.0000 

i 

1.0000 

34.0000 

1 

1.2500 

34.0000 

i 

1.5000 

34.0000 

l 

1.7500 

34.0000 

i 

2.0000 

34.0000 

i 

1.0000 

12.0000 

1 

1.0000 

34.0000 

i 

1.2500 

34.0000 

i 

1.5000 

34.0000 

i 

1.7500 

34 .0000 

i 

2.0000 

34.0000 

l 

1.0000 

12.0000 

1 

1.0000 

34.0000 

1 

1.2500 

34.0000 

1 

1.5000 

34.0000 

1 

1.7500 

34.0000 

i 

2.0000 

34.0000 

1 

1.0000 

34.0000 

1 

1.2500 

34.0000 

1 

1.5000 

34.0000 

1 

1.7500 

34 .0000 

1 

2.0000 

34.0000 

l 

1.0000 

34.0000 

1 

1.2500 

34.0000 

1 

1.5000 

34.0000 

i 

1.7500 

34.0000 

l 

2.0000 

34.0000 

l 

1.0000 

34.0000 

i 

1.2500 

34.0000 

1 

1.5000  : 

34.0000 

1 

1.7500 

34.0000 

1 

2.0000 

34.0000 

1 

1.0000 

34.0000 

i 

1.2500 

34.0000 

1 

1.5000 

34.0000 

1 

1.7500 

34.0000 

1 

2.0000 

34.0000 

1 

1.0000 

34.0000 

1 

1.2500 

34.0000 

1 

1.5000 

34.0000 

1 

1.7500 

34.0000 

1 

2.0000 

34.0000 

1 

1.0000 

34.0000 

1 

1.2500 

34.0000 

1 

1.5000 

34.0000 

1 

1.7500 

34.0000 

1 

2.0000 

34.0000 

i 

4002 

4003 

6027 

1 

4003 

4004 

6028 

1 

4004 

4005 

6029 


12031 

12056 

12081 

0.0000 

12006 

12031 

12056 

1 

456 

0.0000 

1 

456 

0.0000 

1 

456 

0.0000 

1 

456 

0.0000 

1 

456 

0,0000 

1 

456 

.5000 

1 

456 

.5000 

1 

456 

.5000 

1 

456 

.5000 

1 

456 

.5000 

1 

456 

.5000 

1 

456 

1.0000 

1 

456 

1.0000 

1 

456 

1.0000 

1 

456 

1.0000 

1 

456 

1.0000 

1 

456 

1.0000 

1 

456 

1.5000 

1 

456 

1.5000 

1 

456 

1.5000 

1 

456 

1.5000 

1 

456 

1.5000 

1 

456 

2.0000 

1 

456 

2.0000 

1 

456 

2.0000 

1 

456 

2.0000 

1 

456 

2.0000 

1 

456 

2.5000 

1 

456 

2.5000 

1 

456 

2.5000 

1 

456 

2.5000 

1 

456 

2.5000 

1 

456 

3.0000 

1 

456 

3.0000 

1 

456 

3.0000 

1 

456 

3.0000 

1 

456 

3.0000 

1 

456 

3.5000 

1 

456 

3.5000 

1 

456 

3.5000 

1 

456 

3.5000 

1 

456 

3.5000 

1 

456 

4.0000 

1 

3456 

4.0000 

1 

456 

4.0000 

1 

456 

4.0000 

1 

456 

4.0000 

1 

456 

6003 

6002 

4027 

6004 

6003 

4028 

6005 

6004 

402? 

4028CBC 

4029CBC 

4030CBC 


94 


CHEXA 

4004 

4004 

+BC  4004 

4031 

CHEXA 

4005 

4005 

+BC  4005 

4053 

CHEXA 

4004 

4006 

+BC  4004 

4054 

CHEXA 

4007 

4007 

+BC  4007 

4055 

CHEXA 

4008 

4008 

+BC  4008 

4054 

CHEXA 

4009 

4009 

+BC  4009 

4078 

CHEXA 

4010 

4010 

+BC  4010 

4079 

CHEXA 

4011 

4011 

+BC  4011 

4080 

CHEXA 

4012 

4012 

+BC  4012 

4081 

CHEXA 

4013 

4013 

+BC  4013 

4103 

CHEXA 

4014 

4014 

+BC  4014 

4104 

CHEXA 

4015 

4  015 

+BC  4015 

4105 

CHEXA 

4014 

4014 

+BC  4014 

4104 

CHEXA 

4017 

4017 

+BC  4017 

4128 

CHEXA 

4018 

4018 

+BC  4018 

4129 

CHEXA 

4019 

4019 

+BC  4019 

4130 

CHEXA 

4020 

4020 

+BC  4020 

4131 

CHEXA 

4021 

4021 

+BC  4021 

4153 

CHEXA 

4022 

4022 

+BC  4022 

4154 

CHEXA 

4023 

4023 

+BC  4023 

4155 

1 

4030 

4005 

4004 

4004 

1 

4052 

4027 

4028 

4028 

1 

4053 

4028 

4029 

4029 

1 

4054 

4029 

4030 

4030 

1 

4055 

4030 

4031 

6031 

1 

4077 

4052 

4053 

6053 

1 

4078 

4053 

4054 

6054 

1 

4079 

4054 

4055 

6055 

1 

4080 

4055 

4054 

6056 

1 

4102 

4077 

4078 

6078 

1 

4103 

4078 

4079 

6079 

1 

4104 

4079 

4080 

6080 

1 

4105 

4080 

4081 

6081 

1 

4127 

4102 

4103 

6103 

1 

4128 

4103 

4104 

6104 

1 

4129 

4104 

4105 

6105 

1 

4130 

4105 

4104 

6106 

1 

4152 

4127 

4128 

6128 

1 

4153 

4128 

4129 

6129 

1 

4154 

4129 

4130 

6130 

6005 

4030 

4031 CBC 

6027 

4052 

4053CBC 

6023 

4053 

4 054 CBC 

6029 

4054 

4055  CBC 

6030 

4055 

4056CBC 

6052 

4077 

4078CBC 

6053 

4078 

4 07 9 CBC 

6054 

4079 

4 080 CBC 

6055 

4080 

4081 CBC 

6077 

4102 

4103CBC 

6078 

4103 

4104CBC 

6079 

4104 

4105CBC 

6080 

4105 

4106CBC 

6102 

4127 

4128CBC 

6103 

4128 

4129CBC 

6104 

4129 

4130CBC 

6105 

4130 

4131 CBC 

6127 

4152 

4153CBC 

6128 

4153 

4154CBC 

6129 

4154 

4 155 CBC 
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CHEXA 

4024 

1 

4130 

4131 

6131 

6130 

4155 

4156CBC 

4024 

+BC  4024 

6156 

6155 

CHEXA 

4025 

1 

4152 

4153 

6153 

6152 

4177 

4178CBC 

4025 

+BC  4025 

6178 

6177 

CHEXA 

4026 

1 

4153 

4154 

6154 

6153 

4178 

4179CBC 

4026 

+  BC  4026 

6179 

6178 

CHEXA 

4027 

1 

4154 

4155 

6155 

6154 

4179 

4180CBC 

4027 

+BC  4027 

6180 

6179 

CHEXA 

4028 

1 

4155 

4156 

6156 

6155 

4180 

4181CBC 

4028 

+BC  4028 

6181 

6180 

CHEXA 

4029 

1 

4177 

4178 

6178 

6177 

4202 

4203CBC 

4029 

+BC  4029 

6203 

6202 

CHEXA 

4030 

4030 

1 

4178 

4179 

6179 

6178 

4203 

4204CBC 

+BC  4030 

6204 

6203 

CHEXA 

4031 

1 

4179 

4180 

6180 

6179 

4204 

4205CBC 

4031 

+BC  4031 

6205 

6204 

CHEXA 

4032 

1 

4180 

4181 

6181 

6180 

4205 

4206CBC 

4032 

+BC  4032 

6206 

6205 

PLOAD 

130000*00 

4002 

4003 

6003 

6002 

PLOAD 

130000.00 

4003 

4004 

6004 

6003 

PLOAD 

130000.00 

4004 

4005 

6005 

6004 

PLOAD 

130000.00 

4005 

4006 

6006 

6005 

PLOAD 

130000.00 

4006 

4031 

6031 

6006 

PLOAD 

130000.00 

4031 

4056 

6056 

6031 

PLOAD 

130000.00 

4056 

4081 

6081 

6056 

GRID 

8002 

1 

1.0000 

56.0000 

0.0000 

1 

456 

GRID 

8003 

1 

1.2500 

56.0000 

0.0000 

1 

456 

GRID 

8004 

1 

1.5000 

56.0000 

0.0000 

1 

456 

GRID 

8005 

I 

1.7500 

56.0000 
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